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I. Introduction 

This document represents the final technical report of research 
sponsored under NASA Grant NSG 5209, Multidisciplinary Research 
Program in Atmospheric Science. ?.*he present document consists of 
ten individual research reports involving studies of satellite and 
radar based indirect sensing of the atmosphere; aspects of atmoa- 
pheric radiative transfer through cloudy atmospheres; and numerical 
modeling cf large scale atmospheric flow. Combined with a previous 
report*, there are, thus, twenty-two extensive reviews of research 
conducted during the period of this grant. 

In addition to the individual research projects, the grant 
has sponsored joint NASA-U. Maryland Summer Research Colloquia in 
atmospheric science each summer since 1977, and a special working 
seminar on a tmospheric_ blocking during the past 15 months. The 
grant has also supported developmental research involving remote 
computer graphics capabilities. This has led to system software 
development which will support graphics on three separate types of 
plotting equipment for the dozen or so remote university user 
groups of the NASA-Goddard Modelling and Simulation Facility. 


♦Multidisciplinary Research Pro gra m in Atmospheric Science- 
A Quasi-Biennial Report. Department of Meteorology, University 
Maryland, College Park, MD. April, 1980. (O. E. Thompson, 

Principal Investigator.) 
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ABSTRACT 


A theoretical analysis of tha vortical resolving power of 
tho High resolution Infra-red Radiation Sounder CHIR8) and 
the Advanced Meteorological Temperature Sounder (AMTS) is 
carried out. First* the infra-red transmi t tanco weighting 
functions* and associated radiative transfer kernels* are 
analysed through singular value decomposition. The AMTS is 
found to contain several more pieces of independent 
information than H1R3 when the transmittances are 
considered* but the two instruments appear to be much more 
similar when the temperature sensitive radiative transfer 
kernels are analysed. 


The HIRS and AMTS Instruments ar 
thorough analysis using the theoret 
and Gilbert. From this analysis* it 
instruments should have very simll 
power below SOOmb but that AMTS 
resolving power above 200mb. In the 
AMT? shows badly degraded spread fun 
not be a realistic assessment of ve 
there. 


e also subjected to a 
i c a 1 methods of Backus 
is found that the two 
ar vertical resolving 
should have superior 
layer 200-500mb* the 
ction which may or may 
rtical resolving power 
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1. Introduction 


This paper concerns ths vortical resolving powar of 
satellite born* temperature sounding instruments. Specifi- 
cally* wa wish to prasant information on tha capabilities of 
tha High Resolution Infra-rad Radiation 8ounder (HIRS) and a 
proposed sounding instrument called the Advanced Moisture 
and Temperature Sounder (AMTS) designed following tha work 
of Kaplan* Chahine* Susskind and Saarla (1977). In this 
paper and tha one which follows* wa discuss two quite diffe- 
rent methods for assessing tha vertical resolving power of 
satellite sounders. The first is the theoretical method of 
Conrath (1972) which was patterned after the work of Backus 
and Gilbert (1968). In this method* vertical resolving power 
is defined in terms of certain properties of the transmit- 
tance weighting functions corresponding to the sounding 
channels of the satellite radiometric device. Resolving 
power defined in this theoretical approach desls with the 
c harac ter i st ic thickness of averaging kernels which relate a 
temperature estimate* at a given atmospheric level* to the 
distribution of temperature throughout height as it eiists 
to produce the radiation measurements. The relation between 
the two is an integral relation deriving from the atmos- 
pheric radiative transfer equation with a presumed linear 
dependence between estimated temperature and measured radi- 
ance. The Backus-Gi lbert-Conrath (BGC) approach includes a 
formalism for deriving a retrieval algorithm for optimizing 
the vertical resolving powor. However* a retrieval algor- 
ithm constructed in the BGC optimal fashion mag not neces- 
sarily be optimal as far as actual temperature retrievals 
are concerned. This is due to the fact that integral averag- 
ing kernels which have minimum thickness may also have unde- 
sirably large side-lobes which degrade the overall retrieval 
of a profile of temperature. 

For the reasons introduced above* we have developed an 
independent criterion for vertical resolving power which is 
based on actual retrievals of signal structure in the temp- 
erature field. (See the following paper by the present 
author. ) The criterion is patterned after that introduced by 
Lord Rayleigh for testing the resolving power of rptical 
instruments. Specifically* in the paper which follows* we 
define the vertical resolving length of a satellite based 
temperature sounding system as the minimum separation of two 
distinct signals in the "field of view" (that is* along the 
vertical profile) which can just j e r-solved as distinct 
signals by the observing system. 
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In what it to follow in this paper* wo present on analy- 
•it of tho vorticol revolving powor of the H1RS ond AMTS 
instrument* using tho thooroticol critoria dovolopod by Rec- 
tus ond Oilbort ond Conroth. Tho two otudioo toton together 
will provide o comparison of thooroticol ond ompiricol 
results concerning vorticol resolution of sotollito tempore- 
turo sounder*. It will else yi»!d o fairly consistent com- 
parison of tho vortical resolution of HIRS ond AMT8 instru- 
ments* but will also rovool some important discrepancies in 
the two methods of analysis. Since vertical resolving power 
is expected to be o strong function of instrument noise 
charec terist ics« we will carry out our analysis for the most 
representative estimates of instrument noise available for 
these very high technology instruments. 

2. The concept of vertical resolving power 

It is worth a few paragraphs to discuss the general con- 
cept of vertical resolving power of a satellite temperature 
sounding system. Even for radiosonde measurements of the 
temperature profile* it is not immediately clear what the 
vertical resolving power actually is. For a sensor of low 
thermal inertia (quick response to temperature change) the 
resolving power should be closely related to the "least 
count" of the system! that is* the vertical spacing between 
successive measurements. If the sensor has large thermal 
inertia so that its response time exceeds the time to ascend 
to the next measurement level* the resolving length would 
exceed the least count. If the sensor or its attendant 
electronic package were noisy* then the true resolving 
length becomes a statistical function of the instrument 
noise as well. Further* if the telemetering system is not 
perfect* then even the least count distance is imperfectly 
known and the true .esolving length also depends on the 
errors in the telemetering system. 

In a satellite temperature sounding system* the situa- 
tion is even more complex. The integral nature of the radia- 
tive transfer physics of the atmosphere assures that true 
temperature signals are grossly smoothed into bulk radiation 
measurements. (See* for example* Thompson* Eom and Wagen- 
hofer (1*776). ) While narrowing the spectral bandwidth of 
the radiometric instrument will sharpen the transmi ttance 
weighting functions* one can never design an instrument with 
a true Dirac delta-function relationship to the temperature 
profile. Thus* the integral smoothing is always real and 
significant and vertical resolution will always suffer 
because of it. The integral nature of atmospheric radiative 
transfer also makes the least count of tho system (in this 
case* the spacing between adjacently positioned 
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transmi ttance weighting function*) a poor measure* in and of 
itself* of vortical rooolving powor. Wo nooi only comparo N 
cnannol* of tho NIMBUS III 8IRS Instrument with N similarly 
positional channolt of tho HIR8 or AMTS inotrumont* to 
appreciate thi* fact. Tho channel* of the SIR8 instrument •• 
overlapped that it could not resolve vertical structure as 
well as a similar number of channels of the current HIRS 
instrumentation. Thus* the width of the transmi ttance 
weighting functions as well as the least count* or number of 
separately located channels* af*ect the '.ortical resolving 
power of the system. 

Noise in a satellite borne sounding Instrument will also 
degrade its vertical temperature resolving power just as it 
would in a radiosonde system. Such Instrument noise is fac- 
tored explicitly into the Bachus-Oi lhert-Conrath resolving 
power analysis where there is a trade-off between the char- 
acteristic width of integral averaging kernels and the sen- 
sitivity of a retrieval to instrument noise. 

A very Important aspect of vertical resolving power is 
the algorithm for transforming satellite radiometer measure- 
ments into temperature profile estimates. Unlike an jjx gl tit 
direct measuring system, such as the radiosonde* where the 
sensor signal is dir .ctly related to the temperature* the 
satellite radiation measurement is indirectly and Integrally 
related to the vertical distribution of temperature and must 
be unraveled by some sort of Inverse solution to the integ- 
ral relationship. The solution is not unique and there are 
many possible temperature estimation algorithms* each per- 
forming some sort of c omp lei transformation of radiance into 
temperature. The importance of this step is dramatized by 
the fact that one can construct an algorithm which trans- 
forms all radiometife measurements of a satellite sounding 
instrument into a coarse vertical mean value of temperature. 
In such a case* the resolving length is the entire depth of 
the atmosphere. (For example* the inverse solution given by 
Conrath (1772) for the case q**0 does essentially this.) In 
such an example* what inherent resolving power there may 
have been in the satellite radiometer measurements has been 
deadened by the temperature retrieval algorithm. The mini- 
mum information solution of Smith* Woolf and Fleming (1972) 
can be tuned to merely reproduce the first guess field 
itself* thereby completely nulling the resolving power r.f 
the radiometer measurement. To one extent or another* a?l 
temperature retrieval algorithms will modify what may be 
considered an inherent resolving capability of a satellite 
sounding instrument. 
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It it for tho obovo reasons that wo contitfor tho total- 
lito temperature founding problem to involvo two separate 
components of o sounding "system”. Tho firtt it tho intru- 
ment with itt inhoront choroctorittica of tpoctrol bandpass* 
Itott count* «gnal sensitivity* noito* otc. Tho second it 
tho mathematical algorithm utod to map tho measured radl- 
ancot into ottimatoo of tho tomporaturo profile* with itt 
intrinaic smoothing and mathematical stability characteris- 
tics. 

Thoro it ono final aspect of the issue of vertical 
resolving power which it important. That it* the vertical 
resolution of a priori data used in the sounding. Rodgers 
(1976) hat spoken to this point by showing that 2 priori 
data increases the overall vertical resolution of tempera- 
ture profiles derived from satellite measurements. In an 
even simpler eiample* the minimum information method* tuned 
as discussed above* yields a vertical resolution for the 
sy ste m which it just the vertical resolution of the a priori 
data. In practice* it is upon this initial vertical resolu- 
tion that the sounding system may improve* or degrade. Even 
such gross information at correct tropopause level can make 
substantial improvements to retrieval accuracy (Thompson and 
Wolski (19771* Westwater and Orody (1981)1. In our studies* 
the issue of the resolution of a priori data it specifically 
and carefully r emove^ from our analysis so that we may 
address only the question of the latent characteristics of 
the instrument and a retrieval algorithm* which incorporates 
only a gross a priori depiction of the vertical structure of 
the atmosphere. 

3. Channel characteristics of H1RS and AMTS 

Figure 1 shows the transmi ttanc e weighting functions* 
dH <\ )/di* for the HIR3 and AMTS instruments. These func- 
tions were computed for a standard atmospheric temperature 
profile using “iransmi ttance data kindly provided by J. Suss- 
kindi NASA/CLAS* and are normalized by their maximum values. 
Me show here 11 »elected channels of each instrument in the 
4. 3/qm or 13,qm C02 absorption band. The circled labels indi- 
cate 4. 3/4* channels. These >ub-sets of channels were 
selected so that thare are the s^ime number of 4. and 15 m* 
channels in each set and their distribution with height is 
as uniform as possible. The AMTS channels are generally nar- 
rower than the HIRS and there would appear to be more infor- 
mation gathered above lOOmb by the AMTS measurement chan- 
nels. We note the region from about 1 50mb to 300mb in which 
each Instrument a ppears to have lower vertical sensitivity. 
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Figure 2 ihowi a comparison of the radiative traniNr 
kernels. <dl( . TC i > )/dT*< i >><dr< ?i >/dx>. which hovo boon 
computed for o standard atmosphere. and normalised. Thtso 
functions art particularly rolovant to tho ompirical tost of 
vortical rosolving power discussod in tho companion papor 
(Thompson. (1901)) whoro they aro usod in linoariiod lnvorso 
rotrioval algorithms. Tho tomporaturo dopondonco of tho 
transf ormat i on to radiativo transfor kernels causos somo 
significant difforoncos botwoon tho functions in Figures 1 
and 2. Channels 1-6 of tho AMT8 aro only slightly modified 
but 7-11 are drastically changed. Tho levels at which tho 
radiativo transfor kernels peak aro significantly low** than 
corresponding levels of tho transml ttanc o weighting func- 
tions. These latter channels are also sharpened when ren- 
dered as radiative transfer kernels. However, the combined 
effect is to create an information "gap" around 200mb for 
the AMTS Instrument. For the MRS. channels 1-3 are only 
slightly modified, channels 4-9.11 have their peaks signifi- 
cantly lowered and sharpened. Channel 10 completely changes 
character in the transformation. The main peak is lowered 
from near the lOOmb level to 300mb with only a secondary 
peak near lOOmb where the transmittance weighting function 
peaks. Also. the radiative transfer kernel for Channel 10 
shows sensitivity at very high levels (‘Mmb) which adds to 
the information collected in Channel 1 of the H1RS. 

The distribution of radiative transfer kernels for the 
HIRS also shows an information gap around 200mb although 
Channels 4 and 10 help to supply information in this region. 
One of the most significant differences in the sets of ker- 
nels Is that MIRS has 8 channels strongly sensitive below 
200mb and 3 or 4 strongly sensitive above. AMTS shows 6 
channels sensitive to the higher levels and 9 sensitive 
below nOOmb. 

To further analyte the information measured by the two 
instruments. we performed a singular value decomposition of 
the weighting functions and radiative transfer kernels of 
each. Figure 3 shows the cumulative sum cf normalized, 
ranked eigenvalues as a function of indei for the HIRS and 
AMTS. The solid curves refer to c harac ter i s t i c s of the tran- 
smittance weighting functions while the dashed curves refer 
to the radiative transfer kernels. The results show that 
about 99X of the information gathered by a HIRS instrument 
could be explained by only three or four orthogonal eigen- 
functions of the transml t tanc es while six to seven eigen- 
functions of the AMTS would be required to explain the same 
amount of information. This implies that the AMTS instrument 
should possess more* vertical resolving power since the 
information in the weighting functions is spread over more 
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orthogonal dimensions of radiance space. 

The situation with radiative transfer kernels for the 
HIRS and AMTS is quite different. The information distribu- 
tion for the kernels is very nearly the same for each 
instrument. The AMTS instrument shows perhaps only one more 
significant eigenfunction than the HIRS when information is 
rendered in terms of radiative transfer kernels. For each 
instrument* only about 3 or 4 eigenfunctions are necessary 
to account for most the information contained in the radia- 
tive transfer kernels. 


4. Theoretical analysis of vertical resolution 

The work of Backus and Gilbert (1968) opened a theoreti- 
cal approach to the subject of resolving power in inverse 
geophysical problems. Conrath <1972) adapted the work of 
Backus and Gilbert specifically to the satellite temperature 
retrieval problem. I.. this section* we present an analysis 
of vertical resolving power following this theoretical 
approach. 


In the notation of Conrath* the measurements of spectral 
radiation by a satellite radiometric device may be repre- 
sented by the equation 
Xr 


/)!/ =/ K; 


< y ) AT ( y ) dy + 


6 i 


i-1. 2, 


. N 


< 1 ) 


where I ^ is the measured radiance at each of N frequen- 
cies. |)^ « less the surface contribution* which is contami- 
nated by errors €£ » T(y> is the atmospheric temperature at 
some vertical level y* Kj(y) are the radiative transfer ker- 
nels defined earlier for the sounding frequencies t)^. The (A) 
indicates that differences have been taken betw.en true 
and background values* (e. g. ^T < y ) =T < y ) -T*( y ) i j±i = I £ --I* ). 

If one considers the class of linear temperature retrieval 
algorisms. yielding estimates of temperature at levels x* 
expressible by 


✓\ N 

AT< x ) * 2 a £ (x> AI t 
i-t 

then, the retrieval estimate* 
to the ambient profile. 

^*r 

/'N 
AT ( 


iff < x ) , 
by 


< 2 ) 


is integrally related 


•> ■ / £>. 

•i. '*> 


<x> K,-<y> AT<y> dy 


<3> 


In (3)* the estimate of temperature at each level* x* in 
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the atmosphere it given by a weighted average of the entire 
ambient vertical profile* AT(y). The fidelity of this esti- 
mate depends on the narrowness of the averaging kernel# or 
spanning function# A(x.y) ■ ?"• ,(«) K 4 (y>. and the amount 
of transfer of measurement errors# • into the temperature 
retrieval. The error transfer also depends on the narrowness 
of the averaging kernel so that resolution and error trans- 
fer are not independent but# rather# play against each 
other. 


The Bac k us Gi 1 ber t-Conrath ( BGC ) analysis defines two 
properties of the sounding system. The characteristic verti- 
cal width of the averaging kernel# or "spread", is given by 

/r 


six) - 12 


/ 


< x y > <A( x. y ) > dy 


(4) 


measurement# 


The errors of radiance 
min ate the retrieval through <2). 
errors into temperature retrieval 
an error variance# 

N N 


£, < *# in (1) will conta- 
The mapping of measurement 
space may be expressed as 


0-*<.) - <1/M> 6, £*>;<»> 


*?•/ <«l J*l 


(5) 


where we imagine that M soundings have been made with 
the measurement error in the i'th radiometric channel during 
the k'th sounding. In matrix notation# the spread and Tor 
variance may be written 


s(x ) 


T ^ ^ 

(a) <S> (a) 


A 

c r < x > 


( a ) i E > ( a ) 


where (a) is a column vector whose elements are the selected 
retrieval coefficients a*(x)« and matrices Jaix) and are 
given by 

Xr 

Six) (S dJ > « 12 / ( x -y )** ( y > Kjiy) dy 

■l# 




M 

< i/m) yz 





Therefore# the spread function depends on instrument design 
(K^(y>) and retrieval coefficients a 4 ix) whfle the noise 
transfer depends on the retrieval coefficients and the noise 
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covariance matrix of th* radiometric instrument. 

For the theoretical analysis of resolving pouisri ui com- 
pare the two instruments - \.he HIRS and AMTS - using appro- 
priate matrices ?3( x ) and §T for each. 

To proceed with the analysis# we specify the radiative 
transfer kernels* K;(y)» and noise characteristics* <eT)# for 
each instrument. Figure 2 shows the radiative transfer ker- 
nels K^(y) for the infra-red sounding channels selected for 
the HIRS and AMTS test. For the noise character istic s# we 
specify two noise covariance matrices for the AMTS and three 
for the HIRS. All are diagonal and the first in each set is* 
what we shall call the IX set. The noise values are approxi- 
mately IX of radiance values emanating from a standard atmo- 
srhere. The values are actually rounded so that there is a 
single characteristic value for the 1 5/^m channel* and a sin- 
gle value for the 4. 3x^m channels. The second set of noise 
values for each instrument were provided by L. McMillin of 
NCJAA/NESS and represent a more realistic set of noise values 
based on NQAA/NESS investigations. A third set of noise 
values are used in the HIRS test which simulate the noise 
levels of a HIRS instrument assuming its detectors wire 
cooled in a way similar to the proposed AMTS sensors. McMil- 
lin<*) has estimated that cooling the HIRS detectors in a 
manner similar to that ’specified for the AMTS would reduce 
the noise values by a factor of ap pro x imat e 1 y 2. 5-3. 0. We 
use here a reduction factor of 2. 5. Table 1 summarizes the 
noise characteristics of the simulated instruments tested in 
this paper and the paper to follow. 

Figure 4 shows spread and noise variance as a function 
of height for the HIRS1 instrument with IX radiometer noise. 
Spread is shown here in terms of geometric length rather 
than a pressure related variable. Graphs are shown for vari- 
ous values of the retrieval parameter q. For comparison# 
Figure 5 shows corresponding curves for AMTS1. In all 
cases# varying the retrieval parameter has a larger effect 
on reducing noise variance than on increasing the spread 
function. The spread function curve for q=1.000 is a mini- 
mum limit on spread function values. The spread function 
varies from a few kilometers in the troposphere to greater 
than 15 or 20 kilometers in the stratosphere. At 500mb# the 
curves indicate that in order to reduce noise variance to a 
tolerable level -say# less than IK - the spread function 
must take on a value roughly twice as large as the 


(*) Personal communication 
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theoretical ly minimum value. (About 5km for q«0. 990 rather 
than 2.5km for q*1.000 ) 

The 1% AMTS1 instrument shows a better stability of its 
spread characteristics over the range of solution parameters 
used. For example* the degradation in spread function which 
occurs whtrn q is reduced to 0.990 from 1.000 is not star- 
tling even though the associated reduction in noise variance 
is significant. The spread function for AMTS1 is signifi- 
cantly smaller than for HIRS1 at high levels above lOOmb but 
the two instruments show comparable response in the tropo- 
sphere. The most notable feature of the AMTS1 curves is the 
apparent degradation of resolving power in the layer 
200-400mb which is not exhibited in the HIRS1 channels. 

Figures 6 and 7 show similar spread-var iance curves for 
HIRS2 and AMTS2 instruments which incorporate more realistic 
values of radiometer noise. An important thing to notice is 
that reducing instrument noise certainly improves the noise 
transfer characteristics of an instrument but does not lead 
to order of magnitude improvements in the spread function 
characteristics. This is particularly true for AMTS2 as com- 
pared with AMTS1: the spread function curves. if super-im- 
posed. are clustered very close together. For HIRS2 compared 
with H I RSI . there is an improvement of 2 to 3km in the 
spread function between the surface and lOOmb for comparable 
noise variance levels. (E. g. compare HIRS1 with q-0. 990 with 
HIRS2 with q=0. 998. ) As with HIRS1 and AMTS1, the AMTS2 
theoretical spread c harac ter i s t i c s are superior to HIRS2 in 
the stratosphere. comparable to HIRS2 in the lower tropo- 
sphere but worse than HIRS2 in the layer 20Q-400mb. 

At the suggestion of L. McMillin. NOAA/NESS. we also 
tested a simulate^ cooled version of the HIRS instrument 
wherein radiometer noise levels were reduced by a factor of 
2.5 of current noise estimates. The resulting spread and 
variance curves are shown in Figure 8, Such cooling would 
produce a modest improvement in the theoretical spread char- 
acteristics! about l~2km in the stratosphere and less than 
that in the troposphere. By this analysis. the HIRS3 would 
be more competitive with an AMTS2 in the stratosphere and 
comparable in the troposphere. 

As a convenience to the reader, spread and error variance 
for the five simulated sounding instruments are compared in 
Figure 9 for values of the trade-off parameter. q » leading 
to similar traces of instrument noise transfer. 
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3. Summary and conclusions 

In this paper* we have given a discussion of the concept 
of vertical resolving power of a satellite temperature 
sounding system. Me have noted that the vertical resolving 
power should depend on the number* separability* and sharp- 
ness of transmittance weighting functior.s* the noise levels 
of the instrument channels* the resolution of a priori 
information used in the retrieval* and* importantly* on the 
retrieval algorithm used to transform radiance information 
into temperature information. Me suggest two independent 
methods of assessing vertical resolving power* one of which 
was carried out in this paper* the other appearing in a com- 
panion paper. 

tn the present study* w n have obtained measures of verti- 
cal resolving power using th ? theory of Backus and Gilbert. 
That is* we have presented distributions of the spread func- 
tion* corresponding to the Backus-Gilbert integral averaging 
kernel* and the error variance function representing the 
transfer of radiometer noise into the temperature retrieval. 
Such results have been presented for three simulated ver- 
sions of the HIRS instrument and two simulated versions of 
the AMTSi the versions differing only in the assumed noise 
levels of the radiometric channels of the instruments. 

In analyzing the transmittance weighting functions of the 
HIRS and AMTS instruments* we find that the AMTS has one or 
two more pieces of significant information in its C02 chan- 
nels than does HIRS. This conclusion is based on a singular 
value decomposition of the transm i t tanc e weighting func- 
tions. However* when the instrument characteristics are ren- 
dered in the form of temperature sensitive radiative trans- 
fer kernels* we find that the amount of information 
collected by each instrument is very nearly the same* but 
with AMTS still showing a slight advantage. 

In the spread function analysis* we find that the optimum 
parametric condition (q,-l> leaiing to averaging kernels with 
minimum spread generally allows too much noise transfer for 
both instruments. When the trade-off parameter is set to 
reduce the noise transfer to tolerable levels* the corres- 
ponding spread function increases a significant amount. This 
tuning degradation is much less severe for the higher preci- 
sion versions of the sounding instrument where noise levels 
are extremely small to begin with. 

Generally speaking* the theoretical analysis of vertical 
resolving power indicated that the AMTS instrument should be 
superior to a HIRS instrument at high levels above 200mb 
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where the spread function for AMTS is only 50-60% of that 
for HIRS. The two instruments should be comparable below 
500mb having spread function of order 2 or 3 km. In the 
layer from 200-500mb< the AMTS instrument shows a disturbing 
surge in spread function not evident in the HIRS instrument. 
The theory yields a spread function for the AMTS from 3 to 
10 km in this layer* but only 4 to 6 km for the HIRS. 

Now Newman (1979) and Newman and Sagan (1978) discussed 
some deficiencies of the Backus-Gilbert method of analysis 
of vertical resolving power. Ir particular* they have shown 
anomalous behavior of the spread function to the extent that 
as defined here* the spread function may not be an accurate 
measure of the resolving ability of the averaging kernels. 
In the context of our present comparison of HIRS and AMTS* 
this issue is important especially in light of the rather 
large values of spread function for AMTS in the 200-400mb 
layer. 

To help resolve the issue* the author has developed a 
test for vertical resolving power which deals specifically 
with retrievals of temperature signal structure. This empir- 
ical test indicates that the vertical resolving power of 
AMTS around 300mb is much better than indicated by the BGC 
analysis. This test, with a comparison with Backus-Gilbert 
theory* is discussed in the paper which follows. 
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TABLE 1 


Charac t*ri*t ict of the simulated HIRS and AMTS instrument*. 



HIRS 



AMTS 
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— 
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Channe 1 ICenter 

HIRS! 

HIRS2 H1RS3 1 

! Center 

AMTS1 AMTS2 

1 


(Noise Values)* 1 

t i 

(Noise Values)* 

1 (cm* 1 > 

( I*/.) 

(Current) (Cooled) 1 

1 ' ( cm~* ) 

(17.) (Current) 

! 


! 

1 i 




1 

668. 60 

0. 63 

C. 82 

0. 328 

1 

1 

667. 50 

0. 50 

0. 222 

2 

679. 05 

0. 63 

0. 15 

0. 060 

l 

668. 70 

0. 50 

0. 220 

3 

689. 70 

0. 63 

0. 1 1 

0. 044 

1 

667. 00 

0. 50 

0. 220 

4 

703. 80 

0. 63 

0. 08 

0. 032 

1 

666. 00 

0. 50 

0. 222 

5 

716. 70 

0. 63 

0. 05 

0. 020 

J 

652. 75 

0. 50 

0. 222 

6 

731. 85 

0. 63 

0. 06 

0. 024 

1 

646. 75 

0. 50 

0. 250 

7 

2192. 50 

0. 0044 

C. 0011 

0. 00044 

t 

2383. 74 

0. 0025 

0. 000282 

8 

2211. 65 

0. 0044 

0. 0012 

0. 00048 

i 

2386. 10 

0. 0025 

0. 000360 

9 

2237. 35 

0. 0044 

0. 0009 

0. 00036 

! 

2388. 19 

0. 0025 

0. 000293 

10 

2271. 20 

0. 0044 

0. 0007 

0. 00028 

• 

« 

2390. 19 

0. 0025 

0. 000336 

1 1 

2506. 60 

0. 0044 

0. 0005 

0. 00020 

: 

2392. 34 

0. 0025 

0. 000298 
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Figure 1. Transmittance weighting functions for eleven 
comparable channels of the H1RS and AMTS instruments. Cir- 
cled curves refer to channels in the 4. 3mt> C02 absorption 
band# while the others refer to channels in the 15nm C02 


absorption band. 
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ABSTRACT 


An empirical method for assessing vertical resolving 
power of satellite temperature sounders is developed and 
applied to the High resolution Infra-red Radiation Sounder 
(HIRS) and the Advanced Meteorological Temperature Sounder 
(AMTS). The empirical method involves carrying out the 
resolution assessment in a temperature retrieval mode rather 
than in a theoretical mode. Results are compared with the 
standard theoretical Backus-Gilbert analysis. 

Using two well known* physical inverse retrieval 
algorithms* the empirical test demonstrates that the Backus- 
Gilbert spread function does not represent the minimum 
reparation of resolvable signals in the atmospheric 
temperature profile. We find that the HIRS and AMTS should 
!’,ave similar vertical resolution in the troposphere with 
vertical resolving length around 2 or 3 km, quite consistent 
with theory. The AMTS exhibits a slight advantage over the 
HIRS above 200mb, a result which is qua 1 i ta t i ve 1 y consistent 
with theory also. In the layer 200-500mb, the two 
instruments appear to have very similar resolution 
characteristics with AMTS showing a slight advantage over 
HIRS. This is totally contrary to theoretical results which 
would indicate that AMTS should be very poor in this layer. 
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1. Introduction 

In the previous paper (Thompson (1981))# the author dis- 
cussed the vertical resolving pouter of the HIRS and AMTS 
satellite temperature sounding instruments from a theoreti- 
cal viewpoint. The spectral channel and noise characteris- 
tics of each instrument were subjected to an analysis der- 
ived from the work of Backus and Gilbert (1968*1970) and 
Conrath (1972). In this paper> we wish to expand upon the 
study by deriving independent measures of vertical resolving 
power of satellite sounders which are based on the actual 
retrieval of fine scale structure in the temperature pro- 
file. 

Thompson# Eom and Wagenhofer (1976) conducted an empiri- 
cal test of the N0AA/VTPR sounder in a retrieval mode. In 
that study# a single smooth perturbation on a mean tempera- 
ture profile was retrieved in a simulation mode and the 
smoothing ch arac ter i s t i c s of the retrieval system were ana- 
lyzed. There are at least two major limitatiuns to that 
study. First# it involved a sounder whose spectral channels 
and measurement noise c harac t er i s t i c s are far inferior to 
current instruments# such as HIRS and the proposed AMTS. 
Secondly# the retrieval of a single smooth input perturba- 
tion does not tell us whether two or more closely spaced 
p er turbat i ons could be 'sensed by the satellite system. In 
this paper# we will generalize the simulations to come clo- 
ser to the real issue of vertical resolving power. 
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2. Empirical resolving power 

The Ba c kus-G i 1 b er t-Conra th analysis provides very useful 
information about the transmi ttance weighting functions of a 
radiometric sounder. However# Newman (1979) has discussed 
some limitations of this approach and has demonstrated that 
one may obtain spurious values of spread function which do 
not accurately reflect true vertical resolving power. More- 
over# when one uses such a satellite sounding device in 
practice# one normally attempts to retrieve temperature 
structure from the radiance measurements. The fidelity of 
those retrievals certainly depends on the nature of the 
transmittance functions but also is very sensitive to the 
range of difficulties with the inverse problem. Thus# it is 
relevant to develop a resolving power test which is carried 
out in the retrieval# or inverse mode of the problem. 

In developing an empirical retrieval resolving power 
test, we appeal to concepts of Lord Rayleigh. The resolving 
length of an optical system# such as a telescope# is the 
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minimum separation of two point sources (sag* stars) which 
could just be detected as two distinct signals by the tele- 
scope. Adapting this concept to a satellite sounding system* 
one may define the resolving length of a sounding system as 
the minimum separation of two temperature signals along a 
vertical profile which can just be retrieved as two separate 
signals by the retrieval system. It is important to reem- 
phasize that the retrieval system consists of both the 
instrument* with its transmi ttance functions and noise cnar- 
acteristics* and the mathematical algorithm for analyzing 
the measurements. 

In the version of the empirical test of vertical resolv- 
ing power to be discussed here* we wish particularly to omit 
the issue of the resolution of £ priori temperature informa- 
tion which might be used for a retrieval. Although Rodgers 
(1779) has shown that a priori data can be used to further 
decrease the spread function below the minimum values estab- 
lished by the Backus-Gilbert theory* we wish to focus here 
on the resolving power of the radiance measurements alone 
without adding any additional information except a knowledge 
of the mean temperature profile and a good model of atmos- 
pheric transmi ttance. For this reason* we sought a 
retrieval algorithm which does not use a priori statistical 
data but does explicitly use the atmospheric transmittance 
function information. In fact* in our tests we used two 
algorithms meeting these requirements: A Backus-Gilbert-Con- 
rath (BGC) retrieval algorithm and a modified minimum infor- 
mation algorithm* called Moderate Information (MI) here* 
which is patterned after Foster (1961) or Smith* Woolf and 
Fleming (1972). This algorithm uses different noise levels 
for each channel of the radiometric device and a single 
(temperature variance) tuning parameter* s 4 . Further* we 
biased the problem towards the ideal so that we could obtain 
upper limits to system c apab i 1 i t i es. Accordingly* no cloud 
effects or other error contamination of the simulated radi- 
ance meaurements were considered* although the noise charac- 
teristics of the instruments are explicitly included in the 
retrieval algorithms and the convergence criteria. 

The simulations to define a resolving length were con- 
structed as follows. Two carefully defined gaussian temper- 
ature signals separated by a known vertical distance were 
super- imp osed on a standard atmospheric temperature profile. 
Exact radiances for each channel of the HIRS and AMTS 
instruments were computed for this profile. These radiances* 
in turn. were submitted to each of the retrieval algorithms 
mentioned above to produce a retrieved estimate of the dou- 
ble peaked s ignal. 
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For each retrieval* the first guest field was a tempera- 
ture profile everywhere PM cooler than the base profile upon 
which the temperature signals were super-imposed. Thus* the 
first guess profile contains no information about the temp- 
erature signal. Each of the two retrieval algorithms 
involve the transmi ttanc e weighting functions and a covari- 
ance matrix representing radiometer noise information. The 
noise covariance matrix was assumed diagonal with noise var- 
iances for each channel for each particular instrument as 
specified in Table 1 of Thompson (1981) distributed along 
the diagonal. Preliminary tests were run on each of several 
versions of each algorithm using different methods of han- 
dling temperature updating and linearization. Generally 
speaking. the retrieval algorithms Wf>n iterated until the 
root mean square deviation between actual and retrieved 
temperature profiles reached a first minimum. This conver- 
gence criterion is not identic 1 to requiring that radiances 
corresponding to the retrieval converge to the measurements 
to within the instrument noise level for each channel. nor 
to requiring that the radiances converge in a least squares 
sense. When linearized. physical matrix inverse algorithms 
are used with high precision instruments. such as HIRS or 
AMTS, it maj well happen that the algorithm will diverge, in 
the rms temperature sense* before reaching an acceptable 
l%vel of convergence in a radiance sense. This, presumably* 
is due to a cong 1 omeratton of difficulties related to the 
1 i near i za t i on approximation applied to such high precision 
instruments. and to the implicit i 1 1 -c ond i t i on i ng of the 
linearized solution. Although the temperature convergence 
criterion could not be used in practice, we feel comfortable 
in using it in this analysis of optimum resolving power. 
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3. The simulation approach 


Figure 1 illustrates the general nature of a determina- 
tion of empirical resolving length for a given instrument 
(AMTS). a given retrieval algorithm (DOC). at a particular 
atmospheric level (lOOmb). In Figure 1. the base atmosphere 
has been subtracted from both the input and retrieved sig- 
nals so that only the perturbations are shown. The input 
perturbation is a sup er -p os i t i on of two gaussian curves each 
of amplitude 10K and half-width of 8 levels <8 equal incre- 
ments of p**2/7) separated by a distance L. The separation L 
is set at some small value. radiances are computed and the 
retrieval algorithm is iterated until convergence. At that 
point. a determination is made of whether the input signal 
has been resolved. This resolution criteria is as follows: 
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<«> There most be tuio local maxima with an 
intervening minimum in the retrieval "near" the 
input perturbation. 

<bj Given (a)» if themaiima are denoted <%, ) 
and ( t^ ) and the minimum (t # )« then the perturba- 
tion ie consider ed to be resolved if 


(l/2H(t, > + <t A >> - (t. > > (e) 

where <e> is an rms deviation between input and 
retrieved signals "away from" the perturbation. 
The region over which (e) is computed is from the 
surface to lOmb but excluding the region between 
( x T ) and <x B ) where 


< x r > =* ( x^ ) - H 
<x a > =* ( x, > + H 

(c) Given that (a) and (b) are met# a simple 
measure of the quality of vertical resolution of 
the input signal is given by £ =P/Q where P is the 
rms deviation between the input and retrieved sig- 
nals between lOOOmb and lOmb# and Q is the rms 
signal strength of the input signal calculated 
between lOOOmb and lOmb. 

In Figure 1# results for three successively larger values of 
L are shown. It is clear that the input signal is not 
resolved# by the criteria above# for the first two experi- 
ments. In the third# the input signal has been separated by 
a sufficient distance so that the retrieval system resolves 
it although with reduced amplitude# shifted center# and 
separated maxima. 


4. Retrieval method test 


In the course of completing the research reported here# 
we encountered many difficulties in carrying out retrievals 
using the two linearized physical inverse methods applied to 
the high resolution# low noise sounders. The convergence 
criterion mentioned earlier had to be used. instead of one 
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requiring that r.adiances converge to simulated measurements! 
because of certain over-relaxat ion properties of the inverse 
solutions. In some experiments# the last channel to converge 
would converge so reluctantly - or possibly not at all - 
that the algorithm either would diverge altogether! in the 
sense of rms retrieval error# or would swing away only to 
begin a slow convergence again. Twomey (1977# Section 9.1) 
showed a similar example of residual error oscillation! alt- 
hough with a direct rather than physical inverse solution. 
Further# in some cases# a solution satisfying the radiance 
convergence criteria would be worse# from the vertical reso- 
lution criteria# than one which had not converged in the 
radiance sense. The difficulty lies# in part# with the 
1 inear i xation used in the retrieval algorithms applied to 
instruments with low noise levels# and possibly due# in 
part# to the algorithm for temperature updating of the tran- 
smittances. Put simply# a sequence of linearized solution 
iterations incorporating temperature correction of the tran- 
smittance* may not converge# within the very small toler- 
ances of the HIRS or AMTS# to a correct solution of the full 
non-linear problem. Further# a linearized solution which 
does produce radiances within noise levels of the measure- 
ments may not# in fact# be optimal in temperature space. 

For the empirical test of the AMTS and HIRS instruments# 
we began by using iterated retrieval algorithms which start 
with the first guess temperature field and with transmit- 
tances also computed for this profile. At each iteration 
step# transmi ttanc es were recomputed using the NASA/Gt-AS 
algorithms# provided to us by J. Susskind# and the new esti- 
mate of temperature. In this paper# this is called Version 
1 of either alqorithm. In our experiments# we experienced 
inordinate difficulties with this version of the MI 
retrieval algorithm applied to HIRS2 and HIRS3. Of course# 
to obtain the sharpest possible vertical resolving power# 
the algorithm was tuned as close to a high resolution limit 
as feasible. Generally# the algorithm was mostly unstable 
for the higher values of s* tried. When s* was lowered# so 
as to stabilize the solution# the resolving power was suffi- 
ciently reduced so as to produce no meaningful determina- 
tions. The moderate information r e g u 1 ar i z at i on of the 
inverse problem does not seem to be an optimal solution for 
the HIRS when very low noise levels are incorporated into 
the soluti oni at least in the case of these simulations. 
Therefore# we used an alternate version of the MI solution 
for the HIRS2 and HIRS3 in which kernels corrected to the 
true input perturbation temperature profile are used in the 
iterative solution. (This is called Version 2 of either 
algorithm. ) For the AMTS2# we used both versions of the MI 
algorithm for comparison. Version 2 ought to represent an 
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ideal optimum can tine* it contain* correct kernel* from 
the out*et even though the algorithm operate* on the «mooth 
first gu««s temperature field at the first iteration. While 
Version 1 does not achieve these optimum results* it gives 
some indication of how close to optimum one can get operat- 
ing objectively with no 1 priori information whatever. 

In our trials* each solution was iterated until the rms 
deviation between input and retrieved temperature reached 
its first relative minimum. (Generally speaking* this 
"final* 4 solution does not necessarily produce radiances 
within the noise level of the input radiances* but does pro- 
duce the best solution in temperature space. ) Also* optimum 
retrieval parameters were determined for each instrument and 
each algorithm by extensive tuning tests. 

9. Empirical resolving length results 

Figure 2 shows empirical resolving length as a function 
of height for the HIRS2 instrument. The two curves shown 
were obtained by using each of the retrieval algorithms dis- 
cussed above. The dots appearing on the diagram represent 
smallest resolvable differences at each level in our deter- 
minations of empirical resolving length. The successive dots 
at a given level rrpresent successive increases by two lev- 
els (two increment* of ’p**2/7) of the separation of the 
input perturbations. Also* shown as a dotted line is the 
theoretical spread function of Backus and Gilbert* taken 
from Thompson (1*781)* for the same q. parameter as used in 
the BGC retrieval. Though the BGC retrieval results were 
limited to the troposphere* they follow very closely the 
theoretical spread function. The empirical resolving length 
using the BGC retrieval is smaller than the theoretical 
spread function in the layer 200-50Gmb and only slightly 
greater below that. The result for the Mi-Version 2 solution 
shows an empirical resolving length less than the theoreti- 
cal spread function at almost all levels. At 200mb* for 
example* the MI-2 empirical resolving length is only about 
half as great as the theoretical spread. These results give 
us our first indication that the theoretical spread function 
does not represent the minimum resolving length possible for 
a remote sounding device. Even though the MI-2 is an ideal- 
ized algorithm* it is one which should be approachable in 
prac t ice. 

The panel to the right in Figure 2 shows the normalized 
retrieval error ^ as a function of height corresponding to 
the retrievals in the left panel. The MI-2 retrievals were 
very good at all levels producing an rms retrieval error 
from 1000-10mb varying from 10-35% of the rms temperature 
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input signal strength. Thus* not only were th* input pertur- 
bations resolved* but they were resolved with reasonably 
high quality. The quality of the BQC retrievals was somewhat 
lower and the result near 400mb was of particularly low 
quality* even though it also yielded an empirical resolving 
length less than the theoretical spread function. 

Figure 3 shows a similar display of results for the AMTS2 
instrument. In this series of experiments* we were able to 
obtain excellent* and quite comparable results from Version 
1 of both retrieval algorithms. We also show results for 
Version 2 of the MI algorithm so that the most optimum case 
is shown as well. The curves show that retrievals using the 
AMTS2 instrument come very close to the ideal in the lower 
troposphere with resolving length - by any measure - being 2 
to 3 km. Above 200mb# the empirical results fall short of 
the theoretical result. Even the MI-2 retrievals have verti- 
cal resolving length s i gn i f i cant 1 y larger than the theoretic 
cal spread function. The QGC-1 and MI-1 algorithms yield 
quite comparable results throughout the entire atmosphere 
with the largest discrepancy around lOOmb where the BQC 
algorithm shows a decided advantage. Elsewhere* the MI 
algorithm produced sharper resolution with lower retrieval 
error. 

The most striking feature of Figure 3 is that th* large 
peak in the theoretical spread function between 200-400mb 
appears to represent a gross under es t i ma te of the actual 
vertical resolution of the AMTS2. On the other hand. the 
right panel of Figure 3 shows a substantial increase in the 
retrieval error for both the BGC-1 and MI-1 algorithms in 
this same region. This means that while the input signals 
are resolved - by the criteria set out for this study - the 
quality of the retrieval is not as high in this layer as in 
th* troposphere. Note. however* that th* idealized Ml-2 
algorithm yielded higher quality retrievals through this 
same region. This is near th* region where the distribution 
of th* AMTS radiative transfer kernels leaves an apparent 
information gap and th# ability to retrieve fine structure 
in th* profile was expected to be degraded. (*> The large 
discrepancy in theoretical and empirical measures of 


(*) It should be emphasized here that our test of AMTS is 
not altogether fair since we used only 11 sounding channels, 
fewer than are actually available. Susskind* in personal 
communication. has pointed out that an additional AMTS COj 
channel near 300mb would be available to an operational ver- 
sion of the instrument. 
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vertical resolution in this region call* for a deeper look 
into the matter. 

Figure 4 shout the AMTS2 averaging kernels appropriate to 
the 300mb level computed for the DOC and MI retrieval coef- 
ficients . The amplitudes of these curves are arbitrarily 
normalized so that only the shape is important here. The 
kernel for the DOC coefficients is separated into two promi- 
nences* one slightly above 200mb and one between 400-500mb. 
The 300mb level itself is a local minimum for the DOC aver- 
aging kernel. Figure 5 shows the seri»« of DOC retrieval for 
perturbations centered at 300mb# and for separations ranging 
from 6 levels (1.94km) to 16 levels (9. 16km). The resolution 
criteria is met for a separation of 16 levels* nearly 
exactly the separation of the prominences in the DOC averag- 
ing kernel. This suggests that this experiment may only be a 
"resonance" phenomenon where perturbations are only ret- 
rieved when they become synchronized with prominences or 
side lobes of the averaging kernels. However* the averaging 
kernel in Figure 4 for the MI retrieval coefficients is 
sharp* at 300mb* with no such prominences. The corresponding 
series of retrieval experiments are shown in Figure 6 where 
the input signal is also resolved at separations of 14 and 
16 levels: even better than for the DOC retrievals. We con- 
clude* therefore* that the resolution of these signals by 
both methods is an event more fundamental than resonance 
with kernel side lobes. That is* the empirical resolving 
lengths determined between 200 and 400mb are just as signi- 
ficant as those determined at other atmospheric levels. 

The empirical resolving length results for the HIRS3 
instrument are shown in Figure 7. Once again* we had diffi- 
culty in producing stable* highly resolving retrievals with 
the MI-1 algorithm. Even with the BGC algorithm* the inter- 
play between the high sensitivity tuning of the algorithm 
and the temperature correction of the kernels produced 
unstable retrievals at many of the levels investigated. 
Those BGC-1 results obtained are consistent with the Mi-2 
results and we take the pattern of both to be representa- 
tive. By comparing Figure 7 with Figure 2* one sees that 
lowering the noise level of the HIRS instrument should pro- 
duce an increase in the vertical resolving p ower* e v e q 
though the theoret ical s pread f unc t i ons s hown for HI RS2 anjJ. 
HIRS3 in. t hese fi gures are virtual It! indist ingu ishable . The 
empirical resolving length for HIRS3 is significantly better 
than HIRS2 above 200mb and slightly better in the tropo- 
sphere. As with HIRS2* HIRS3 shows consistently better MI-2 
results than are predicted by Backus-Gilbert theory. Even 
the few BGC-1 results obtained suggest that the high empiri- 
cal vertical resolution may not be solely due to the 3 
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priori information in the kernels for method MI-2. 

For convenience to tho reader* tho MI-2 results for oach 
instrument totted horo aro grouped together in Figure 8 for 
instrument comparison. The three instrument* behave compa- 
rably in the troposphere with the HIRS3 generally showing 
the highest quality resolution of input signals and AMTS2 
generally outperforming HIRS2. Between 200mb and 30mb* the 
HIRS3 appears to have a decided advantage over AMTS2 even 
though the theoretical spread function results would lead 
one to believe just the opposite. The present version of 
HIRS ( HIRS2) would appear to be inferior to the AMTS2 in 
this layer. Above 30mb* the AMTS produces higher resolution 
than either version of H1R8. 

6. Conclusions 

In this paper* we have developed and presented a method 
for assessing the vertical resolving pjwer of satellite 
borne temperature sounders which is carried out in a 
retrieval mode. The method is conceptually similar to meth- 
ods of defining resolving power of telescopes and micro- 
scopes. The advantage of this method over the theoretical 
method of Backus and Gilbert is that the retrieval algorithm 
itself plays a strong role in the determination of resolu- 
tion properties. The me'thod could be used for any retrieval 
algorithm* including direct relaxation methods and non-li- 
near methods* since assessment is carried out in temperature 
retrieval space. Moreover* in light of certain known limita- 
tions of the theoretical Bac k us-Gi 1 b er t approach* the empir- 
ical method shows directly - instead of i nf er ent ia 1 ly - the 
minimum separation of signal perturbations on the tempera- 
ture profile which car be resolved by a temperature sounding 
system. 

T l .> most significant disadvantage of the method is that 
the simulation used to define resolving length i i decidedly 
unatmospher ic. On the contrary* isolated signals are ret- 
rieved centered around specific levels instead of a more 
natural super-position of irregular signals such as might be 
encountered in nature. 

Nevertheless* the approach developed here is reasonably 
objective and does yield a baseline measure of the vertical 
resolving power of satellite spectral radiation measurements 
absent the expectation of structure derived from independent 
a priori information. 

The empirical method has been applied to two high tech- 
nology sounding instruments* the HIRS and A' iS* using each 
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of two standard* well known physical invars* retrieval solu- 
tions. These results* in turn* have been compared with the 
standard theoretical assessment given by the Backus-Gilbert 
method. 

W* have found that the empirical resolving length of a 
sounding system may be either smaller than or larger than 
the Backus-Gilbert spread function* depending on instrument* 
retrieval method* and atmospheric level. That the empirical 
resolving length may be smaller than the theoretical spread 
gives evidence that the latter may not comprise the best 
test of vertical resolving power. That the empirical measure 
may often be larger than the theoretical measure gives evi- 
dence that a given retrieval algorithm will not always 
extract the most temperature structure information incorpo- 
rated in the radiation measurements from satellite sounders. 

Based on our empirical tests of vertical resolving power* 
we have found that the HIRS and AMTS instruments should have 
nearly comparable vertical resolution in the troposphere — 
AMTS shou-lng a slight advantage — with vertical resolving 
length varying around 2 or 3 kilometers. The AMTS instrument 
also exhibits a slight advantage over the present version of 
HIRS* from the resolution standpoint* above 200mb with maxi- 
mum advantage occurring near 20mb. If* however* the sensors 
of a HIRS instrument were cooled in a fashion similar to the 
proposed AMTS* our results show that the vertical resolving 
power above 200mb would be substantial ly improved* perhaps 
exceeding the capabilities of the AMTS. The effect of sensor 
cooling for improving tropospheric retrievals by HIRS does 
not appear to be dramatic. In the layer 200-500mb* we find 
a significant difference between our empirical results for 
the AMTS and the theoretical* Backus-Gilbert results. The 
empirical resolving length is much less than the Backus-Gil- 
bert spread function in this layer and is also generally 
smaller than the resolving length for the HIRS. We conclude 
that the spread function for the AMTS is misleading in this 
layer and that the AMTS instrument is superior to HIRS bet- 
ween 200 and 500mb» rather than inferior as would be infer- 
red from Backus-Gilbert theory. 

Finally* it is important to note that we experienced many 
difficulties in conducting empirical tests of vertical 
resolving power using the high precision HIRS and AMTS 
instruments. Some of these difficulties can be traced to the 
linearization of the f undamenta 1 1 y non-li^ecr problem which 
is incorporated into the physical inverse retrieval solu- 
tions used in our tests. The present author* therefore* con- 
siders the comparisons between HIRS and AMTS shown here to 
be tentative and very much in the context of the approximate 
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r*tr i eval methods used. We suggest that a testing approach 
similar to ours be adapted for more sophisticated retrieval 
algorithms before final assessment of the merits of each 
instrument are made. 
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Figure 5. BGC retrieval experiments at 300mb for the 
AMTS2 corresponding to Figure 3. 
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Figure 6. MI retrieval experiments at 300mb for the AMTS2 
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Figure 7. Empirical resolving length and quality for 


HIRS3 
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Figure 8. Comparison of empirical resolving length an 
quality for HIRS2* HIRS3# and AMTS2. 
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Information concerning the ambient atmospheric thermal 
structure can be severely distorted by satellite temperature 
retrieval systems and subsequent analysis algorithms. 
Generally speaking* we know that a finite set of spectral 
radiance measurements bears only a coarse. integrated 
relationship to the ambient temperature structure which 
emits the radiation. Temperature retrieval algorithms may 
smooth* add* or further distort information depending on how 
stable is the algorithm, and how heavily influenced it is by 
i priori data. Further analysis of retrievals in 
preparation for a final meteorological study will generally 
distort the information further. 

In this paper* the transfer and distortion of the basic 
ambient temperature information is discussed. The approach 
is to examine the variance spectrum of information as it 
takes different forms in the retrieval-analysis cycle. 

It is shown that the (total) adjusted variance of 
radiative brightness temperatures measured by the HIRS 
satellite instrument is only about half the variance of 
corresponding atmospheric temperature profiles in middle 
latitudes* and only about 30% in the tropics. For 
statistical (regression) temperature retrievals* about 
43-49% of the i nf ormat i on content is attributable to & 
priori data while for physical retrieval methods* only about 
25-27% is attributable to prior information (such as 
transmittancei), Other retrieval methods lie between these 
extremes. The distortion of information is explicitly shown 
by examining empirical orthogonal functions of corresponding 
RAOB and retrieval profiles. 


1. Introduction 
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Tht impact of remote sensing data on global and mesoscale 
models of atmospheric behavior depend* fundamentally on the 
quantity and quality of new observational information 
provided by tho remote sensing device. In particular* 
satellite based soundings of temperature and moisture should 
have an Impact on diagnostic and prognostic models 
proportional to the amount of new observat i onal information 
imbedded in the retrievals. However* owing to the well known 
difficulties with the inverse problem* satellite derived 
soundings are most always a complex blend of the new 
rod i or r tr i c measurements with a priori information on the 
atmospheric structure. This a priori information may involve 
c 1 ima t o I og i cal mean structure* statistical or probab i 1 i stic 
information* in situ measur ements nearby in time and space* 
numerically forecasted structure* and so on. 

In this paper* an analysis of information transfer for 
the satellite sounding problem is presented. Of particular 
interest is the distortion of the atmospheric information - 
the information content of the ambient atmospheric thermal 
field - at it naturally produces the outgoing radiance 
field* and the further distortion of this information as the 
radiation field is transformed back into thermal structure 
information through the retrieval step. The amount and 
quality of new information about the ambient atmosphere 
provided by this satellite sounding process is not always 
clear but is of fundamental relevance to the meteorological 
impact of that observing system. 

Figure 1 schematical ly illustrates how information 
embodied in* say* the ambient temperature structure is 
transformed into different forms of that information in a 
satellite temperature retrieval and analysis system. The 
variables enclosed in ovals in Figure 1 represent various 
versions of the original temperature information which art 
expected to be moderately to highly distorted versions 
compared with the original information in the ambient 
temperature structure. In the r e tr i eva 1 -ana 1 y s i s system 
illustrated in Figure 1* the original information has been 
transformed into various kinds of radiance ensembles or 
temperature estimate ensembles by the various steps in the 
re t r i e va 1 -anal y s i s procedure. In order to know how accurate 
is the depiction of a meteorological phenomena described by 
some variate based on remote radiation measurements* one 
must know how the original information in the ambient 
tempera! ure field has been degraded by the natural radiative 
transfer integration, by the error contamination of radiance 
measurements* an* 1 by subsequent processing and analysis of 
that radiance dat (See. for example* Thompson* Eom and 
Wagenhofer (1976)* Thompson (1982a. 1982b). ) The radiative 
transfer physics and the finite* spectral resolution of a 




radiometer assure that substantial vertical averaging has 
been done to the actual temperature data in transferring 
that informaticn to the radiance domain. Temperature 
retrieval algorithms! often involving a priori or forecasted 
temperature information which is independent of the radiance 
observation* supply their ouin information transfer 
properties so that the retrieval is a distorted version of 
the ambient signal. The distortion is partly due to bias by 
the a priori information used and partly due to the inherent 
mathematical properties of the retrieval algorithm. Further 
analysis of the retrieved temperature profiles contribute 
even more distortion of the information. 

The information can be examined at each step of the 
process illustrated in Figure 1 - either in the temperature 
domain or the radiance domain - and the quality and quantity 
of information existing at each step compared with the 
original ambient temperature information. A systematic 
approach to this comparison is given in section 3. 

2. Previous relevant work 


Satellite meteorologists have sought to retrieve the 
structure information of the ambient t emp er a tur e field by a 
variety of techniques. Retrieval algorithms based on the 
methods of Backus and Gilbert (1963) seek to optimize the 
s'ertical resolution of structure in the temperature field by 
sharpening the integral everaging kernels so as to extract 
the structure information directly from the radiance data in 
some optimal sense. Another method seeks to recover 
structure by optimizing the construction of a first guess 
Meld using# say# a numerical forecast model# and then 
constraining the retrieval to deviate minimal amounts from 
that field# (Smith# Woolf and Fleming (1972) for example). 
Other methods constrain the retrievals statistically so t n a t 
they have the structure of an independent# but presumably 
representative# set of temperature data derived from RAOBs 
(Smith# Woolf and Jacob (19701# Smith and Woolf (1976) for 
example)# or by performing a final i n t er p o 1 a t i on or fitting 
of radiance based temperature information by empirical 
orthogonal functions derived from an independent set# (J. 
Susskind# personal communication). In methods utilizing a 
priori information independent of the radiance measurements# 
it is sometimes difficult to know how much of the final 
information has actually been synthesized from a priori data 
and how much is new information contributed by the radiance 
measurements themselves. 

The issue of information content has been studied by 
several i nve s t i g a t or s . Crosby and Weinreb (1974) analyzed 
information content of satellite retrievals based on 
s ta t i s t i c a 1 1 y optimized solutions and have further shown how 
incorrect independent statistics incorporated in a retrieval 


53 

solution can distort the information produced from radiances 
by that solution. Spankuch# Timofayav and Guldner (1977) 
also showed tha distortion of measurad radianca information 
that can occur if ona usas incorract statistics to produce 
ratr iavals. 

In a discussion of early TIROS-N temparatura retrievals. 
Phillips# McMillin# Gruber and Wark (1979) shouted# among 
other things# that tha variance of satellite derived 
temperatures about their mean is often significantly less 
than the variance of co-located RAOB temperatures about 
their mean for most layers in the atmosphere. The ratio of 
variance of TIROS-N soundings to the variance of RAOBs are 
shown to be as low as 0. 2-0. 3 in the higher layers and 
generally range somewhere between 0.5 and 0.9. This implies 
that actual temperature information gets filtered by the 
satellite sounding system, often severely. 

Schlatter (1981) showed the results of comparisons of 
TIROS-N temperature soundings with an NMC (optimum 
interpolation) analysis of RAOB data. The TIROS-N soundings 
were done using procedures described by Smith. Woolf# 
Hayden. Wark# and McMillin (1979). The NMC analysis# done 
without satellite soundings during the period of Schlatter's 
comparisons# is described by Bergman (1979) and McPherson# 
Bergman. Kistler. Rasch and Gordon (1979). Schlatter's 
comparisons showed a significantly different pattern of 
error between TIROS-N soundings and analyzed RAOB data than 
one normally sees in comparisons of satellite soundings 
directly with co-located RAOBt. The rms differences between 
TIROS-N soundings and NMC analyzed temperatures shown by 
Schlatter did not exhibit the larger values at tropopause 
level that many other studies have shown. (Eg. Fig. 6 from 
Schlatter as compared with Fig. 4 from Smith. Woolf. Hayden# 
Wark and McMillin (1979)). Schlatter showed rms errors of 
about 2K from 700-70 mb with no apparent strong maximum near 
the tropopause. Smith, et al (1979) showed rms deviations 
between TIROS-N sounuings and co-located RAOBs which 
achieved local maxima around 200mb of 2.5- 3. OK. While 
somewhat mixed, the results of Phillips# et al (1979) also 
show local maxima in errors at tropopause level when 
considering deviations between TIROS-N soundings and co- 
located radiosondes. Halem. Ghil# Atlas. Susskind and Quirk 
(1978) examined NIMBUS-6 soundings during two periods of the 
Data Systems Test and also found prominent errors at 
tropopause level ranging to as much as 3K. 

Schlatter's result. when compared with the other studies 
cited above# illustrates the smoothing of the RAOB 
information which is effected by the NMC analysis and which 
leaves a result involving the comparison of two fields each 
of which are smoother than the original RAOB data. The two 
smoothed fields (TIROS-N and analyzed RAOB) seem to agree 
with each other better# most likely# than either would 
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compare with the original RAOB information. That it* the 
original information in the ambient temperature field hat 
been transformed into the smoother information contained in 
the analysis. Further* the radiative transfer physics and 
the retrieval algorithm also transforms the ambient 
temperature structure information into smoothed temperature 
retrieval information. The comparison of these two fields* 
then* is that of two smoothed and distorted versions of the 
original temperature information possessing much more 
struc ture. 

In an earlier study* Schlatter and Branstator (1979) 
showed* in fact* that the degree of accuracy in comparing 
NIMBUS-6 temperature retrievals to analyzed radiosonde 
profiles depends significantly on the amount of smoothing 
performed on the radiosonde data during the analysis cycle. 
If satellite retrievals are subjected to a further analysis 
routine - alone or mixed with radiosonde data - the 
information content of the radiance measurements would be 
further distorted. 

The issue of the transformation of atmospheric 
temperature information through the various stages of 
distortion which it encounters in a particular application* 
is quite important to meteorologists. Satellite temperature 
soundings will continue to be used in observing the 
structure of atmospheric systems. It is important to know 
rather accurately what the relationship is between the 
retrieved temperature structure and the ambient temperature 
structure. Sounding systems such as VAS will be providing 
data from which meteorologist will hope to infer the thermal 
structure of mesoscale phenomena. The inference of structure 
snd variation on that scale will undoubtedly pose an even 
more difficult problem than the laige scale weather 
prediction problem. Moreover* one may generally expect that 
satellite soundings constructed to analyze mesoscale 
systems* or to initialze models thereof* will need to 
possess more detailed* accurate vertical and horizontal 
structure than is new obtained for the global problem in 
order to ha\e a positive impact in that arena. 

In the section to follow* a systematic approach to 
investigating the transfer and distortion of information in 
the t emp era t ur e retrieval problem is discussed. 

3. The flow of information 


To quantitatively analyze the flow of information* it is 
necessary first to define information for this study. A 
large ensemble of RAOB profiles and corresponding satellite 
T-ad i ometer observations will form the basis of the data to 
be analyzed. The RAOB data have been assembled under the 
guidance of N. Phillips* NQAA/i-IMC to serve as the basis of 
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several tests of satellite sounding systems. Using this 
data* we have synthesized radiance observations that would 
be gathered by the HIRS and A(1TS(«) sounding instruments. 
Information for this study is defined to be the variance 
spectrum of empirical orthogonal functional representation 
of the ensembles of either temperature profiles or radiance 
values in each oval in Figure 1. While the data in each oval 
in Figure 1 may be taken as a representation of the data in 
any other oval* it is certainly expected that the 
corresponding variance spectra will be quite different. For 
example* Figure 2 shows plots of cumulative percent 
explained variance plotted as functions of the number of 
empirical orthogonal functions used in representing the 
ensemble of RAOBs and corresponding (synthesized) satellite 
radiances for 200 wintertime observations between 30 and 60 
north latitude. The broad spectrum of RAOB information is 
compressed into fewer significant orthogonal functions of 
radiance. 

Clearly* Fig. 2 demonstrates the significant smoothing of 
information uhich occurs even before a satellite measurement 
is made. A temperature retrieval method based only on the 
radiative transfer physics (that is* not incorporating a 
priori information on the atmosphere's thermal structure) 
should exhibit the same drastic smoothing of information. 

Now * in fact* satellite temperature sounding algorithms 
generally add structure information to the radiance 
information in producing retrievals. For example* highly 
statistical temperature retrieval methods* such as the 
regression method of Smith* Woolf and Jacob (19/0) or its 
structured variant (Smith and Woolf (1976)) will add a good 
deal of structure information by using historical knowledge 
of the atmosphere's thermal structure and its relationship 
to satellite radiance measurements. In fact* since the 
physics of radiative transfer is not explicitly included in 
these algorithms* such methods produce retrievals with 
atmosphere-like thermal structure even though they operate 
on the much smoother radiance fields. In the middle of the 
spectrum of retrieval methods* there are physical- 
statistical methods* such as minimum rms (Foster (1961)* 
Strand and Westwater (1968))* known also as maximum 
probability (Rodgers (1970))* statistical regularization 
(Tichonov* as quoted by Kondratiev and Timofeyev (1970))* 
which blend historical t h er ma 1 information with the 
radiative transfer physics to produce a retrieval possessing 
more structural detail than is provided solely by the 
radiances. Finally. there are methods which are mostly 
physical in the sense that no a priori statistics are used 
at all. Such physical methods include the minimum 
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information invers* method (Foster (1961)* Smith* Woolf and 
Fleming (1972))* and direct* relaxation methods such as 
those by Chahine (1968) and Smith (1970) or Smith and Woolf 
(1981). These physical methods* as with some of the others 
above* do make use of a priori transmittance data* and all 
methods use some a priori first guess field. 

Statistical methods (pure or mixed) add historical 
thermal information to the radiance information in producing 
an estimate of the atmosphere's thermal field. If the impact 
of the satellite observation on some analysis or forecast is 
to be fairly assessed* then it is important to distinguish 
how much of the retrieval is really due to the radiance 
observation and how much is due to the climatology of the 
3 tmo sphere. 

4. The method of analysis of information flow 

Let us imagine an ensemble of RAOBs which we consider to 
defi ne 100*/i of the thermal structure information for a given 
situation. For the present study* we will restrict our view 
to only the vertical dimension so that this original 
information may be described by the empirical orthonormal 
functions of the ensemble (eigenfunctions of the covariance 
matrix of the RAOBs) and the corresponding eigenvalues which 
give the variance described by each orthonormal structure 
function. We consider that the eigenvalues are normalized by 
the total variance of the ensemble (the trace of the 
covariance matrix) so that they represent fractions of total 
variance explained by each eigenfunction. If these 
eigenvalues are ranked from highest to lowest value and 
accumulated* one obtains curves such as those in Fig. 2 
labeled TEMPERATURE. 

Now* suppose radiances for some given satellite sounding 
instrument are synthesized using the radiative transfer 
equation. Then* that ensemble may be similarly decomposed to 
produce a variance spectrum such as that labeled RADIANCE in 
Fig. 2. The variance spectra for the ensemble of retrievals 
may be compared with the spectra for the RAOBs to assess the 
distortion of information as it passes through the radiative 
transfer equation and through the retrieval algorithm. 

In the remainder of this paper* four temperature 
retrieval algorithms distributed between the extremes of 
pure statistical to pure physical will be tested to discover 
their information transfer c h ar ac t er i s t i c s . The methods are 
the Regression Method (Smith* Woolf and Jacob (1970))* 
Statistical R e g u 1 ar i z a t i on (Foster (1961))* Minimum 
Information (Foster (1961)* Smith. Woolf and Fleming (19/2)) 
and the iterative* physical method of Smith (19/0). as 
modified by Smith and Woolf (1981). The Smith iterative 
algorithm is generally a numerica lly stable* physical 
method. It involves no sensitive matrix inverse step* and no 
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a priori statistics except for a first guess profile and a 
priori estimates of atmospheric transmi ttances. The 
regression algorithm is also ma th ema t i ca 1 1 g stable but 
possesses the greatest statistical influence of a priori 
data with virtually no explicit radiative transfer physics 
involved. The statistical r egular i i a t i on method involves a 
matrix inversion of the linearized radiative transfer 
equation - hence* has the possibility of numerical 
instability through i 1 1 -c ond i t i oning — but is efficiently 
regularized through the use of a priori statistical 
information derived from historical temperature data. The 
minimum information special case of statistical 
regular i zation also involves matrix inversion of the 
radiative transfer physics with r eg u lar i i a t i on obtained 
through a minimal amount of a priori statistical 
information} namely# a single# mean value of atmospheric 
temoerature variance. 

In the analysis to be presented here# several important 
simplifications or ap p r o x imat i ons are involved. First# 
satellite measurements were synthesized from RAOD profiles 
without any moisture# ozone or cloud effects. For some 
experiments# random errors were added to these simulated 
measurements while other experiments were conducted with 
"perfect" measurements. Thus# our simulation of satellite 
measurements is much cleaner than actual measurements would 
be. Also# only results for the HIRS sounding instrument are 
included in this paper. Retrievals involving transmittance 
weighting functions made use of exactly the same functions 
as those used to synthesize the radiance measurement. Thus# 
even in iteration cycles of a retrieval algorithm# "correct" 
transmittances were used throughout. Next# the surface 
contribution to the radiance measurements is subtracted 
exactly from the problem so that the integral term in the 
radiative transfer equation is isolated perfectly. The 
minimum information and statistical regularization methods 
apply to linearized versions of the radiative transfer 
equation in which radiative transfer kernels and temperature 
deviations about the current temperature estimate replace 
the transmittance weighting functions and Planck function 
profiles. <See Conrath (1972) for the linearized version of 
the direct problem used here.) These two methods are also 
iterated with the required statistical information 
recomputed at each step of the iteration (either the whole 
temperature covariance matrix about the current estimate# 
for statistical regularization# or the average diagonal 
element for minimum information). Finally. eleven C02 
sounding channels of the HIRS instrument are used in this 
study — the same used by Thompson (1982a). The instrument 
weighting functions are computed for each object profile 
rather than being computed for the current (or first) 
estimate. The noise covariance matrix for the HIRS was 
formed as a diagonal matrix using values supplied by L. 
HcMillin of NQAA/NFSS and are explicitly given by Thompson 
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(1982a) as HIRS2. Finally# in Fig. 2# th* eigenvalues were 
normalized by the total variance of the appropriate 
ensemble. In considering sn ensemble of temperature 
retrievals# it is more interesting to normalize on the total 
variance of the original RAOB information so that one can 
discuss the relative amount of original information which is 
explained by the ensemble of retrievals. Thus# if P is the 
total variance of a given ensemble of retrievals# and Q is 
the total variance of the original RAOB ensemble# then an 
adjustment factor F=P/Q may be applied to each eigenvalue to 
convert explained test ensemble variance to explained RAOB 
ensemble variance. It should be explained that a singular 
value decomposition is performed on each retrieval ensemble 
independent of the RAOB information# except that the 
covariance matrix contains variances and covariances about 
the RAOB mean rather than the mean of the given retrieval 
ensemble. The scaling performed by the multiplication by 
(P/Q) merely transforms the retrieval covariance matrix 
eigenvalues to numbers which reflect a percentage of total 
variance as though the total variance was that of the RAOB 
ensemble. The ensemble covariance matrix and the structure 
of the eigenfunctions# however# reflect the variability of 
the retrievals and not the variability of the RAOB 
information. Thus# to speak in terms of how much RAOB 
variance is "explained" by a given e i g enf unc t i on of the 
retrieval covariance matrix is. strictly speaking# an 
incorrect phrasing but will be used here anyway. 

5. Results 

RAOB data from the Phillips data set collected over ocean 
surfaces in January were partitioned into a tropical set 
(30S.30N) and a mid-latitude set <30N,60N>. Radiance 
measurements were synthesized and each of the retrieval 
algorithms were applied to the result. For this study# a 
priori statistical information was derived using one-half of 
each ensemble# and retrievals were performed on the 
independent remaining half. There were 200 mid-latitude 
profiles divided into dependent and independent sets of 100 
each# and 218 tropical profiles divided into sub-sets of 109 
each. As described in the previous section# a singular 
value decomposition was performed on each ensemble# and 
eigenvalues were F-adjusted. 

Fig. 3 shows a spectra of F-adjusted# cumulated 
eigenvalues of the information for the mid-latitude data. 
The RAOB curve shows a rather broad spectrum of structure 
information in which 99V. of the variance information is 
spread over about 15 orthogonal eigenfunctions. Radiances 
were converted to brightness temperatures using the Planck 
relationship# and the spectrum of this information is 
labeled BRIGHT. Since there are 65 levels of RAOB 
information but only 11 channels of radiance information# 
the total variance of the brightness temperature ensemble 
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was further adjusted by the factor ( 65 / 11 ) in order that 
results could be compared to RAOBs and retrievals. Comparing 
the RAOB and BRIGHT curves* one sees the substantial 
smoothing of information effected by the radiative transfer 
physics of the atmosphere. The integration of the detailed 
thermal structure of the atmosphere into 11 channels of 
radiometric information is* indeed* a drastic process. The 
brightness spectrum is assymptotic to the value 51.5 which 
means that the total variability of information in the 
radiance data is only about 51 . 5 % of the total variability 
of the information in the RAOB data. This is a significant 
reduction of information content. Temperature retrieval 
algorithms can extract no more information than this from 
the radiance data alone. The corollary is* of course* that 
any extra information in an ensemble of retrievals has to 
come from other sources* such as a 1 priori data. 

The curve labeled SMITH represents the spectrum of 
information in the retrieval ensemble using the Smith 
physical iterative method. As discussed earlier* the Smith 
method does not use a prio-i statistics as such but* as the 
other methods* does use a first guess profile* which is the 
RAOB mean profile here. Further* the transmi ttance 

information* which here has been computed exactly for each 
RAOB profile to be retrieved* is a form of a priori 
information and helps to expand the spectrum of information. 
Hence* the SMITH spectrum is assymptotic to 69. 8% of the 
total variability of RAOB information* which is somewhat 
more than just the brightness information alone. 

On the other end of the range of retrieval methods, 
regression retrieval results form the spectrum labeled 
REGRESSION in Fig. 3. The regression method does not use a 
priori estimates of atmospheric transmi t tanc es but does use 
a priori statistical information in the fullest possible 
way. The information spectrum is the broadest of any shown. 
In fact* in this case the regression retrievals exhibit 
slightly more total variance than the RAOB data. The 

difference between the^ REGRESS I ON and 3RIGHT spectra 
reflects the amount of a priori information added to the 
satellite measurement in producing the temperature 

retrieval. Evidently* this added information is a 

significant amount of the total information in those 

retrievals. For example* the total RAOB variance is 2340(K ) 
The adjusted total variance of the brightness temperatures 
is only about 51. 5% of that. Adding a priori RAOB 

statistics through the regression approach brings the total 
variance of the regression retrievals to 2350(ft* ). Thus* 
the a priori statistics imbedded in the regression method 
contributes about 1145(K*) to the total variance of the 

retrieval ensemble which means that the regression 

retrievals are about 51.3% satellite information and 48.7% a 
priori radiosonde information. One should keep this in mind 
when assessing the impact of "satellite temperature 
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retrievals'* on* say* numerical prediction models since the 
satellite data seems to have only a limited impact on 
"satellite temperature retrievals" depending on the method 
used. 

The information spectra curves for statistical 
regulari rati on and minimum information retrieval methods are 
also shown in Fig. 3. The minimum information algorithm 
yields a similar* although slightly broader information 
spectrum in comparison with the Smith method. The minimum 
information method is almost a pure physical method except 
for the mild influence of the a priori vertical mean 
temperature variance included ' n the retrieval matrix. Thus* 
the spectrum is assymptotic to a slightly larger percentage 
(73 0‘/.) of total RAOB variance than the Smith method. The 
statistical r eg u 1 ar i r a t i on which includes both the physics 
of radiative transfer and the covariance statistics of the a 
priori radiosonde information gives a broader spectrum than 
either the minimum information or Smith method* but one 
which is flatter than the regression results. The total 
variance of the ensemble of statistical r e g u 1 ar i i a t i on 
retrievals is 2085(K**> » about 89. 17. of the total variance 
of the original RAOB set. 

The stability and reliability of a retrieval method can 
be assessed by examining the error characteristics of the 
method. Fig. 4 shows vertical profiles of the root mean 
square retrieval errors for the results in Fig. 3. The 
error traces are similar with regression and statistical 
r eg u lar i z a t i on methods showing better accuracy than minimum 
information and Smith iterative methods. Thus* a 
combination of results* such as contained in Figs. 3 and 4* 
are necessary to make a reasonable assessment of the 
fidelity of information transfer. 

Figs 5(a)-(f) show the actual empirical orthogonal 
eigenfunctions of the various covariance matrices for the 
mid-latitude data. These figures clearly show the smoothing 
and distortion of the structure information due to the 
radiative transfer integration and the retrieval* or inverse 
step. The Smith and minimum information methods show the 
most smoothing of the vertical structure. Once again* these 
two methods make minimal use of a priori data from past 
radiosonde observations. Thus* the structure shown* even 
though drastically smoothed* exhibits the vertical thermal 
structure information contained in the radiance 
measurements. On the other hand. the regression and 
statistical regularization retrieval methods recovers more 
accurately the vertical structure of the thermal field* but 
clearly this is due to the incorporation of historical data 
on that structure. By the third or fourth eigenfunction* the 
structure of the original RAOB data is very poorly 
represented by any of the retrieval methods. A curious 
feature of these results should be pointed out. As the 
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spectra in Fig. 3 become flat* there is an ambiguity in 
ranking the eigenvalues. That is* successive eigenvalues are 
small and not significantly different from each other in 
this limit. Thus* objective ordering of the eigenvalues may 
actually lead to an incorrect ordering of the structure 
functions. In Fig. 5(d)* for example* the fourth structure 
function of the Smith retrievals is shown in comparison with 
the third structure function of the RAOBs. In Fig. 5(e)* the 
third Smith function and the fifth minimum information 
function bore the closest resemblance to the fourth RAOB 
function. In Fig. 5(f)* the sixth minimum information 
function is shown in comparison with the fifth RAOB 
function. These events are not particularly troublesome but 
indicate that care should be taken in truncating sets of 
eigenfunctions in representing atmospheric soundings. Of 
course* since degradation of the structure functions begins 
with the first function itself. one expects that it grows 
continually worse through higher order functions. 

Figs 6 and 7 show results for the mid-latitude case when 
random errors are added to the synthesized measurements. 
Errors were drawn randomly from a population whose standard 
deviations were equal to the expected noise levels of the 
channels of the HIRS instrument. These noise values were the 
current estimates for HIRS and were provided by L. HcMillin 
of N'QAA/NESS. For the results in Figs 6 and 7* the same 
error contaminated radiances were presented to each 
a 1 g or i thm. 

Now the variance of the set of random errors added to the 
synthetic measurements is significantly less than the 
variance* over the mid-latitude ensemble. of the synthetic 
measurements themselves. Thus* the impret of this noise on 
the information spectrum for brightness tem p era tur ef is 
rather minor. The error contaminated brightness temperatures 
have total variance of 1214(14** )» only slightly larger than 
the 1205(K*) for the unc ont am i na t ed measurements. Yet* the 
temperature retrieval algorithms respond to these errors in 
various ways. The regression* statistical r eg ul ar i za t i on* 
and minimum information variance spectra are flattened in 
this experiment in comparison with Fig. 3* while the Smith 
variance spectrum is slightly broadened. This result was 
unexpected for it implies that* except for the Smith method* 
the temperature retrieval algorithms actually 
overcompensated for the radiometer noise and produced 
retrievals with less overall variance than the retrievals 
using error-free measurements. On the other hand, the error 
profiles shown in Fig. 7 for these experiments clearly show 
a degradation of accuracy for retrievals using error 
contaminated radiances. Curiously, the Smith method behaved 
very well in this experiment with only a very minor 
degradation in overall accuracy. Thus. in analyzing the 
information transfer and distortion in sounding systems* it 
is of importance to first examine noise-free cases in order 
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to establish a baseline which involve* only tha radiative 
transfer integration and the retrieval method. 

It was also of interest to test the effect of radiometer 
noise level estimates on the information transfer. The 
statistical regular! ration# minimum information* and Smith 
methods were iterated until either the root mean square* 
instrument noise scaled deviations between estimated and 
actual radiances fell below 1* or until the rms temperature 
retrieval error reached a minimum. In addition* the 
statistical regular! ration and minimum information methods 
contain the radiometer noise covariance matrix explicitly in 
the inverse solution. Thus* assumed radiometer noise values 
will have an impact on the results. 

A series of experiments for the mid-latitude band was 
performed in which the radiometer noise values were 
increased to values approx imately equal to IX of spectral 
radiances produced by a mean atmospheric profile. (*) These 
noise values are used both in the noise covariance matrix 
and in the convergence criteria. Figs 8 and 9 show the 
results* for noise-free radiances* which should be compared 
with results in Figs 3 and 4. The RAOB* BRIGHT and 
REGRESSION results are not changed by the change in 
radiometer noise characteristics. The remaining three 
retries'al methods each produced smoother variance spectra 
with less total variance when the expected noise values were 
increased. 

Figs 10 and 11 show results for the tropical latitude 
tand using noise-free synthetic radiance measurements. The 
total RAOB variance for the tropical data set was 3810(K^-> 
as compared with 2340<K*) for the mid-latitude data. Notice 
also that the brightness temperature variance was a much 
smaller fraction of RAOB variance in the tropics: only 27. 3X 
or slightly more than one-fourth. Consequently* all 
retrieval methods in the tropics recovered smaller fractions 
of the RAOB information than in mi d-1 a t i t ud es because the 
information was more drastically filtered by the radiative 
transfer physics. A comparison of Fig. 11 with Fig. 4 shows 
that the tropical RAOB variance was smaller in the 
troposphere but larger in the stratosphere when compared 
with mid -latitudes. Generally speaking* the retrieval errors 
exhibit this same behavior with the regression retrievals 
showing the most reliability* even in the tropical 
stratosphere. 


<*) The noise values are rounded so that there is a single 
value for the 15i*.C02 channels and a single value for the 
4. 3u C02 channels which are* respectively* 0.63 and 0.0044 
mw/(m -sr-cnf >. 
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Fig* 12 and 13 show tropical results with synthesized 
radiances which have bran c ontami natad by random irrori. At 
in tha mid-latitude test# all mathodt except tha Smith 
method thowad a flight flattening of tha information 
spectrum and an incraata in ovarall rma arror. Tha Smith 
mathod thowad an incraata in ovarall varianca and handlad 
tha radianca arrort reasonably wall. 

Experiments in tha tropic* using tha largar "15i H 
instrument noita valuat behaved in a manner similar to those 
exparimantt in mid-latitude*. The result* are not shown 
hare. 

To summarize our experimental Table 1 shows values of tha 
lO-lOOOmb mean value of overall rms retrieval errors for 
each latitude bandi each simulated version of tha HIRS 
instrument# for both noise-free and noi s^fcontaminated 
measurements. As seen# tha regression retrieval method gave 
the best overall results. This is to be expected whan tha a 
priori statistics at hand are truly representative of tha 
profile* to b* retrieved. The Smith iterative method 
exhibited highly favorable results for all categories. This 
is a particularly important finding because that method 
depends least on a priori data and# hence# should suffer 
least if object profiles are not wall represented by a 
priori data. Tha two matrix inverse methods exhibit less 
stability to radiometer errors# although the statistical 
regularization method operated best for error-free 
measurements and "correct" a priori statistics. 

Table 2 shows what we shall call the impact of satellite 
data on satellite temperature retrievals. Tha entries are 
tha ratio of the (adjusted) total varianca of brightness 
temperatures to the total variance of various temperature 
profile ensembles. The brightness temperatures represent 
the information measured by a satellite device. As seen in 
Table 2# this information is only about one-half the RAOB 
i nf ormat i on in mi d -1 a t i t ud es. and only about oi. e-fourth tha 
RAOB information in the tropics. Retrieval methods using a 
lot of a priori information (Eg. regression) get much of 
their variance structure from this a priori information. 
Thus, the satellite data contributes only a lelatively small 
amount of the total retrieval information in this method. 
Physical methods (Eg. Smith iterative) are more faithful to 
the satellite measurement information - hence# the satellite 
data contributes a relatively large fraction of the total 
xe*-rieval information. In mid-latitudes# the regression and 

+ methods show least sensitivity to measurement noise 
uni le ali but the Smith method show mild sensitivity in the 
tropics. The Smith and minimum information methods would 
exhibit the least sensitivity to "incorrect" £ priori 
statistics. 


6. Conclusions and discussion 
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In this paper* we have had an eitensive look at the 
transfer and distortion of information concerning the 
atmosphere's thermal structure as this information takes on 
the different forms relevant to the satellite temperature 
retrieval problem. The contribution* here* has been to 
quantify this transfer process since most uiorkers in the 
field understand that distortion takes place. Our analysis 
method has been to perform singular value decompositions of 
the various ensembles of information and then to compare the 
variance spectra of each version of the information and to 
compare ensemble averaged errors between representations. 
Spec i f i c a 1 ly * the goal has been to quantitatively assess the 
distortion of atmospheric thermal information as it acts 
through the physics of radiative transfer to produce the 
outgoing radiation which is measured by a satellite 
radiometer* and the further distortion caused by various* 
classical temperature retrieval algorithms operating on 
these radiation measurements. 

We have shown that the process of radiative transfer in 
the atmosphere — being an integral process transforming 
atmospheric temperature into outgoing radiation — causes a 
drastic smoothing of information. That is* the total 
variance of outgoing radiation (rendered as a finite set of 
brightness temperatures) is half or less of the total 
variance of the atmospheric thermal field* even when values 
have been adjusted for the different numbers of measuring 
intervals for the two quantities. 

Temperature sounding systems which produce retrievals 
with a large amount of variance. and with greater apparent 
vertical resolution than simple profiles of brightness 
temperature* generally do so by adding a priori information 
(virtual measurements in the language of Rodgers (1976)). 
taken usually from past observations of the temperature 
profile. or from a numerical forecast model. While passing 
no judgement on the usefulness of this technique in any 
given case. we have attempted to quantify how much a priori 
information is added in each of several types of temperature 
retrieval scheme. Thus, for mid-latitude retrievals using a 
statistical correlation-regression technique* with 
r ep r e s eri t a t i ve i priori statistics* about 4S-49X of the 
total variance in an ensemble of retrievals is attributable 
to the a priori data. The total variance of the regression 
retrievals themselves is about the same as the total 
variance of the atmospheric thermal field producing the 
radiances upon which the retrievals are based. Using a 
physical retrieval technique. such as the iterative method 
of Smith (1970). about 73-75X of the total variance of the 
retrievals is attributable to radiance information. while 
on iy 25-277. is attributable to a priori information such as 
atmospheric transmi ttance functions. However. the total 
variance of these retrievals is only about 70V. of the total 
variance of the original thermal field. Retrieval algorithms 
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such *i minimum information* or statistical regular i zat ion 
product result* between these two extremes* and reflect the 
effective mil of physics and a priori statistics of those 
methods. 

In examining the average temperature retrieval errors* we 
find results consistent with other investigators. Namely* of 
the four retrieval methods tested* the ranking by overall 
accuracy and stability to measurement error* is regression 
(highest)* Smith iterative* statistical regular i tat ion* 
minimum information. 

The vertical profiles of the empirical orthogonal 
eigenfunctions of each retrieval ensemble were compared 
against the same for the RAO0 data defining the original 
information. The smoothing and distortion of these functions 
has been specifically illustrated. Generally speaking* we 
have shown that fidelity of the vertical structure of 
satellite retrieved temperature profiles is rather low for 
physical retrieval method*. The increase in fidelity for 
physical -statistical or pure statistical retrieval* is 
attributable* therefore* to the inclusion of a priori data. 
As suggested by the work of Crosby and Weinreb ( 1974 ) and of 
Spankuch. Timofeyev and Gtildner ( 1977 ). the validity of this 
process of "fleshing out" the vertical structure in 
satellite temperature retrievals is only as good as the 
correctness of the atmospheric statistics for each object 
retr i eval. 

We end this paper with a few philosophical remarks. There 
is* and should continue to be* a question concerning the 
impact of satellite observing systems on the analysis and 
prediction of atmospheric behavior. The Data Systems Test 
presented our first comprehensive look at the impact on the 
large scale numerical prediction problem. This test has 
produced a mixture of results (See Tracton* Desmaris* 
VanHarren* and McPherson (1980)* Ohil* Halem and Atlas 
(1979)) which indicates that the i inpact is a complicated 
function of initialization and assimilation procedures and 
the numerical forecast model itself (Atlas* et al ( 1981 )). 
Our results here show that even the question of the impact 
of satellite data on satellite temperature soundings is one 
which produdes a mixture of results depending on the type of 
retrieval algorithm used and the representativeness of a 
priori data. Eecause of the nature of the ata set used 
here, we were not able to simulate the i ni t i a 1 i z a t i on- 
assimilation step to measure the further smoothing of the 
original measurements in preparation for the large scale 
problem. We strongly recommend that this be done. The 
results shown here indicate that s ta t i s t i c a 1 1 y optimized 
satellite temperature retrievals ar if very heavily influenced 
by the statistics so that impact te vH may be assessing the 
effect of a priori statistics on wo c.e r forecasting as much 
as they are assessing the impact of satellite observations 
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on forecasting. It would be interesting to conduct an impact 
test using* sag* a physical method such as the Smith 
iterative method. Alternatively* one could seek a wag to 
conduct the assimilation step in r ad i anc e space instead of 
temperature space* perhaps avoiding the inverse problem 
il together. 

We expect that these issues will become even more pointed 
as onb delves deeper into mesoscale applications of 
satellite soundings. To achieve mesoscale observations from 
satellites in the horiiontal or temporal sense suco as 
through VAS* (Smith* Suomi* et al (1981)) is significantly 
different than to achieve the correspondingly relevant 
v er tical. scales for these me teor o 1 og i cal systems. 
Presumably. the interesting vertical structure of mesoscale 
systems will have to be gotten from 5 priori data* from 
independent observations such as discussed by Westwater and 
Grody (1980). or. possibly. from a significantly different 
incorporation of radiosonde observations (Fritz (1977)). To 
the extent that atmospheric mesoscale systems are not merely 
miniatures of the global scale system* we suggest that 
global scale satellite sounding procedures are not 
necessarily the optimum techniques for sounding mesoscale 
weather systems. 
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TABLE 1. lO-lOOOmb mean of root mean square retrieval 
errors. 


30N-60N 

RAOB MEAN 

REGRESSION 

STAT. REO. 

MIN. INF. 

SMITH IT. 

HIRS2 

no error 

6. 042(K) 

1. 730 

1. 625 

2. 882 

2. 702 

HIRS2 

errors 

— 

2. 073 

3 517 

4. 047 

2. 703 

HIRS1 

no error 

6. 042 

1. 730 

2. 208 

2. 747 

3. 049 

HIRS1 

errors 

— 

2. 440 

533* 

5- 3ee. 

z. 

30S-30N 

RAOB MEAN 

REGRESSION 

STAT. REG. 

MIN. INF. 

SMITH IT. 

HIRS2 

no error 

4. 885 

1. 627 

2. 896 

3. 656 

2. 703 

HIRS2 

errors 

— 

1. 679 

3. 141 

3. 896 

2. 647 

HIRS1 

no error 

4. 885 

1. 627 

2. 006 

3. 190 

3. 088 

HIRS1 

errors 

— 

2. 176 

z s-yf 

3- 3?* 

3- z'? 


TABLE 2. The impact of satellite data on satellite 
temperature soundings. Ratio of adjusted total variance of 
satellite brightness temperature ensemble to total variances 
of other temperature profile ensembles. 




C TOTAL 

VARIANCE 

X 3/C TOTAL VARIANCE 

Y3 



MID-LATITUDES 

TROPICS 




30N - 60N 

30S - 3 ON 


CX/Y3 


No Errors 

Errors 

No Errors 

Errors 

BRIGHT/RAOB 


51. 47. 

! * 

I 

1 to 

1 

1 

27. 37. 

27. 17. 

BRIGHT/REGRESSION 

51. 4 

53. 5 

28. 8 

28. 7 

ERIGHT/STAT. 

REO. 

57. 7 

76. 2 

36. 7 

38. 4 

SR IGHT/MIN. 

INF. 

70. 5 

96. 8 

64. 1 

65. 0 

DR IGHT/SMI TH 

IT. 

73. 8 

73. 9 

84. 3 

66. 1 
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Figure 1. The f 1 ou of information through a satellite 
temperature retrieval and analysis system. Ovals represent 
different versions of information concerning the atmospheric 
thermal structure. 
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figure 2. Singular value decomposition of temperature and 
corresponding satallite measured radiance information. Shown 
is the spectrum of cumulative percent explained variance in 
the empirical orthogonal functional representation of the 
information of each type. 
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Figure 3. Cumulative percent explained RA03 variance for 
various satellite temperature retrieval methods in mid- 
latitudes (30Ni60N) for January. Also shown are RAOB and 
satellite brightness temperature information spectra. 
Radiance measurements are synthesized for the H1RS sounder 
with no errors. (See text for explanation. ) 
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ERROR PROFILES 
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Figure 4 . 
retrieval 


Vertical profiles of root mean square 
errors for mid-latitude cases included 


temp erature 
in Fig. 3. 


*c. 


KIGKNI- UNCTION - I 


EIGENFUNCTION • 1 


Dr.GRF.ES KELVIN 



EIGENFUNCTION »3 

EIGENFUNCTION «4 

KIGKNFL'NCTION -5 


Figures 5(a) — (f). Mean profile and first five empirical 
orthogonal functions of RAOS temperature profiles and 
corresponding satellite temperature retrieval ensembles in 
mid-latitudes correspond ing to Figs. 3 and 4. 
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Figure 6. Information spectra for mid-latitudes as in Fig. 3 
except using error contaminated HIRS synthetic radiance 
mea s ur erne n t s . 



ORIGINAL PASf J* 
OF POOR QUALITY 


78 


30 S LATITUDE -S60 



Figure 8. Information spectra for m i d- 1 a t i tud es» as in Fig. 
3/ except using radiometer noise values equal to 
approximately IV. of standard radiances. 
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igure 9. Vertical profiles 
n mi d - 1 a t i t ud e s for the 


of temperature retrieval 
experiments in Fig. 8. 
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10. Same as Fig. 
C !) in January. 


tropical latitudes 
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Figure li. 
(30S, 30N) 


30S-30N, MIRS1 , NO ERRORS 
ERROR PROFILES 



RMS TEMPERATURE DIFFERENCE 


Same as Fig. 4 except for tropical latitudes 
in January. 
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Figure 12. Samt ts Fig. 10 except using error contaminated 
synthetic HIR5 radiance measurements. 
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Figure 13. Same as Fig. 11 except using error contaminated 
synthetic HIRS radiance measurements • - 
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V. Satellite Temperature Soundings in Cyclones 
(Microwave Satellite Observations) 

by 


Sigmn- 
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ABSTRACT 


Satellite T1R0S-N microwave radiances were used to derive 
temperatures and thickness over a cyclone in the Gulf of Alaska. 

The radiances were used in conjunction with the closest radiosonde 
data to the cyclone. 

The method employed solves three equations in three unknowns. 
The unknowns are the three coefficients which are needed at each 
pressure level, to multiply the measured radiances in order to 
retrieve the vertical temperature distribution. From the derived 
vertical temperature distributions, the thickness between 700 and 
300 nib were computed. Maps of the thicknesses were compared with 
NMC maps and with thicknesses derived from NESS operational sound- 
ings. 

From a few soundings the tropopause height is seen to have a 
minimum near the Gulf of Alaska . /clone. 
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Satellite Temperature Soundings in Cyclones 
(Microwave Satellite Observations) 
by Sigmund Frit* 


1 . Ijil r o duc t i on 

Cyclonic areas over the Pacific Ocean are among the most 
significant regions for weather analysis and forecasting for 
North America. 

However, in. such important regions, it is difficult to make 
accurate temperature soundings with infra-red satellite measure- 
ments because clouds greatly influence the emission to space. Con- 
sequently, satellite measurements do not represent the effect of 
atmospheric temperatures alone and it becomes difficult to extract 
the atmospheric temperature from the satellite observations. 

We have therefore restricted our investigation to the use of 
microwave satellite observations. At some frequences microwave 
observations are usually hardly influenced by the presence of 
clouds. Nevertheless, recently, Phillips (v80) has raised a caution- 
ary flag regarding the influence of precipitation areas on the measure- 
ed microwave data; we shall address this problem later. 

On April 6, 1979 a cylcone was located over the Gulf of Alaska 
(Fig. 1) and TIROS-N microwave observations were available. From 
the observed radiances, temperature soundings tan be derived by many 
methods. This case was also studied by Brodrick (1980). Brodrick used 
NESS operational soundings which depend on certain parameters; these 
parameters are derived from statistical relations between satellite 
soundings and radiosondes which are mainly located far from the Gulf 
of AlacVj, and also proceed the April 6, 1979 date by a week or more. 

Yet according to Grody (1980) "temperature profiles, having large 
amplitude small scale fluctuations relative to the a priori mean, 
result in the largest retrieval errors". This leads to his further 
conclusion that the "a-priori" vertical temperature sounding should 
be based on a particular synoptic situation rather than on "a-priori" 
statistics according to three basic latitude bands. It is partly 
for that reason that Fritz (1976), and Hillger and Vender Kaar(1977,1979) 
used radiosonde data together with satellite data to derive vertical 
temperature structure in specific synoptic situations. 
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The method we have selected does not depend on any statistics. 
Instead the available radiosondes closest in space and time to the 
cyclone are accepted as basic true imput data. The satellite data 
are then used to interpolate the temperature soundings between the 
radiosonde stations. 

In this method, there is no error at the radiosonde stations 
because we have assumed that the radiosondes are perfectly correct. 
The retrieval results can be compared everywhere with the results 
from operational methods discussed by Brodrick (1980), and with 
the data used in NMC analyses. 

2. Thejlethoti 

Solving 3 equations in 3 unknowns^ The method used here has 
been described earlier (Fritz, 1976, 1977(a), (b) , 1978). The 

method assumes as do most linear methods, that the temperature, at 
any height or pressure level in the atmopshere, can De expressed as 
a linear combination of satellite observations. 

Thus, we compute a ‘ emperature deviation, AT, from the re- 

J ation 

AT^(k) =* E F(k , j) • AR(j) (1) 

and find the temperature, T, from 

T (k) - T (k) +AT (k) (2) 

r S 

Here T is temperature, R is radiance measured by the satellite, and 
the F' s are a set of coefficients. The subscript, r , means the 
retrieved value, "g" means the "guess" or initial value. The letter, 
"k" refers to a pressure level in the atmosphere, and "j" refers to 
the spectral f i equency of the radiation measured. 

A R( j) = R m (j)-R g (j) 

where the subserpt "m" means the measured value. 

The problem is to select the "best" values of F(k,j). 


( 3 ) 


88 


In general the values of F(k,J) depend on the vertical temper- 
ature structure of the atmospheres from which they are derived. Thus, 

In statistical methods the F's depend on the lapse rates of the sample 
of soundings which are used In the statistical sample. However, if 
the F* 8 could be related to vertical temperature distlbution In the 
cyclone under investigation, the temperature retrievals might be more 
accurate than the retrievals which depend upon a relatively unrelated 
statistical sample. 

Therefore, in the study we selected the nearest radiosondes avail- 
able. Figure 1 shows a picture of the cloud distribution over the 
cyclone; the dots locate the radiosonde stations (Table 1) used to 
derive the F's. The cyclone center was located on the surface synoptic 
map at 1200Z at about 54N, 167W. The center of the cloud spiral seems 
to be at about 53N, 165W. 

The Microwave Sounding Unit (MSU) on TIROS-N has four channels 
Those are: 


Channel No. 

1 

2 

3 

4 


Microwave Sounding, Unit Frequencies 
Frequency (GHz) 

50.30 

53.74 

54.96 

57.94 


The radiance in channel 1 is strongly influenced by clouds. Therefore 
channel 1 was not used to sound the atmosphere. Channel 2 may be affected 
somewhat by clouds where rain is occuring. (Phillips, 1980). It may have 
been raining at Ship Papa, when TIROS-N viewed the area, however we shall 
show that channel 2, radiance did not seem noticeably affected in this case. 


The remaining channels, channels 3 and 4, are presumably affected 
only by changes in atmospheric temperature. 


We therefore have radiances available from 3 MSU channels namely, 
channels 2,3, and 4. 


The method used here is therefore limited to the solution of 3 
equations in 3 unknowns. The equations are 3 equations like Eq.(l). 

To generate the AT(k) at 3 radiosonde stations, we select one radio- 
sonde station as the "first guess -tation. We selected station //70308 
(St. Paul, Alaska). We then selected 3 additional stations "surrounding" 
the cyclone. These three stations are listed in Table 1. 
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Table 1 - "first Guess' 1 and dependent radiosonde s tations . 


Station 

Station name 

Lat. 

Long. 

70308 

St Paul 

57.15N 

170. 22W 

70454 

Adak 

51.88 

176.65 

70316 

Cold Bay 

55.20 

162.72 

C7P 

Ship Papa 

50.00 

145.00 


(C) 


(G) =* "First Guess" Station, It 70308 

TIROS-N viewed the ar>a at about 1 320Z and 1500Z on April 6, 1979 
the observed radiances over the four stations (Table 2) were selected 
from maps of MSU radiances for each of the channels. No adjustment 
was made for the difference in time from the radiosonde data, nominally 
at 12 00Z • 

o . 

Table 2 -MSU radiances over Radiosonde stations. ( K) 


Station 


Channel N 


2 

3 

70308 

245.5 

229.5 

70454 

244.0 

228.5 

70316 

245. 3 

230.6 

C7P 

2^8.7 

225.5 


4 

222.0 
222.1 
221.5 
215. 3 


(G) 


(G) 3 "First Guess" stations 


Table 3 -- Values of AR(°K) 


Stall on 


(Oh^U 

70316 

CiP 


Channel K umber 


2 

■ 1.5 

- 0.2 

3-2 


,_3 

- 1.0 
1.1 
-U . 0 


h 

0.1 

-0.5 

-6.7 
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The radiance values, expressed In equivalent degrees K, were 
subtracted from the values at the "First Guess" station , to obtain 
the values of AR(j) shown in Table 3. 

The radiosonde temperature data were available at standard levels 
up to 100 mb, except that at Ship "Papa", the data were given to 30 mbs. 
In order to use a program already available, fictitious temperatures 
were added at the 5 mb and 0.1 mb level. The radiosonde temperatures 
were then interpolated to 100 pressure levels from 1000 mbs to 0.01 mb 
level according to a p2/’ law. 

This Interpolation retains the measured temperature values at the stand- 
ard levels exactly. Tncidently, the surface pressure was generally low- 
er than 100 nibs;e.g.,it was 967 mb at station #70316. 


Having now obtained temperatures at 100 pressure levels (actually 
only those at 50 levels, from 1000 to 104 mbs were meaningful) it was 
a simple matter to solve the 3 equations like Eq(l), separately at each 
pressure level, k, to obtain 3 values of F(k,j) at each pressure level. 

An example of the distributions of F(k,j), is given in Fig. 2. 

Below about 400 mbs, the coefficients are fairly small never exceeding 
an absolute value of 2. However wear the 300 mb level , the coefficient 
for channel 3 reaches a value of 1-5. This may be too large to give stable 
values of retrieved temperature for independent radiance measurements, if 
substantial random errors exist in the radiance observations. For that 
pressure level, errors fit the radiance measurements would be amplified 
more because they would be multiplied by the large values of F(70,3); 
(k-/0 corresponds to 299 mb). However, when we derive the thickness 
over the layer 700-300 mbs, as we shall do later, the coef f iecients are 
quite small when averaged over that thickness layer. 


In spite of the variations of the F(k,j) in Fig. 2., obviously 
when the F(k,j) are multiplied by the AR(j) at the radiosonde stations 

(Table 3), and then added to "guess" temperatures, we obtain exactly 
correct temperatures at all levels at each radiosonde station. With 
this method there are no errors in the temperature at the dependent 
radiosonde stations. Since there are no errors in the temperatures, 
there are no errors in the heights of the pressure surfaces or of the 
layer thicknesses at those dependent radiosonde stations. 

3. Application to independent data. 


Having obtained values of F(k,j), we are now ready to apply them 
to all the TIROS-N MSU measurements over the cyclone. We derive temper- 
a t u res at all 50 levels from 1000 mbs to 100 mbs. From those derived 
temperatures,^ we can compute the heights of pressure surfaces, and th - . 
compute the thickness between any two pressure levels. 
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This was done for an area larger than the region surrounding tha 
flulf of Alaska cyclone (Fig. 1). For that, data from two T1R0S-N orbita 
were used. The times of those observations were approximately 1320Z and 
1500Z. 


The thickness of the 700-300 nib layer is shown in Figures 3 (1320Z) 
and Fig. 4. (1500Z). The dots in Fig. 3 show the locations of the HSU 
observations at about 1320Z; in Fig. 4., the dots show the locations at 
about 1500Z. Along the latitude 50N, the thickness increases from 175E 
eastward. The thickness is about 5850 m at 175E, about 6000 o at 156W, 
and reaches a maximum of 6230 ra, at 140W. At 53N, 170W, Just west of the 
cyclone center, as shown in Fig. 1. , the thickness gradient is quite small/ 
in fact a low thickness of 5925 ra was located there, surrounded by some- • 
what larger values of thickness. 

How do these values compare with other results? Mr. H Brodrick of 
NESS, kindly supplied values of 700-300 mb thicknesses derived from 
"operational" retrievals. Unfortunately those retrievals where not made at 
each MSU observation spot. So Brodrick's values were interpolated between 
the thickness shown in Figs. 3., and 4. NMC's thicknesses, also supplied 
by Brodrick, were also compared in a similar manner. 


Comparisons with NMC (thickness 700-300 mb) and between NESS and 3 Eqs 
retrievals 


It is interesting to compare the 700 -30C mb thickness retrievals 
with NMC's thicknesses. This is done in figs. 5(a) and 5(b), for 1 500Z 
observations, and in fig. 6(a) and 6(b) for 1320Z observations, bet s 
call tire methods described above the "3-F.qs" method. Fig. 5(a) compares 
the TIROS -N 1 500Z observation results of the "3 Eqs" method with NMC; 

Fig. 5(b) compares the results of the NESS statistical method with NMC 
(Brodrick, 1980). The largest differences from NMC occur near 45N* 155W. 
Near lhat location, both comparisons would reduce the NMC thicknesses. 

But the "3-F.qs" method would reduce the thickness by about 80 meters, 
whereas the NESS method would reduce the NMC values by about 150 meters. 
Thus while both methods would intensify the NMC trough at Lhat location 
the "3-F.qs" method would do so less drastically. That result is doubtless 
influenced hj the utilization of the sounding at Ship Papa (50N, 145k) In 

the "3-Fqs" method. 


Figs. 5(a) and 5(b) are based on comparisons with interpolations by 
NMC and NFSS to the NMC polar-sterographic grid-point locations. Moreover, 
in these figures, the TIROS -N observations were made near 1500Z. Mr. Brodrick 

kindly supplied also the NESS retrievals at specific locations. 

These are NESS retrievals made over an average of several satellite obser- 
vation : therefore the location of the NESS retrievals, unfortunately, also 
do not correspond exactly to the retrievals at the MSU observation locations 
indicated by the dot in Figs. 5. Nevertheless the NESS retrievals were 
compared with the "3 -Eqs" retrievals; the comparisons are shown in 5(c). 

ORIGINAL PAGE IS 
OF POOR QUALITY 
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Fig. 5(c) shows that differences of about +100 10 were computed near 
155W from about 49N to 56N. Elsewhere the differences were smaller; near 
50N, 170W the largest difference was -AO m. Doubtless some of the val- 
ues in Fig. 5(c) are Influenced by the interpolation made necessary be- 
cause the retrievals were not made at exactly the same locations; but the 
patterns of positive differences near longitude 155W, and negative value 
near 170W is probably valid. This would of course change the gradient 
of thickness and, therefore, the thermal wind. 

Another factor which may have affected the results, is the fact 
that 45N , 1 55W is near the satellite viewing horizon. The observed 
radiances are made at a large slant angle there; corrections are made 
for limb effects in an attempt to calculate what the radiances would 
have been if the satellite observations had been ir. the nadir. It was 
these . itnb-cor rocted" radiances, kindly supplied by Dr. Grody, which 
were used in the retrievals with the "3 Eqs" method. Perhaps some vari- 
ation in the "lint-corrections" made in the operational retrievals may 
account for some of the differences in the NESS and "3-Eqs" results. Of 
course near the center of Fig. 5 and Figs. 6, the satellite observes 
near nadir, and the differences in Figs. 5 and 6 would be due to other 
factors. 


TIROS-N, 1320Z observations 

Tn Figs. 6(a) ind 6(b), the retrievals Horn the "3 F.qs" and from 
the NESS methods are compared both NMC for satellite observations made 
ut about 1320Z. Near 50N, 155W, the results are similar to the ones in 
Figs. 5(a) and 5(b) at that location. In Fig. 6(a) the results are some- 
what different than in Fig. 5(a) because no adjustment was made to make 
the retrievals equal at both observation times. Tn Fig. 6(a) near 55N, 

135W, the "3-Eqs" method differs from NMC by about 100 rn whereas the 
NESS retrievals agree with NMC almost exactly. At 57N, 1 30W NESS retrievals 
are about 28 m larger than NMC, while the "3-Eqs" method shows about 
+70 m. The reason for these discrepancies are not obvious, and should 
he investigated fuither. However, it should be noted that in an opera- 
tional application o' the "3 Eqs" method, a radiosonde (Annette Island), 
near 56N, 132W would have been utilized to retrieve the correct thickness 
at that location. 

The difference between the "3 Eqs" and NESS retrievals for 1320Z, 
is shown in Fig. 6(c). Here again we note the relatively large differ- 
ences near the horizon. Near the center of Fig. 5(c), that is, near 
satellite nadir, the differences are smaller although they still reach 
about +50 m. 
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Which result Is better? It does not seem possible to decide from the 
thickness field alone. If the radiosonde measurements are accurate, 
obviously the results of Fig 3 and 4 must be better at the radiosonde 
stations. But the NESS results are also in approximate agreement with 
the Aleutian Is. radiosonde stations. Near Ship Papa the NESS results 
are not in such good agreement. 

Thicknesr was also derived tor the layer 1000-500 mb, but the re- 
sults are not r.eported here. The average values of the coefficients 
which yield average temperatures over the layers are shown in Table 4. 


Table 4 

Average coefficients, F(k,J) averaged over pressure 
Pressure Intervals Channel No. 


(mbs) 

1000-500 

700-300 


2 

3 

1.16 

-0.48. 

+0.86 

+0.01 


4 


- 0.21 


-0.99 


Intervals 


4. Tropopause 

Since the retrievals yield temperatures at 50 pressure levels 
between 1000 and 100 mbs, it was interesting to see what the vertical 
temperature structure was. 

It should be remembered that the vertical temperature structure 
at the radiosonde stations (Table 1) are accepted as true. Therefore, 
the tropopause at these stations, and all the temperatures at every 
level are assumed known. The satellite measurements merely serve a3 
a rational way of interpolating between the radiosonde stations. 

With that in mind, we see in Fig. 7 two radiosonde temperatures, 
and also a set of retrievals at five locations along a line irom Ship 
Papa to Adak (// 70454). It is interesting to note that the lowest trop- 
opause height occurs just south of station II 70316 which itself had the 
lowest tropopau'-'e height of the four radiosonde stations; namely at 394 mb. 

The lowest tropopause height judging by Fig. 7, seems to have 
occurred just to easjt of the low level cyclone. 


• *4 
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5 • On ths question of the value of channel 2 radiance at Ship Papa 


Phillips (1980) has shown that when precipitation is cccuring in 
the atmosphere under the satellite, the observed radiance in channel 2 
may be too low. Under such circumstances the cold cloud is partially 
substituted for the warmer atmosphere which lies below. 


In order to detect situations when precipitation is occurlng, 
Phillips suggested comparing channel 1 with channel 2. Without pre- 
cipitation, over the ocean, the radiance R. in channel 1, tends to 
be much smaller than the radiance R 2> in channel 2, but when precipi- 
tation Is occuring, R^ approaches R* in value, and R^ may even exceed 

r 2 . 


Fig. 1 shows that Ship Papa was located under a frontal band; 
it may have been raining there. Is the radiance R^, observed over 
Ship Papa too low? Should the radiance be increased somewhat? To 
test this a map of R 2 - R^ was prepared; the data arc shown in Figure 8. 


The value of (R 2 “ Rj) at Ship Papa is about 10; Phillips suggested 
that when (R 2 - Rp is less than 12, then R^ might be affected by the 
clouds. To test this further, and to determine by how much R- at Ship 
Papa should be increased the observed values of Rg where compared with 
theoretically computed values. 

The value of R 2 iB' given by mainly 


R 2 " J 0 T dT 2 * (4) 

In Eq.A, "T" is the temperature obtained from radiosondes and 
t 2 is the atmospheric transmittance for channel 2. Values of t 2 
were kindly supplied by Dr. N. Grody. 

The radiosonde data were mostly available up to 100 mbs although 
some stations supplied temperatures to greater heights. The surface 
pressure at most stations was lower than 1000 mbs. Equation (A) was 
therefore applied for pressure levels from 970 mbs to 100 mbs at 
several stations. The results are shown in Figure 9. From Fig. 9 it 
appears that no correction should be applied to the observed value 
of R_ at Ship Papa. A line fitted by eye through the Hawaiian stations 
and the Alaska stations passes through the point for Ship Papa. Another 
line which takes account of Midway Island, would Indicate that the 
radiance at Ship Papa should be r educ ed, instead of increased. 

Because of this analysis, the thickness shown in Fig. 3 and A 
were computed with the observed value of R 2 at Ship Papa. Computations 
were also made with the value increased from the observed value of 


0<L 


•fc * 
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of 248.7 K to i fictitious value of 250.5 °K. Of course, because 
of the method of solving 3 equations in 3 unknowns to derive F(k,J), 
no errors appear in the temperatures of the radiosonde stations, with 
either value of R^. It is not until the retrievals are made at a sub- 
stantial distance south of the radiosonde stations that a significant 
difference appears in the 700-300 mb thickness between the retrievals 
with R 2 * 248.7 and R 2 =“250.5. 

The coefficients F (k,j) for the average temperature from 700- 
300 nibs were: 


Channel No. 


R 2 - 248.7 

L. 

0.86 

0.01 

R 2 = 250.5 

0.82 

0.08 


4 

-0.99 


-0.83 


Since AR is small near the Aleutian Island stations, and the 
temperature errors are zero at the all dependent radiosonde stations, 
it is only when the AR's change appreciably from the value at the 
radiosonde stations thai significant changes in the thickness can 
occur. 


6. Further Discussion 


As Grody (1980) suggests, the details of the vertical temperature 
structure may introduce errors in temperature retrievals from satellite 
data. Fritz (1969) also indicated that true deviations from assumed 
average temperatures make it impossible to derive error free retrievals 
at a given height from satellite data alone. Grody suggests that the 
vertical temperature structure in every synoptic situation needs to be 
considered separately. 

Retrievals, using both radiosondes and satellite radiances, attempt 
to employ both types of data as a consistent data set. There are however 
questions raised when such methods are used. The present method depends on 
radiosonde data. If any radiosonde sounding is in error, that error 

will be propogated into the retrievals in the vicinity of the erroneous 
radiosonde. Moreover, if the satellite data are not made at the same time 
as the radiosonde, in a changing synoptic situation the two data sets 
would not be consistent. Therefore, the two sets of data should really 
be colocated in space and in time, although no attempt was made in the 
present study to account for the time differences between the radiosonde 
ascents and the satel .te observations. The radiosonde data are taken, 
nominally, at 1200Z, although some allowance should be made for time 
taken by the ascent and deviations of the release time from 1200Z. The 
satellite data were taken in 1320Z, in Fig. 3 and near 1500Z in Fig. 4. 
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Doubtless, somewhat different retrieval results would have been 
obtained if the radiosonde data had been adjusted to the satellite times, 
by for example, using forecast temperatures made at hourly intervals. 

In the course of the work, it wan noted that the satellite radi- 
ances were substantially different during local daytime from the values 
during local night-time. If this is due to some Instrumental error in 
the satellite, obv'ously the coefficients derived for 1200Z would not 
be applicable to any other times. Space Interpolation might still be 
permissible, but extrapolation in time would not be valid. This day- 
night variation in the MSU data needs to be investigated further. 

Another item which needs to be looked into further, is the reason 
for the 100 m discrepancy in Fig 6(a) between the thicknesses derived 
by the "3-Eqs" method and NMC's result, in the vicirity of Annette Island. 

Further test's of the "3-Eqs" method can readily be devised. For 
example it wc.dd be interesting to substitute another station for Ship 
Papa, and to 'ompute thickness retrievals again. Instead of Ship Papa 
we might employ a radiosonde in the Hawaiian Islands, or another one 
of the Alaska coastal stations. 
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Flguro Legends 


Figure 1 


Figure 2 : 


Figure 3 : 


Figure A 


Figure 5(a): 


Figure 5(b): 


Figure 5(c): 


Figure 6(a): 


Figure 6(b): 
Figure 6(c): 


Picture of spiral cloud band over the Gulf of Alaska from 
TIROS-N, April 6, 1979. The heavy line, surrounding the 
spiral , encloses the area emphasized in the report. The dots 
in the enclosed area show the locations of the radiosonde 
stations used as dependent stations. "Ship Papa" is denoted 
by "P," the "first guess" station by "G." Midway Island is 
labelled "M." Ship Papa appears to be under n frontal cloud 
band. 

Vertical distributions of the coefficients F(k,J), for 3 
microwave channels. The tick marks on the pressure scale 
are for every 5 units of "k." k * 50 refers to p » 97 mb, 
k ■ 100 refers to p - 1000 mb. 

Thickness (meters) of 700-300 mb layer, computed on the 
basis of TIROS-N microwave radiance measurements and radio- 
sonde data. Microwave measurements were made at about 
13202. April 6, 1979. The dots show the locations 
of the microwave measurements; temperature and thickness 
retrievals were computed at each dot. The symbol at 50N, 
1A5W shows the location of Ship Papa. 

Same as Figure 3, except for microwave measurements made at 
about 150GZ. The symbols denote the locations of the 
dependent radiosonde stations. The storm symbol denotes 
the approximate center of the cloud spiral In Figure 1. 

Comparison with National Meteorological Center (NMC) 
thicknesses (meters). Thickness from Figure A (1500Z) 
minus thickness from NMC. Symbols denote locations of 
dependent radiosonde stations. Note about -80 m 
difference near A5N, 155W. 

Same as Figure 5(a) except difference between NESS 
operational thickness and NMC thickness. Note about 
-150 m near <*5N, 155W. 

Comparison of thickness (meters) in Figure A with 
thickness from NESS operational soundings; Figure A 
thickness minus NESS thickness. Symbols show locations of 
dependent radiosonde stations. 

Same as Figure 5(a), except for 1320Z TIROS-N observations 
(Figure 3). 

Same as Figure 5(b), except for area further East. 

Same as Figure 5(c), except for 1320Z TIROS-N observations 
(Figure 3). 
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Figure 7 : Temperature soundings. Soundings no. 1 end no. 7 are observed 

from radiosonde stations. Soundings no. 2 through no. 6 are 
retrievals from the "3-eqs." method, along latitude 50-51N, 
between the two radiosonde stations. Note the low tropopause 
near the center of the longitude interval. 

Figure 8 : Radiance, R^, minus radiance, R. , for T1R0S-N observations at 

1320Z. Note tight gradient along coastal areas. 

Figure 9 i Comparison of partial radiances computed from radiosonde 

temperature (eq. 4) versus MSU channel 2 radiances observed 
at the same radiosonde stations. Values "First Guess" 
station (computed or observed) subtracted from corresponding 
station value. 
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1 QUILLAYUTE (72797) 

2 ANNETTE (70398) 

3 YAKUTAT (70361) 

4 FAIRBANKS (70261) 

5 ANCHORAGE (70273) 

6 KODIAK (70350) 

7 McGRATH (7023l) 

8 BETHEL (70219) 

9 COLD BAY (70316) 

10 C7P 

SI ADAK (70454) 

12 MIDWAY IS. 

13 LIHUO (91165) 

14 HILO (91285) 
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VI. Northern Hemisphere 700 mb Heights and 
Pacific Ocean Temperatures 
in Winter 

by 

Sigmund Fritz 


Abstract 


Horel and Wallace (1980) show the time variations of Pacific Ocean 
sea-surface temperature (SST) patterns for the winter months (average 
for December, January, and February). The 10 years with the warmest 
equatorial Pacific SST and the 8 years with the coldest SST were selected 
from their date. 

Northern hemisphere maps of the 700 mb heights, everagea i or the 10 
SST warm years, were computed. This was done for the mean monthly height 
at each 10 latitude-longitude intersection; it wa« done separately for 
November, December, January, and February. Similar maps were prepared 
for the average of the 8 SST cold years. Then the difference in the 2 
averaged height fields was computed. Maps for the difference are shown 
as Figures 1-4 in this report. 


In February and in January the distribution of the height field 
agrees most with the correlation field shown by Horel and Wallace. 
Especially in February, a ring of relatively low heights surrounds the 
North Pole, with the ring displaced southward along 90W. This agrees 
also with theoretical findings cited. In January the pattern is almost 
as coherent, but in December and even more in November, the pattern 
deviates more from the "ring pattern" evident in February. 

Fourier harmonic analysis shows that wave numbers 1, 2, and 3 con- 
tribute substantially to the pattern, especially in middle and high 
latitudes. 

Finally, an interesting 10 year repetition occurs in the 700 mb 
heights at 60N, 90E. The cause is unknown; but the possibility is 
raised that SST in the Pacific and Atlantic Oceans may respond to gyres 
with different periods that are harmonics of 10 years. 
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Northern Hemi sphere ’OO ah Heights and 
Pacific Ocean Temp ratures In Winter 

by Sigmund Fritz 


Introduc tion 

Horel and Wallace (1980) (abbreviated H-W) show that a correlation 
pattern exists between the 700 mb heights over the northern hemisphere 
and certain patterns of the sea-surface temperatures (SST) of the Pacific • 
Ocean. For the SS7, they use an empirical orthogonal function (EOF) of 
Pacific SST given by Weare, Navato and Newell (1976). H-W show a time 
series graph of normalized t'ST from 1949-1979. They indicate by shading 
their "subjective definition of major warm SST episodes in the equatorial 
east Pacific." (It should be noted, however, that the EOF SST patterns 
characterize the SST over the whole Pacific Ocean. On the average, when 
the equatorial east Pacific Ocean is warm in winter, the SST is warm 
also in the central Equatorial Pacific; but in mid-latitudes at about 
40N, the central Pacific is cold then. It is really this total pattern 
which characterizes the "warm" SST years selected by H-W.) 

Horel and Wallace analyzed their data for winter seasons; namely for 
December, January and February grouped together. It is interesting to 
see how consistent their results are for all the warmest SST years which 
they selected and also for the coldest years which can be selected from 
their diagram (their Figure 2). They found that the largest correlation 
coefficient between SST and 700 mb heights for the winter season occurred 
near 60N, 160E. The correlation coefficient there was -0.66. To test 
for consistency, the departure from normal (DN) of the 700 mb heights at 
that location were tabulated in Table 1. There were 10 SST-warm years and 
8 SST-cold years (i.e., years when the Equatorial Pacific Ocean SST was 
anomalously warm or cold.) 

We note in Table 1 that in all 8 Febru ary m onths with SST cold, 
the 700 mb height was above normal at 60N, 160E; also in 8 out of the 10 
SST warm years, the 700 mb height was below normal. In January , almost 
equal consistency was found. In 7 out of 8 SST cold years the 700 mb 
heights were above normal. In 7 out of 10 SST warm years the 7.00 mb 
heights were below normal. P< rsistence of the DN sign, from January to 
February, was more consistent in the SST cold years than in the SST warm 
years. 


However in December, consistency was not evident. Half the 700 mb 
heights were above normal and half were below normal in both the SST 
warm years and in the SST cold years. 

Therefore, the results shown in Table 1 are consistent with H-W 
as to the sign of the 700 mb DN, but mainly in January and February. 
Presumbably H-W might have found an even larger correlation coefficient 
at 60N, 160E, if they had confined their study to January and February, 
although the magnitude of the DN's are also a factor. 


| 
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Table 1 


Depa rture f roa Norm al 
Me an Mon thl£ 700 wb Helghta (a) 
60 N , 160E 


S S T w a r m_ yea r a 


Years 

Dec 

Jan 

Feb 

1951-52 

-15 

-28 

- 21 

57-58 

-16 

-31 

+ 20 

58-59 

+ 1 

+13 

- 55 

63-64 

-55 

-57 

+ 11 

65-66 

-19 

-61 

-114 

68-69 

+ 5 

-32 

- 37 

69-70 

-60 

+80 

- 34 

72-73 

+15 

-49 

- 47 

76-77 

+40 

+14 

- 59 

77-78 

+48 

-34 

- 20 


numbers of negative cases 
(out of 10 possible) 


SST cold year*. 


Years 

Dec 

Jan 

Feb 

1949-50 

-60 

+ 61 

+19 

54-55 

+ 59 

- 20 

+29 

55-56 

+45 

+ 70 

+68 

64-65 

+29 

+ 21 

+ 9 

67-68 

- 8 

+ 43 

4-86 

70-71 

-28 

+ 91 

+81 

73-74 

-12 

+143 

4-39 

75-76 

+35 

+ 2 

4-25 


number of positive cases 


(out of 8 possible) 
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Years 

1951-52 

57- 58 

58- 59 
63-64 
65-66 

68- 69 

69- 70 
72-73 

76- 77 

77- 78 


Table 2 

Departu re f rog No real 
Me an Mont hly 700 ro b H eight! (a) 
MNjlJIOW 


SST wa™ /ears 


Dec 

Jan 

Feb 

-51 

-46 

+10 

-36 

+59 

+32 

-*•12 

-11 

-8 

+48 

-28 

-9 

-15 

+16 

-20 

+32 

+4 

+50 

+58 

+ 7 

+23 

-29 

+17 

+ 39 

-3 

+33 

+50 

-17 

+89 

+101 


number of positive cases 
(out of 10 possible) 



SST cold /ear* 


Years 



Dec 

Jan 

Feb 

1949-50 

-8 

-68 

+36 

54-55 

-11 

+10 

-8 

55-56 

-1 

+ 74 

-27 

64-65 

•37 

-4 

-68 

67-68 

-5 

+11 

+78 

70-71 

-22 

-17 

-0 

73-74 

+32 

-54 

-16 

75-76 

-14 

+16 

-41 


number of negative cases 
(out of 8 possible) 



< 


lie 


Horsi and Wallses found a correlation coei'ficlent of 40.45 at 60H, 110U* 
Tabla 2 shows tha 700 mb DN aa In Tabla 1. Hart tba results ara not aa con- 
sistent aa they ware In Tabla 1. In February, 6 out of 8 700 mb DM' a ware 
negative for the SST cold years, only 4 out of 8 negative DN's occurad In 
January. In the SST warm yeara tha number of poaitive DH'a were 7 out of 
10 in both February and January. In December, tha SST cold yeara ahowad 
consistency with 7 out of 8 negative DH'a. But for the warm SST yeara, 
December had only 4 out of 10 positive DN'a; i.e., more than half with tha 
"wrong" sign. 

• 

So over all, January and February seem to show greatest consistency 
and agreement with H-W as to sign of the DN, but December seems mora 
doubtful in that respect. 

. # 

Because of the month to month variation of the results at the two 
locations in Tables 1 and 2, we decided to examine the 700 mb height 
variation over most of the Northern Hemisphere Repara rely, for each 
calendar month. A magnetic tape had already been furnished by the N0AA 
Long Range Predict ion CV.'up. This tape Included mean monthly 700 mb 
height data from January 1947 through May 1979, and had been used exten- 
sively to study the 700 mb variation, especially at 60N. (Fritz, 1980). 

The tapes contained data from 15N to 90N. 

We accepted the SST character ! '-at ion as given by Horel and Wallace for 
the years listed in Tables 1 and ?. For the 10 SST warm years, the average 
700 mb height was computed at each 10 latitude and 10° longitude grid 
point from 20N to 80N. The same computation was made for the SST cold years. 
Then the following difference, AH, was computed. 

AH - H - H 
w c 

where H is the 10 year average for the SST-warra years and H is the 8 
year average for the SST cold years. This was done separately for each 
calendar month. And because Table 1 suggests that the February results 
are most consistent we shall describe the maps of AH In the order February 
January, December. We shall also discuss the results for November. 

February AH t 

The map of AH for February is shown in Fig. 1. The largest value of 
AH was -94 m at 60N, 170E. In addition there was a belt of negative heights 
surrounding the pole in a ring which was displaced' southward approximately 
along longitude 90W. Along this circumpolar belt were several negative 
height centers. The most extensive one was' the area centered near 55N, 175E; 
another center, with magnitude about -35 m was located just off the U.S. 
east coast, centered at almost 35N, 70W. A third center, with magnitude 
about -35 m, lay near 60N, 10E; and finally a center of about -57 ro appeared 
near 80N, 90E. Angell (1980) also found that warm sea-surface temperature 
in the tropical Pacific was related to "displacement of the polar vortex . 
toward the region 180-90° west". 


An area of positive heights appeared over Northern Canada and the adjacent 
Arctic Ocean; its central value was about 4-36 ■ and it was located near 65N, 
125W. This is not far f roa the center of positive correlation found by H-W 
for the winter months. 

The striking feature of the circumpolar ring of negative AH, was its 
similarity with theory and with observation. H-W found a ring of centers 
of negative correlation, (as described above) lying also in a quasi-circum- 
polar ring with four centers of negative correlation coefficients. Two of 
these centers lay near the centers In Fig. 1; namely, 
the one already mentioned near 60N, 170E, and the one off the U.S. east 
coast. They show also another center in the Pacific (r - -0.59) near AON; 
1A0W, and one in the Atlantic (r » -0.30) near 50N, 20W. 

This agreement is, of course, to be expected because both studies use 
the same 7 00 mb data as the ones used in Fig. 1. However, their results 
are for the whole season, December, January, February, and they presumably 
used all years, not only the "warm" ana "cold" years. So some of the 
disagreement between the height distribution, AH, in Fig. I and the distri- 
bution of H-W's correlation coefficients, is doubtless due to these differences 
is the data base used. 

Horel and Wallace also find a region of high correlation, fO.63 near 
25N, 155E. In Fig. 1, a region of positive values of AH appears thers 
also. Other areas of positive AH are evident over the Mediterranean and 
over Asia. 

Turning to theory, Opsteegh and van den DjoI (1980) have computed the 
geopotential height response to an SST anomaly in the tropical ocean. They 
also show comparison with Fcowntree's GCM for 300 mbs. Both theoretical 
results contain the ring of low heights displaced from the pole along a 
meridian which Is somewhat downwind from the heating area. In Fowntrce'i 
resjlts, which refer to a tropical heating in the Atlantic, the largest 
negative height anomaly appears near latitude 50N, somewhat upwind (west) 
of the longitude of tropical heating. This too is in qualitative agreement 
with Fig. 1. 

Thus, the February results agree with H-W's seasonal results, and 
with theoretical results too. In Jnnu try, many of the same features 
appear; but differences are present also. The map for AH in January is 
shown in Fig. 2. 

•• 

January AH. 

In January, the largest negative value of AH was -83 ra, also at AON, 

170E. This time, the band of negative AH, does not completely surrdund 
the North Pole; it Is Interrupted by a positive area which stretches from 
Greenland to Arabia. However, the band of negative AH, as in February, 
does reach from Kamchatka across the southern U.S. and Into the Atlantic. 

Again a center of negative AH, value about -35 o, lies along the southeastern 
United States; this time, however, the center seems to be somewhat further 
west than the corresponding low center in February. The low in the Atlantic 
is considerably further south than in February; in January, it was located 
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near 40H, 25W with a valua of about -25 •. Tha negatlva araa which waa 
present in February at BON, 90B, la not present in January. Instead tha 
region of tha whole Arctic basin shows positive values of AH, with a 
maxi nun appearing over Greenland (value about +45 ■) instead of over Canada 
where it appeared in February. The positive area in the Western Pacific 
is again present as it waa in February. A high *, +4 5m, appears near 
55N, 40E, somewhat to the northwest of the weaker high which appeared in 
February near the Aral Sea. 

In summary, AH in January appears quite similar to AH in February 
over the Pacific Ocean and over the southern United States, Elsewhere, 
the values of positive AH seem to have been displaced eastward from * 

Western Canada to Greenland and perhaps have expanded Into the Arctic 
Ocean bordering Siberia'. And as we shall see (Figs. 10 and 11) the 
weak cental in the Atlantic and its southward displacement, give more 
prominence tr> wave no. 1, while wave no. 2 is suppressed in amplitude. 

December AH^ 

In Iocember, further changes f r ora February were observed. The map 
of AH for December is shown in Fig. 3. In December a center of negative 
AH appears over the Pacific near 45N, 1.60W, Its values was -48 m. This is 
considerably weaker than the values in January or February, and the center 
was farther to the southeast. The band of negative AH, again stretches 
across the Pacific Ocean, then weakly across Mexico. This time the negative 
area "trough” passes well inland over the eastern United States and eastern 
Canada, then to a center, (AH * -95 ra) , near 50N, 20W. It then continues 
across Europe to Turkey. Was the large negative AH in the Atlantic related 
to the SST's in the tropical Atlantic Ocean? It would be interesting to 
investigate that possibility. 

The positive values of AH now again appear over western Canada and 
the adjacent Arctic Ocean (AH » +45 m) . The high AH which had appeared in 
February near the Aral Sea and in January south of Finland, now appeared 
still further northwest ' A H -* +60 ra) centered over Scandinavia. 

To some extent then the negative pattern In tha Pacific appears again, 
although weaker and further east. The positive area in the southwest 
Pacific seems to have shifted towards the northwest, with a high valua of 
+20 m appearing over the south China Sea. 

Of significance for the eastern United States is the location of the 
lowest AH there. In January and February, centers of negative AH lay 
along the east coast. Therefore, from geostrophlc wind consideration* 
one might expect below normal temperatures in SST-warm years along much 
of the U.S. east coast in those months; the reverse would be expected in 
SST-cold years. A preliminary review of the "temperature departure from 
normal" maps in the Monthly Weather Review Indicates that in February 
this expectation was often fulfilled. 
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November A H. 

__ From Pacific tropical SST kindly supplieJ by Dr. E. Rasmusson (Climatic 
Analysis Center, NOAA) , it. "appears that SST's in November in general, had 
the same sign' as the SST anomalies for December, January and February. It 
is therefore interesting to examine the November values of AH. 

Fig. 4, shows the map of AH for November. Again we find a low center 
in the central Pacific with a value of -59 m. This center was located near 
50N, 170W. There also was a weak negative area over the southern United 
States, with a center of a’lout -20 m near southern Arizona. Although 
negative values still appear over the southeastern U.S., the values are 
very small. Another negative center appears over the Atlantic Ocean, 

(AH " -46 m) , near 50N, 20W. The negative zone, although weak, continues 
across central Asia. 

Several centers of positive AH appear around the hemisphere near latitude 
and 40N; a zone of positive AH lies over Greenland and the Arctic Ocean 
north of Siberia. One of the positive AH centers appears near the west •«** 
coast of North America. Interestingly, the axis of positive AH, stretching 
from Greenland across the Arctic Ocean to Northern Siberia, is essentia lly 
perpendicular to the axis of positive AH across the Arctic in December. 

In November (Fig. 4) again the negative band of AH is not continous 
around the Earth. If the SST was producing the negative area over of the 
North Pacific, in November and December the effect does not seem vigorous 
enough to produce a continuous ring of negative AH as it did in February. 
Still there Is a suggestion of a continuous ring or spiral of negative AH 
that stretches from North Canada, across the Pacific into Asia, thence 
across Scandinavia and England into the Atlantic; after a weak break 
(AH = +5), that negative band continues across the southern United States, 
to Northern Mexico and into the tropical Pacific. 

Summary of Pacif ic^ and Atl ant ic Ocea n nega tive centers 

In all four months, the negativ ■ AH centers in the Pacific ar *d Atlantic "v 
are summarized in Table 3. 


Tabl e 3 

Locat io ns of negative AH c enters 


Pacific Ocean . Atlantic Ocean 


Nov 

Lat. (H) 
50 

Long. 

170W 

# 

AH (in) 
-59 

Lat. (N) 
50 

Long. (W) 
20W 

AH (m) 
-46 

Dec 

47 

170W 

-48 

* 50 

20W 

-95 

Jan 

55 

180 

-83 

40N 

20W 

-25 

Feb 

55 

175E 

-94 

60N 

10W 

-35 


* • • I 
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The patterns for all four months show persistent negative AH centers in 
the central Pacific near longitude 180*. In the Atlantic Ocean, centers 
of negative AH also appear near longitude 20*W in all four months. This 
suggests that wave number 2 plays a prominent role in the response of the 
atmosphere to the Influence of the Pacific SST anomalies. However, because 
the magnitude of the AH centers in the two oceans varies, wave number 1 and 
other Fourier components also affect the results. To examine this further 
the Fourier components for February and for December were computed. 

Fourier Wave C o mpone nts. 

T he "normal he ight f i el d" , 

Fig. 5 shows the phase and amplitude of the trough of the "normal" 700 
mb heights for wave no. 1. The "normals" are the average heights for the 
period 1951-1973. The phase varies very little from October through March, 
from 20N to 60N; Ellasen (1958) had already found similar results. At 
latitude 70N and 80N more variation occurs, but the amplitudes are relatively 
small there. Fig. 5 is also in agreement with the values of van Loon, Jenne 
and Labitzke (1973). 

The amplitudes are largest near latitude 50N, varying from about 70 m 
in December to 59 m in February. In October and March the amplitudes seem 
substantially smaller; namely, about 48 m. Therefore unless interaction of 
the SST with the mean circulation is very sensitive to small changes in the 
"normal" height field, we probably cannot depend on interactions between the 
SST effect and the "normal" height field to "explain" the observed difference 
between the December AB-.(Fig. 3) and the February AH (Fig. 1). Opsteegh and 
van den Dool (1980) found that the mean zonal wind variations between winter, 
fall and spring, do not produce weaker effects in the response of the SST. 

If anything they found even more pronounced response in fall and spring than 
in winter. 

Fig. 6 shows the "normal" phase and amplitude for wave no. 2. Again 
the phase wave varies little from month to month except at latitude 30N, 
where the amplitude is small. It is interesting to note that the slope 
of the phase line in Fig. 6 is from east to west with increasing latitude 
this is opposite to the slope for wave no. 1 (Fig. 5), a result also found 
by Eliasen (1958) and emphasized by Austin (1980). 

The amplitude in wave no. 2 is again a maximum at 50N. Moreover, the 
amplitude is similar in December, January and February, i.e., between 50 m 
and 59 m. Now however, the amplitude in November is substantially smaller 
and even smaller in October; in March the amplitude is also smaller than in 
the winter months. So again, it might be difficult to look for differences 
between December and February (Fig. 1 and 3) based on variations of the 
"normal" height fields in those two months. 

Finally, Fig. 7 shows the phase and amplitude for wave no. 3, of the 
normal 700 mb height field. Again there was almost no variation of the 
phase from month to month for October through March, this time at all 
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latitudes. The amplitude was somewhat smaller in December than in February 
(46 m verses 59 m). It is interesting to note that at latitude 50N, the 
amplitude for wave nos. 2 and 3 was about the same in March as in December; 
and the amplitudes were also mainly similar for those months at all latitudes. 

Compon e nts of the AH field. 

What about the variation of the wave number component in AH; I.e., 
in the heights for the SST warm years minus the SST cold? Fig. 8 shows 
the phase of the minima for AH for wave numbers 1 through 4, for February. 

Fig. 8 shows that the slope of the phase with latitude was from east to 
west with Increasing latitude. This is the same sense as the slope of wave 
no. 2 in the "normal" height field; only now it is observed even in wave 
no. 1. At latitude 60N and SON, especially, all the wave number components 
contributed to the negative values of AH. In fact all the phase lines 
appear to intersect at about 45N. Depending on the amplitudes of the Fourier 
components, we might find a large negative AH at 45N. However, the amplitude . 
(Fig. 11) of wave no. 2 (40 m) at 60N. was a maximum, decreasing in magnitude 
both north and south of 60N. 

Figures 10 through 14 show the amplitudes of the Fourier components of 
AH for November through March. In December, January and February, the ' 
amplitudes of wave numbers 1, 2 and 3, are usually the largest at latitudes 
50 to 80. In January, as noted earlier, wave number 1 dominated. Still, 
depending on the phase relationships shown in Figs. 8 and 9, contributions 
from wave numbers 4 and_ 5 may also be important. 

The sum of the amplitudes at 60N is not sufficent to yield the observed 
large negative values of AH, namely -94 m in Fig, 1. To attain such large 
values we need to include the zonal average value of AH; i.e., wave no. 
zero. The zonal mean values are shown in Table 4. 

Table 4 - Zonal average value of AH (m) 

Latitude Decem ber January February 


80N 

+ 4 

+25 

- 6 

70 

+ 2 

+ 5 

-10 

60 

-10 

- 2 

-10 

50 

-18 

-11 

-11 

40 

- 6 

-11 

- 6 

30 

+ 6 

0 

+ 4 

20 

+ 7 

+ 7 

+10 
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Table 4 shows that for latitudes 40 to60N AH was negative In all three 
winter months. At latitude 20 and 30N, there was a tendency for positive 
AH to appear. This effect would increase the west winds near 40N in the 
SST warm years relative to the SST cold years. A tendency for the subtropical 
jet stream to strengthen significantly in response to sea surface temper- 
ture changes in the equatorial Pacific was also noted from numerical 
experiments at GFDL (1980), 

The slope of the phase lines in Fig. 8 suggest that one of the effects 
of the SST anomalies might be to strengthen the influence of wave no. 2 in 
the "normal" 700 mb height field. In the normal field, wave no. 2 is 
prominent because of the effect of mountains and continent-ocean contrast on 
the. -atmospheric circulation (Manabe and Terpstra, 1974). 

Discussion 

The "normal" 700 mb height field does not seem to be sufficiently different 
in December from February to explain the difference between the AH fields 
in those months (Fig. 1 and 3). What other effects might have done that? 

Of course the possibilities are very large. Among the possibilities, are 
the actual detailed variations in SST not only in the tropical Pacific but 
also in mid-latitude SSI variations. Numerous empirical studies (e.g. Namias, 
1979) and many theoretical studies discuss relationships between central 
Pacific SST and aspects of the atmospheric circulation. For example Fritz 
(1980) has found that the correlation between the SST in the central Pacific 
and the mean monthly 700 mb height at 60N, 170W was positive in all winter 
months and in two sepaf&te decades. Therefore it might be useful to examine 
the SST in the central Pacific for the "warm years" and the "cold years" in 
Table 1. We note again that, on average, the SST in central Pacific would 
tend to have opposite signs from the SST in the equatorial Pacific (Weare, 
et al, 1976). 

Another factor might have been the SST in the Atlantic, Fowntree (1976) 
had shown that the SST in the Atlantic would affect the circulation of the 
atmosphere over the Atlantic. Perhaps a large SST anomaly in the Atlantic 
produced the observed large negative value of AH (-90 m) in the Atlantic 
in December (Fig. 3). Depending on whether anomaly in the Atlantic were t . 
in phase or out of phase with the SST in the Pacific, the Fourier wave 
number components of the AH pattern might be affected in different ways. 

Other factors might be the snow and ice distributions. These would not 
only vary from November to February but also from one year to the next. 

The southward extent of snow and ice, and the thicknesses would of course 
influence the surface temperatures. Can differences in the surface 
temperatures, due to ice and snow, also influence the general circulation 
of the atmosphere, in addition to sea surface temperature effects? 

In December (Fig. 9), the slope of the phase line of wave no. 1 is 
different from the slope in February (Fig. 8). Although, there is some 
ambiguity about the slope direction between 60N and 50N, from 50N to 30N,* 
the slope is opposite in December from the slope in February, In December 
for w f / e no. 1, the slope is more like the slope of the "normal ,, wave 
no. 1 phase. The shift of the wave no. 1 minimum in December from the 
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Pacific north of 60N to rhe Atlantic sector at 40N and SON* is connected 
with the large negative value of AH (-95) in the Atlantic in December. 
Whether this shift was caused by a large SST anomaly in the Atlantic 
Ocean is unknown at present. 

The possibility that the SST in the Pacific and Atlantic might 
affect the large scale circulation simultaneously raises many possibilities. 
For example, the Pacific and Atlantic Ocean gyres have different periods 
(Reiter, 1979). If the gyres carry with them SST anomalies, then observed 
cyclic variations in the atmospheric circulation might be related to 
combinations of two or more periodic variations of SST in the two oceans. 

Earlier Frits (1980) showed that the lowest 700 mb height anomaly at 
60N, 90E tended to occur at 10 year intervals. Specifically, in the sum 
of wave no. 1 and wave no. 2, the lowest 700 mb height at 60N did occur 
at about 90E in 1948, 1957, 1958, 1967, 1968, 1977, 1978. To see what 
happened in all the years from 1948 to 1978, the 700 mb height departure 
from normal (DN) , is shown in Fig. 15. The similarity in the DN's for the 
decade 1948-1957 and 1968-1977, is striking. Moreover even in the decade 
1958-1967, the years 1958, 1959, 1966, 1967 follows the trends in the other 
two decades. Fourier harmonic analysis of the data in Fig. 15 show that the 
period 10/4 =» 2.5 years had the largest amplitudes in the decades 1948-57 
and 1968-1977. The time resolution of the data by Fourier analysis is not 
very detailed. It is nevertheless interesting lo note that 2.5 years is 
not far from the quasi-biennial period; also 2.5 years is not very far 
from 28 months, twice .the Chandler wobble period . of about 14 months. * 
Maksimov (1958) claimed that the surface pressure near 90E in high latitudes 
exhibited a periodicity which he attributed to the Chandler "polar tide" 
and to the non-spher ical shape of the Earth. Whether periodic SST 
variations in the Earth's oceans, or the Chandler wobble, or some other 
effect is the "cause" of the regularity in Fig. 15 is at present uncertain. 

In summary, the circulation change differences from November to 
February (Figs. 1-4) in relation to SST seem consistent in same respects 
but different in other respects. The cause for the difference are unknown; 
perhaps surface effects or atmospheric «ffects such as vertical stability 
might be causative factors. It \<rill also be interesting to extend the . 
analysis to other months, pe< baps for the whole year July to June, extending 
Figures 1-4 in both time directions. It should also be interesting to 
examine maps such as Fig. 1-4 for sea-level pressure data. 
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Figure 1 : Mep for ft£ruarg of AH for 700 nb height.. AH - H - fi . 

where Hi. the average 700 mb height for the 10 year” wheS 
he tropical Pacific SST we. anomalously warm, and ii wa. the 

**4 When 1 th * SST WM cold - Note the circumpolar 
ring of negative value., displaced southward along about 90°W. 

Figure 2 : Map for Ja nua ry of AH. (see Fig. 1). 

Figure 3 : Map for December of AH. (see Fig. 1). 

Figure 4 : Map for November of AH. (see Fig. 1). 

Figure 5 : Phase and amplitude (ra) of the "normal" 

700 rob height field. The phase is the longitude of the 
minimum height for that Fourier component. 

Figure 6 : Wave number 2. (otherwise the same as Figure 5). The bar 

shows the range of phases at latitude 30N. December phase was 
at the western edge, January at the eastern edge. 

Figure 7 : Wave_numjber _3. (otherwise the same as Figure 5). 

Figure 8 : Phase (Longitude of minima) of AH, for February; wave 
numbers r l through 4. 3 

Figure 9 S ^hasej (Longitude of minima) of AH for December; wave numbers 

Figure 10s Wavermmber 1. Amplitude (m) for November through March 
from latitude 20N to 80N, of AH. ‘ 

Figure 11: Same as Fig. 10, exc ept wave num ber 2. 

Figure 12: Same as Fig. 10, except wajve number 3. 

Figure 13: Same as Fig. 10, except wave number 4 . 

Figure 14: Same as Fig. 10, except wave number 5 . 

Figure 15: DNWIW2 is the sum of Fourier components for wave number 1 
plus wave 2 of 700 mb height departure from normal (DN) for * 
ecember, at 60N, 90E. The abscissa show the years arranged 

and ’ 1948 - 1 ? 57 * 1958-1967, 1968-1977. TT,e years 1958 

and 1968 are repeated in order to include 1978. 

Figure *16: "Normal" 700 mb height map for February after latltu<H„„l 

heights sho™ t L‘TeetangiIarho 1 ,L r0 ”l“ ean ,'’ Cf8htS - Z ° nal avers S' 

22 'Ai™ H < ’ 1SCUSSed *«. hu t n fnc luded“f or^compar is2 ' 
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1. Introduction 

The research of this project la directed towards seasonal climate 
variations; that la, fluctuations In the meteorological variables on tlaa 
scales of about a month and on spatial scales of 1000* s of kilonetara. 

It is an underlying assumption that the "climate" on these scales can ba 
resolved if the causes of planetary-scale, quasl-statlonary disturbances are 
understood. Then, since there Is a correlation between the amplitude and 
phase of the planetary waves and the location and Intensity of synoptic- 
scale disturbances (e.g., Blackmon, 1976, Lau, 1979), it may be possible to 
predict the fluctuation in regional climate through the use of empirical 
statistics. Although the prediction, per se, Is not the object of our 
project, it Is the underlying motivation for the study of planetary-wave 
dynamics. 

The causes of the planetary wave-patterns have been recognized for many 
years since the rather successful studies of Charney and Eliasuen (1969), 
Snagor lnsky (1953) and Doos (1963). Using linear steady-state models, they 
were able to explain the existence of the stationary, planetary-scale waves 
In middle latitudes. From these early studies the Importance of topographic 
forcing and thermal heating were recognized. 

It Is somewhat of a surprise, then, that the prediction of the long waves 
in forecast models has been so difficult. Apparently, the nonlinear dynamics 
of the long waves on a sphere are a more complex process than was visualized 
with early models for stationary flow. Recently, there have been several 
developments which have illustrated the Important dynamics of planetary waves. 

* Spherical geometry is important in the description of ultra-long 
waves (Kallen, 1980; Sommervllle, 1980). 
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* Thera arc Important "talacannactiona" batwaan the forcing in the tropica 
latitudes and the raaponaa In middle latitudes (Moral and Wallace, 1980; 
Opateegh and van dan Dool, 1980). 

* Multiple equilibria In a nonlinear dynamical system may exist for 
soma conditions of mechanical and thermal forcing (Cherney and DaVore, 
1979; Vlekroy end Dutton, 1979; Lorenz, 1980). 

These developments give some Insight Into the possible prediction of 
planetary-scale waves and their seasonal variations. That le. It may be 
possible to simulate the ultra-long waves beyond the time period of deter- 
ministic forecasting if we can correctly model the dynamic response of the 
atmosphere to slowly changing, but persistent thermal and topographic forcing. 
The minimum requirements of such models are: 

* Global in lateral extent, two layers. 

* Non-linear dynamics. 

* Resolution of planetary scales. 

* Balanced conditions for pressure/velocity determination. 

* Capable of long-term integrations; l.e., energetically consistent 
and able to explain long-term changes in the distribution of mass 
and kinetic energy. 

A low-order, spectral model with these characteristics has been developed 
and is operating on the GLAS computer (Section 3). A steady-state, geostrophlc 
version of this model has recently been constructed also by Asha (1979). This 
type of ondel recof 'zee the collective forcing of all scales less than 
planetary size only implicitly through an equivalent imposed heating function. 

Of course, it is a serious limitation to ignore the explicit interaction 
of cyclone-scale eddy fluxes with the planetary scale steading waves. However, 
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if the forcing by synoptic, barocllnlcelly unstable disturbances is really 
the most Important forcing function of planetary- scale dynamics, efforts to 
make dynamic, seasonal climate forecasts are seriously Impaired, since the 
forecast will depend on statistical products from an Instability process which 
cannot itself be forecast more than a week or so in advance. In fact, we 
have some evidence that this pessimistic view is not warranted; reasonable 
planetary wave structures have been determined by many authors from topo- 
graphic and external heating sources without resort to explicit forcing by 
the transient eddy fluxes. However, this question is by no means settled 
(Gall, et al. , 1979). To predict the variations in seasonal climates, it 
may be necessary to have some explicit knowledge of the transient eddy fluxes. 
In that case, two choices are apparent; the eddy fluxes must be parameterized 
or specified from climatic data or contemporary analyses of synoptic-scale 
observations. 

Another important question for a model of planetary waves concerns the 
mechanism of "blocking" which is usually associated with stationary waves of 
synoptic scale. To explain this phenomenon there are several theories which 
have been proposed which use physical mechanisms that are also important in 
the dynamics of planetary waves: 

NON-LINEAR THEORIES 

* Multiple equilibria (Chamey and DeVore, 1978; Hart, 1979; Kallen, 

1980; Lorenz, 1980; Vickroy and Dutton, 1979). The possibility of 
the existence of multiple equilibria in a nonlinear dynamic system 
suggests that the "blocking" phenomena may simply represent dynamic 
solutions in a different attractor set of a stable equilibrium point 
which is distinct from the usual case of quasi-periodic flow. 
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* Wavt; wav Infraction resonance (Egger, 1978). Blocking action 
occurs In preferred regions of tha standing wave forced by planetary 
scale heat ? 'i and topography. Travelling waves near resonance (C s 0) 
Interact nonlinearly to create a smaller-scale, persistent distur- 
bance resembling a blocked flow condition. 


LINEAR THEORIES 

* Local resonance of Rossby lee waves caused by topography (Kalnay-Rivas, 
1980; Faller, 1980). The response of the atmosphere to topographic 
forcing depends on the local flow condition at the mountain ridge and 
the existence of a disturbance upstream to deflect the flow. 

* Planetary-wave resonance controlled by mean flow conditions, thermal 
and mechanical forcing (Tung and Lindzen, 1980). Medium-scale, propa- 
gating waves are excited by fluctuations in thermal forcing or changing 
flow conditions over variable topography. These waves are easily 
trapped (in the vertical) by the zonal wind field, and under the proper 
mean flow conditions, may be near the resonance for stationary waves, 
thereby producing a "blocked" condition. Blocking by longer, planetary 
waves is more unusual, but its influence is felt for a longer period 
and over a greater area, even penetrating into the stratosphere. 

* Unusual thermal forcing drive both exceptionally large amplitude 
stationary and propagating planetary waves to produce an (apparent) 
large-scale blocked flow (Austin, 1980). Subsequent steering of cyclonic 
storms by the pit nitary mean flow reinforce, intensifies and maintains 
the blocked flow. Conditions may permit the development of inter- 
mediate wave numbers (e.g., wave number 4) which are easily resonant 

to the usual mean atmosphere flow conditions; this is the usual scale 
of blocking anticyclones. 

In general, theories of blocking depend on the strength of th«! planetary 
thermal forcing or the state of the atmosphere (resonance) or the nature of 
the quasi-equilibrium and its dynamic fluctuations. But these same mechanisms 
are responsible to some degree for the existence of planetary waves and their 
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scale flow. 
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truncated apactr.l model already developed for thl, p„rp„... 


2. Review 


Recent diagnostic studies! e.g. Van Loon (1973), Leu (1979) end 
Boettger (1979) Indicate that the seasonal climate, as defined In the 
introduction, is determined largely by the long-wave components (wave 
number leas than or equal to A). Van Loon (1973) finds, for example, 
that 962 of the variance of the winter normal 300 mb geopotentlal 
height is contained in the first four wave numbers at midlatitudes. 
Blackmon (1976) performed an analysis of the same field in terms of 
Fourier series in time and Spherical Harmonics in space, and found 
most of the variance concentrated in the long-wave, slow moving (quasi- 
stationary or stationary) part of the spectrum. Lau (1979) finds that 
the contribution of the stationary components of the local vorticity and 
heat transport budgets' exceeds that of the transient component by a fac- 
tor of 2 to 5. Similar conclusions are reached by Wallace (1978) in his 
discussion of the role of stationary versus transient components in 
maintaining the general circulation. 

The transient components are present at all scales. The short scale 
features (baroclinlc waves) contain most of the variance associated with 
the transient components. Statistically, these waves are correlated with 
the stationary ones. Boettger (1979), among others, clearly shows in- 
creased cyclogeneais east of stationary long wave troughs. t 

The transient long waves have been analyzed by Madden (1979). They 
are propagating Rossby waves. They are of some importance to the 
stratosphere (Charney and Drazin, 1961; Matsuno, 1970), and some slowly 
moving modes near resonance may contribute to important tropospheric 
anomalies (Lindzen, 1980). 
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The stationary long waves arc a reaponae to tha large-scale atmospheric 
forcing flelda: 

(a) orography, 

(b) diabatic heating by turbulent heat tranafer from the ground, 
release of latent heat, direct radiatlonal heating and cooling, 
and radiative transfer within the cloudy atmosphere. 

(c) internal mechanical dissipation, and 

(d) energy transfer from an ensemble of cyclone-scale disturbances. 

The response of the stationary component to forcing was first investi- 
gated using linear models on a beta-plane. Charney and Eliassen (1949) 
considered the effect of orography using a linearized equivalent baro- 
tropic model. The 500 mb. geopotential field was decomposed on a latitude 
circle into a large-scale stationary component and a shorter, transient 
component. The predicted height field was accurate enough to establish 
the importance of orographic forcing. 

Subsequently, Smagorinsky (1953) treated the effect of sinusoidally 
distributed heat sources on a simple baroclinic model with constant lapse 
rate and wind shear, including friction. He found the response of the 
perturbation height field, with a trough east of the point of maximum 
heating. These results agree qualitatively with observation. 

Since the basic state used by Smagorinsky was not entirely realistic, 
other investigators have refined his results. Doos (1968) used a lapse 
rate and wind shear varying with height. His results are not markedly 
different from those of Smagorinsky. Saltzman (1965) considered heating 
and orography simultaneously, using a linearized model. Recently, Tung 
and Lindzen (1979) performed a similar analysis, but using a multiple 
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perturbation technique. They have proposed an explanation for rare but 
large anomalies in the planetary waves by a "blocking" mechanism due to 
a resonance of a nearly stationary Rossby wave. 

All these analyses rely on linearized dynamics. However, they are 
limited by (a) the necessity of prescribing an equllbrlum basic state, 
which is usually assumed to be consistent with some imposed forcing, and 
(b) by the resonance effects which may be unrealistic. These shortcomings 
are reviewed by Ashe (1979). Furthermore, the studies are carried out in 
Cartesian geometry which is not appropriate for planetary waves and which 
may provide spurious reflections of energy by the lateral boundaries. 

General circulation models have also been used to analyze the effects 
of the forcing fields. Manabe and Terpatra (1974) computed the Influence 
of orography in the GEDL model and found significant differences in the 
distribution of energy between stationary and transient components. The 
response of GCMs to thermal forcing at locations in middle latitudes has 
not been very strong (Kutzbach et al. 1977; Houghton et al. 1974). However, 
there is some evidence that the thermal forcing from the subtropics strongly 
influences the atmospheric circulation in the GCM (Chervin, 1980). 

TheGMCsare, as a rule, unable to reproduce planetary scale, persistent 
features (Sommerville, 1980), or persistent "blocked" flow, although 
Miyakoda (1980) has been able to show improvement in this respect. There- 
fore, it is apparent that studies of simplified models are needed to 
clarify the important dynamic processes in complex GCM experiments. 

Simplified models attempt to lower the order of the system by 
expansion in series of orthogonal functions or by asymptotic expansion. 

In some cases analytic solutions are possible (Baer, 1970, 1971); however. 


eh* truncation i* necessarily sever*. An alternative approach is to 
us* the methods developed for analyzing systems of nonlinear ordinary 
differential equations. These methods are treated in Mlnorsky 
(1962). The methods emphasize the study of the system in phase space in 
the vicinity of the points where time derivatives vanish (critical points, 
or steady states). In the case of the atmosphere or simple fluid flows 
each steady state is Identified with a particular flow regime. 

Lorenz (1962a, 1962b, 1963a, 1963b) carried out studies motivated by 
laboratory experiments with rotating flows, where the fluid exhibits 
several flow regimes determined by the imposed forcing (in this case a 
thermal gradient and the rotation rate). He found different, steady 
state solutions as the forcing was changed, i.e., the Hadley and Rossby 
regimes. The Hadley regime contained a number of subclasses identified 
by wave numbers. The behavior of the system near critical points was 
examined to show that transitions between states could occur. For other 
values of forcing the trajectory in phase space formed by the expansion 
coefficients is closed, surrounding two steady states (vacillation). For 
still other values of forcing the trajectory in phase space is chaotic. 

Veronis (1963) performed a similar anaylsis for a barotropic, wind- 
driven ocean in a square basin, with rather realistic results considering 
the approximations made. Stable steady states were found. Under some 
conditions of forcing, the system goes into a limit cycle. 

These studies suggest that even though the forcing is fixed, a system 
can execute transitions from one regime to another. This vacillatory 
behavior is suggestive of the variability of seasonal climate. 

Charney and De Vore (1979) applied the methods to an equivalent baro- 
tropic model forced by orography and a vorticity source, using a channel 
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model with high truncation. For certain combinations of forcing, multiple 
equilibria are found. Three steady states appear for reasonable values of 
forcing although only two of these are stable. One solution suggests a 
"blocked" flow pattern, while the other corresponds to more normal condi- 
tions of small-amplitude waves. This work is suggestive of some atmo- 
spheric conditions although lateral boundaries, truncation and the baro- 
tropic model limit the application of these results. 

Hart (1979) was able to obtain a set of Identical equations using an 
asymptotic expansion method. This result and others indicate that the 
results of Lorenz and of Charney and DeVore are not a fortuitous consequence 
of a particular truncation procedure. 

Furthermore, Ashe (1979) has found a steady state for the Lorenz four- 
variable model in spherical geometry, using modes symmetric about the 
equator. However, Ashe did not investigate the existence of multiple steady 
states nor the nature of solutions near equilibrium. The forcing function 
was obtained from observations, and there was some success in predicting 
the mean fields. 

Charney and Straus (1980) have extended the work of Charney and DeVore 
to a baroclinic, two-layer channel model. Multiple equilibria were confirmed. 
The orography is necessary for the equilibria to exist, but the energy 
comes from the potential energy of the mean flow. Blocked flows, however, 
require unrealistically large thermal forcing. This may he due to the 
extreme truncation, to Cartesian geometry, to quasi-geostrophic dynamics, 
ot to idealized forcing. 
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Spherical geometry eppeers to be mandatory In the study of planetary 
waves. Somerville (1980) shove that the use of a global model produces 
a significant Improvement In the predictability of the long waves, even 
compared to a hemispheric model. A study of the observed global anomalies 
by Wallace (1980) and Horel and Wallace (1980) finds that tcleconnection 
patterns in the 500 mb. geo^otentlal fields are global in extent. A 
dynamical experiment by Hoskins (1978), who investigated the response of 
a linear barotropic model to long-term forcing showed similar results. 

The applicability of quasl-geostrophic dynamics to the long waves 
was reviewed by Kasahara (1976), who notes from scaling arguments that 
the divergence is of the same order as the vorticity for long waves. 
Kasahara suggests the use of the primitive equations, but an alternative 
approach is also possible. Lorenz (1962a) has simplified the primitive 
equations based on energetic consistency, rather than on scale analysis. 
This set of equations is superior to the quasi-geostrophic approximation 
for planetary wave studies. 
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3* 0 Modil Equations 

The equation* of the model are those of the two-layer energy-conserving 
model presented by Lorens (1960). The model consists of prognostic equations 
for i> and r, representing respectively the stream field and its vertical 
shear, and 0 and o, which represent potential temperature and its lapse 
rate. A diagnostic equation is employed to obtain the divergent part of 
the motion field. We employ the notation of Lorenz (I960). 

The model levels are numbered from 0 to 4 with 0 representing the top 
of the atmosphere (p ■ 0) and 4 representing the ground (p - p„ * 1000 mb). 
Variables are subscripted according to which layer they refer. With this 
notation, the winds at the upper and lower layers are, respectively, 


^ = !K xVIlp + T) + O.l) 

- IK X V ( tp'E ) ~ (3.2) 

The corresponding potential temperature fields are 

^ * & t(T 0.3) 

= d~ & • (3.4) 


Thus, i/' and 0 represent the stream field and temperature, respectively, 
at the middle level (500 mb); x and 0 represent their vertical shear. In 
terms of these variables, the equations for x, 0, and o are: 

— ( v Vj = - J’ (ty, t (3.5> 
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= -J(r y f 3 r T /) t J(^ t v 7 r) _ ^ 7 J r 

v f 2 vV + v.(7VX)_' ?v */ 7W (3 -‘> 

Tt = --r< ^<?)-cr fr(<rJr V*-» r «V*ftfr»*j (; j. 7) 

-J7 ^«r; - t; *> * h 2 (e-&*) -+L) fC r f (3>8) 

C^t A 

The constant* and account for frictional dissipation (exclusive of 
Ekman layer forcing which is included in (3. 9)); and the constants h 2 and 
h^ for heat transfer. The function 0* represents the surface temperature 
and therefore controls the magnitude and phase of dlabatic heating. 

These equations will be solved on a sphere of radius a. Therefore, there 
are no lateral boundary conditions. The vertical boundary conditions will 
be specified as: 

(a) At the upper boundary, the vertical velocity is zero. 

(b) At the lower boundary the vertical velocity is the sum of that 

caused by Ekman pumping and of flow over opography of prescribed height 
V This implies that the term x n is specified by: 

<3.,, 

i/2 

Where D is a parameter representing the depth of the Ekman layer, (2K/f) 


A diagnostic equation for the velocity potential x is obtained from the 
thermal wind relation 


h c p = y 0 ( /V V) 


(3.10) 


Where b - .124 


Froa equation* (3.6), (3 .7), and (3.10), wa can derive an aquation for x* 

V.(FV)C V, fv(v 2 (.v.tv?)) = 

( 3 . 11 ) 

i ie f>' v*V f v (Y v _I c <Vc; * <>1 ) 

eft 

Equation* (3.5) through (3.8) arc solved prognostlcally j (3.11) provide* 
the equation* for x* 

3 . 2 Method of Solution 

Equations (3.5) through (3.8) and (3.11) are solved by the spectral 
method. The prognostic variables are expanded In a series of spherical 
harmonics. When the series expansions are substituted into the original 
equations, we convert the partial differential equations into a set of 
ordinary, coupled nonlinear differential equations, one equation for each 
coefficient in the series expansion. In practice, the series must ba 
truncated at some point. The method is well-described elsewhere (Platzman, 
1960, Merilees, 1968) . 

In recent years, the transform method has been used to compute 
nonlinear interactions because of its computational efficiency. The current 
model uses the interaction coefficient method (Merilees, 1968). This is 
done for a number of reasons: 

(a) In systems with few modes (low-order systems) the computational 
panalty is not excessive. 

(b) The interaction coefficients are essential to analytical work. 

(c) The interaction coefficients can be used to filter out certain 
types of interactions, for instance, wave-wav* interactions. 

Thus, it is possible to linearize the model without reprogramming. 



The system of ordinary diffarandal aquatlona which raaulta fro* tha abova 
procadura must ba integrated nuaarlcally. In tha currant varalon of tha 
program, a fourth-ordar Runge-Kutta-Gill procadura takan froa tha IBM 
Sclantlflc Subroutlna Packaga la uaad. 

Tha modal haa baan programmed In auch a way that any combination of 
modaa may ba apacifiad. Tha modal uaar la frea to aalact triangular 
truncation, trapezoidal truncation, or almply a collection of modea. 


3. 3 Model Verification 

A number of preliminary runs have been made with the spectral modal 
described in tha preceding sections. These runs are by no means exhaustive, 
as their purpose is to serve as a demonstration of tha feasibility of tha 
model. The runs all use the same truncation, time step, and initial con- 
ditions; but the forcing is different in each run. We now will describe 
these common elements, in addition, a number of ographic resonance experiments 
have been performed. These are described in section U. 


3.3.1 Truncation 

The modes selected are given in Table 3.1. (called tha "12 mode 6et" 
hereafter). Tha numbers are, respectively, the zonal index l and the longitudinal 
index m of the spherical harmonic 


y 


m 

t 


tl 


v 


p 


Ol 

e 


(f,X)C 



(3.12) 


Whe re N , is the normalization factor (Platzman, 1960) 
ra,l 

2 

The modes selected are a "minimum triangular truncation” plus Yj. This 
mode was included (a) to give some detail to planetary wave 2, and (b) 
to provide a nonzero result for the operator yielding Yj. It 
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It should b« noted that the components of ths stress field (♦, t) ere odd 
for symmetric flow (with respect to the equator, whereas the components of 
temperature are even under the same conditions. Therefore, equatorially 
symmetric conditions require both even and odd modes to be present. 

Table 3-1 

Modes Used In all Runs 


Zonals Sectorlals Tesserals 

°» 0 1.1 2 , 1 

1. 0 2, 2 3, 2 

2. 0 2, -2 3. -2 

3. 0 1 , -i 2 , -1 


3.3.2 Initial Conditions 

The zonal components of the streamfunction t|» and its shear t 
wer® chosen so as to yield a reasonable value of the zonal wind at the upper 
and lower levels. The model then calculates an initial temperature field 
0 from the thermal wind equation (3.10). The static stability fiexd (a) 
was initially chosen to be a constant 30K, representing a vertical strati- 
fication of 30K/500mb in the middle troposphere. Note that the stability 
is not constant, as the model employs the prognostic equation (3.8) to 
calculate stability. 

3.3.3 Time Step and Run Duration 

The time step was chosen by experimentation as a compromise between 
accuracy (Judged by energy conservation properties) and computational load; 
it was fixed at 1.0 hours. It should bo mentioned thac the Runge-Kutta 
algorithm will eventually be replaced by a more suitable mathod. It is used 
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f 

because it Is self-srartlns, stable and accurate, but it is time consuming and 

requires a short time step. Run duration for ail verification runs was 20 days. The 

constants K^, h^, and h^ were chosen to give a relaxation time for 

the model of about 5 days. The run duration is thus a factor o Z 4 larger 

than this relaxation time and permitted many runs to be made with a 

reasonable expenditure of computer resources. 

3.3.4 Results 

A series of runs were made using the initial conditions described 
in 3.3.2. These runs differ in the forcing used in each case, described 
below: 

Case 1. Unforced (Benchmark run). No friction, no heating, 
no orography. 

Case 2. Orography only. The orographic heights are from Tung and 
Lindzen (1979) 

Case 3. Heating and Friction only. 

Case 4. Heating response runs. 

Case 5. Heating plus orography (Phase Plane runs). 

Based on the output of the unforced run, it was decided to analyze 

the output as a time mean plus a transient component. This follows 

standard practice in general circulation and permits one to examine the 

effects of forcing on the stationary components. Presumably with the low 

wave numbers employed there is little if any baroclinic activity so that 

the model is essentially barotropic in the behavior of the low wave numbers. 

In the absence of forcing we should observe, basically, free Rossby-Haurwitz 

M 

, with phase speed -2nM/L(L+l) for the coefficient of iIl . In the lowest-order 


waves 


mode* ^•*^2* ^2 rhi * i8 ex * ctl y wllat wa *••» with «n»ll amplitude 
and phase speed modulations. These waves will have almost no stationary 
component. On the other hand, we expect orography or diabatic heating 
to set up some kind of stationary pattern. This is in fact observed in 
the runs. 

In addition, phase-*plane diagrams are obtained for some 2-dlmensional 
subspaces of the 12-dimensional phase space representing this model. The 
results of the runs are described below. 

(1) Case 1, Unforced Case 

In this case we can observe the appearance of the abovemen tioned 
Rossby-Haurwitz waves after a few days, following the adjustment period 
in the model. The run verifies conservation of energy and enstrophy to 
better '.an . 1 % in the- 20-day period. No stationary components develop. 
Amplitude and phase speed plots for the nonzonal modes versus time would 
be a straight horizontal line and are not reproduced. 

(2) Case 2. Orography Only 

Orography is introduced through the term in psi and tau equations 
(35) and (3.6), of the form 


L 7 , (fV%o) 
2 


(3.13) 


where y Q is obtained from equation (3.9). In these runs, the effect of 
the Ekman layer is not included, i.e. D * 0 In (3.9). Orographic forcing 
then reduces to 

± £ V* ( fv f (h S) ty-V)) (3.14) 

2 


where h is the dimensional height of orography and C is a constant, 
s 

The value of this constant was -5.0 x 10 ^ in a »uns. The order 
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of magnitude of the coefficients of h is 400 m. There is a nonzero 

8 

coefficient for each mode within the truncation limits. 

The effects of orographic forcing oc M« modes and ^ 2 are 
most pronounced (Figures A. 2. a and b) . A stationary part appears. The 
magnitude of the stationary part is about 20Z of the maximum value of 
and 5X of the value of These coefficients now exhibit both ampli- 

tude and phase speed modulations, with the maximum of phase speed often 
coinciding with the minimum of phase speed. 

(3) Case 3. Heating and Friction Only 

Diabatic heating is modeled with a term proportional to the tempera- 
ture difference relative to a fixed "ground temperature" ©*, in equations 
(3.7) and (3.8). In the run described here, the function 0* was chosen 
to be identical to 0 (t * 0). Thus, the initial contribution of the 
heating term is zero. As the temperature 0 and a changes, the temperature 
differential increases (or decreases) and tends to restore 0 to its 
initial value. The constants h 2 and h^ were chosen to be 10 6 sec 1 . 

Friction is introduced into the and T equations (A. 5) and (A. 6) 
by linear terms. The proportionality coefficients used are K 2 - - 

10 ^ sec 

1 2 

The effect of heating and friction in <|>^ and ij> 2 are shown in Figures 
3. 3. a and 3.3.b and on © 2 in Figure 3.3.c. For comparison, the unforced 
behavior of 0* is shown in Figure A. 4. Except for the initial, transient 
behavior the effect of heating in and \p 2 is small. The effect on the 

temperature mode © 2 is to smooth out the fluctuations in the amplitude 
and to introduce nonlinearities in the phase curve. These effects are 
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reasonable, inasmuch as the thermal forcing Is very small and not chosen 
for realism. 

(4) Case 4. Heating Response Runs 

2 

In this case, the heating was Introduced only through the mode © 2 *, 

i.e. all other modes of 0* were set to zero. Twelve separate runs (three 

seta of four runs each) were performed for this case. In each run the 

2 

amplitude and phase of 0^* remains constant throughout the run. However, 

2 

the amplitude and phase of the forcing 0 2 * changes from run to run. Each 

2 

set of runs holds 0^* at constant phase and varies amplitude. For the next 

set, the phase is changed from the previous set's value and the amplitude 

2 -2 

varied. The stationary component of 0 2> 0 2 , is calculated for each case. 

The results are summarized in Table A-2. As one would expect from linear 

—2 2 

forcing, the stationary component © 2 approaches 02*" l ar 8 er the ampli- 

tude of forcing, the larger J^j* The P hase adjustment of 0* after 20 
days depends on the amplitude of forcing; for large amplitudes of the 
forcing, the phase of the stationary component almost coincides with that 
of the forcing. 
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Table A-2 

Response of Stationary Component to Forcing 
(After 20 Days) 


0* * 

Arg. © 2 * 

5 '! 

► 

00 

<DI 

2 

-45° 

2.09 

-55.4° 

7 

-45° 

3.88 

-51.5° 

14 

-45° 

6.40 

-49.8° 

21 

-45° 

8.91 

-49.0° 

2 

45° 

1.41 

-30.5° 

7 

45° 

2.59 

11.5° 

14 

45° 

4.92 

26.9° 

21 

45° 

7.38 

32.6° 

2 

135° 

0.72 

-72.1° 

7 

135° 

1.20 

147.5° 

14 

135° 

3.69 

137.2° 

21 

135° 

6.21 

135.7° 
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(5) Cast 5. Heating and Orography 

In this run, the effect* of heating and orography were combined. 

The results are presented in the form of phase subspace plots in terms of 

the amplitudes of selected coefficients. In Figure 3. 5. a we see the 
2 2 

amplitudes of and plotted with time as a parameter. The model is 

obviously approaching equilibrium, possibly a tight limit cycle as can be 

seen from the spiral behavior of the phase plane plot. 

2 1 

Figure 3.5.b shows similar behavior for (ij^, ) subspace. 



Figure 3.1. a The initial zonal wind field at the upper and lower level 
as a function of latitude. 
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time (DAYS) 


Figure 3. 2. a The response of the nonstationary part of amplitude and 

phase speed to orographic forcing only. The stationary part is 
7 2-1 

1.02 x 10 MS at a phase of 309°. The barotropic phase speed is 
180°/halfday . 
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Figure 3.2.b The response of the nonstationary part of amplitude and 

phase speed to orographic forcing. . The stationary component is 
7 2-1 

1.12 x 10 MS at a phase of 324°. The barotropic phase speed is 
120°/halfday . 




Figure 3. 3. a The response of the 
phase speed to heating. The 


nonstationary part of ** amplitude and 
stationary part is negligible. 


Amptoud*, 



figure 3.3. b The response of the 
phase speed to heating. The 


nonstationary part of v* amplitude 
stationary part is negligible. 


and 
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Figure 3.3.C The response of the amplitude and phase of the nonstationary 
2 

parr, of to heating. Solid lines labeled 1 and L are the amplitudes 
of the forcing and of the stationary component, respectively. The 
broken lines labeled 2 and 3 are the phases of the forcing and of 
the stationary component, respectively. 
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4.0 Resonance Experiment! 

4.1 Motivation 

The existence of multiple equilibria in approximate, but non* 
linear, barotropie systems (Charney and DeVore, 1979; Hart, 1979) is 
intimately connected with the existence of a linear resonance condition. 

In fact, the equilibrium associated with a large wave amplitude "blocking" 
flow is located close to the linear resonance point. Consequently, as 
we will show below, the study of the resonance condition is important 
for understanding multiple equilibria. 

However, it is instructive to begin with a much simpler system 
which exhibits this type of equilibria, even though it will be in an 
ad hoc way. Consider the barotropie vorticity equation with forcing 
on a beta-plane: 


^ (J/ 

at 




(i) 


where V is the stream function and F some as yet arbitrary forcing 
function. We can find a solution to (1) on a (3-plane of the form 
(zonal flow plus a sir. b le wave): 

V' = - 0u +- y> '(*,&) <2) 

Which reduces (1) to the linear partial differential equation 


ok v 






-h i~ 


(3) 
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The plane-wave solution for the homogeneou. equation 1. 


With c given by the usual Rossby wave dispersion formula. 

c= 


(5) 


We can find a particular solution by setting SH'/St - 0 m eq. (3). 

Then, if ^is the particular solution 

# = J— 

* (6> 

and wa ... hots that, as k » 8 /U aquation (6) becomes rasonant , Thi. 
is a well-known condition discussed, for cample, in Dickinson (1980) 
with respect to orography. In general the solution to (3) is the 
su« of the stationary wave plus a traveling wave with phase speed given 
by (5). Resonance corresponds to the case where the traveling wave 
slows down, and interferes constructively vi -h the stationary wave (as 
suggested by Tung and Lindzen, 1979). 

The unbounded response predicted by eq. (6) is not expected to 
occur in the atmosphere. (Note that the inclusion of friction would 
prevent the unbounded resonance in eq. (6)). Instead, nature will act 
in such a way as to prevent the resonance. This can be accomplished 
by non-linear interactions between the different waves which will be 
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present which will modify the zonal flow away from resonance. Thus 
there will be a nonlinear wave-zonal flow feedback mechanism inhibiting 
resonance. Therefore, we introduce an ad hoc linear feedback 
mechanism between the single wave of eq. (2) and the zonal flow, i.e. 

0 - A 9o + (3 rn 

then we can have multiple equilibrium solutions to the system (6) and 
(7) as shown in Fig. 4.1. Here, the solution of the system (6), (7) is 
given by the intersection of the straight line (7) with the resonance 
curve (6) for a particular value of k. 

Fig. (4.1) has been drawn with B>0 and A>0 such that there are 
three possible solutions (1,2,3). These correspond to the multiple 
quilibria found by the previously cited authors, who, however, have 
introduced more than one wave into the problem and physical mechanisms 
representing the ad hoc constants, A and B. The existence of the 
equilibria depends on the values chosen for A and B. 

Thus, the existence of multiple equilibrium solutions to forced 
barotropic (and presumably also to baroclinic) models is closely 
connected with the existence of a resonance condition in a forced wave, 
and to a feedback mechanism of the wave with the zonal wind. It is 
of interest, therefore, to use the model described in the previous sections 
and to attempt to induce a resonant state In this model, both for the 
baroclinic and barotropic. conditions. This work is the subject of 
the following sections. 
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4.2 Barotropic Reaonanca Experiments 

W« consider the 12-aod. a.t de.crib.d in ,.ctlon 3. I„ t„l. ,. t of 

■Oden. w. restrict oora.ive. to the berotropic c ... , T . o . s . 0) 

.tapllfy th. 8.C by assuming hemispheric lyi » atry . ni . acco „ pluh , d 

by setting eil even node, of no aero. Tb. training mo de. are given 
below • 

Mode number 1357 

M 2 10 0 
L 3 2 3 1 

We retain the node (2,2) for the orography, if ^ . i^ei, then t „. 
following eigenvalues are equal: 2 


10 

-1 

2 


12 

-2 

3 


C 1 ’ c 5 1 
C 3 * c 10 


“12 


The spectral form of the barotropic vorticity equation (1) is the 
following set of ordinary differential equations. 


(8)a 

(8)b 


= ^// h 2Si ‘-2 % +JL, % 


(9) 


% ” hsz +■ ?3?3 ^ 3^7 i- f-k 


3 h Z ^ 


( 10 ) 


* 

^6 - Hjj I2 f r h 2 if,- + l-l 1L > J L J ^ 


(ii) 
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( 12 ) 


Where is a mode in the stream field (see table) and 

■V - ^ v 

c 4r 


(13) 


C L - t-i Ul +L ) 


is the eigenvalue of ith mode (14) 


J ij k 18 the JacoblaR Interaction coefficient, or "coupling integral." 


^ - 


- 2 JZ. Ml 


( L < +L) 

L ‘*e Rossby-Haurwitz phase speed. 


(15) 


H ijk is the coefficient of orographic interaction. The notation 
— Mn 

>Fn means (-) fn*. In general, (8)-(ll) would contain other terms 
(see e.g. Kallen (1979)) but these drop out because of the Eigenvalue 
relations (8). Thus, in this case we have no wave-wave +zonal flow 

interactions. Of course, in a two-wave set we cannot have wave-wave-wave 
interactions. 
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W« can obtain a resonanca condition by linearizing (9)-(12), 
assuming the zonal flow is constant: 


t = % 


( 16 ) 


h = <!>, 


( 17 ) 


We can then solve for a steady state, analogous to eq, (6), by taking 




at 


1 ■ 0 in eq. (9). The result is 


Where H, 


^7 ~~ ^ 

T in 


(18) 


H 251 h 2 * 


Equation (18) is the analog in a spherical spectral system of eq. (6). 
Since 4^ is the earth's angular momentum, it is very unlikely that 
we will ever encounter this particular resonance in practice. However, 
a different set of modes would exhibit similar behavior (Kallen, 1979). 
The resonant value of 4^ (zonal flow) is given by eq. (18): 




1, 


ih 


A- 6^* 2 y tO 
4T 


(19) 


using a reasonable value for the orographic mode h 2 (Tung and Lindzen, 

1979) we can plot eq. (18), which is shown as the broken line in Fig. (4J.), 


It li interesting to observe the behsvior of the system when the non- 
linesr aquations ere used. Accordingly, « series of experiment, were 
performed with the model to determine the shape of the resonance curve 
for the nonlinear case of eq«. (9)- (12). Since the numerical model 
cannot solve for the steady state, the steady state must be found by 
averaging. Fig. (4.2) shows the results of these experiments. The 
solid curve represent, found by averaging over a period of 30 days. 
The ordinate is the initial value of 'f ? , which for practical purposes 
is the average value of since it does not change very much over the 
run. It can be seen that the response curve is broader than that 
predicted by (18) but that resonance occurs in about the right place. 

As predicted by (18), the traveling portion of the disturbance was found 
to slow down as resonance was approached. 
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4.3 Baroclinic Resonance Experiments 

We consider the effect of baroclinicity in the resonance conditions. 
We will use the 12-mode set, restricted to odd modes in f and t and even 
modes in 0 and X . The stability o i. taken to be constant. In order to 
simplify the analytical work we will consider only the f and t equations. 
This will bypass analysis of the "omega equation” for the variable X ; 

whi-h is not -^gorously justified. However, the principal results will 
be obtained anyway. 

Using the same mode designation as before, we have the following 
equations : 


f 


T .l = 


- r /7/ ft ^ + J th U ' 5 < ( St, ^ (20) 

r iv VC} t r i?i Cf ^ -h i fs/ r, <P S t £ s „ +2,1;, c 2 j if ( 2 i) 

V (Vi-T/) + W <»> 


r 3 - H 


' 7 - 


2,10, 3 6} * } + i&jVz 

= ree (xl) 


7 zz 


5 = 


5 » 


;ce (tz; 

/ 


123) 

(24) 

(25) 

(26) 
(27) 


The linearized analysis proceeds by taking the zonal flow and its shea 
to be constant: 


t.< ** s 
% - % 


XT _ T 
S ~ L S 

^7 - XL 


(28) 

(29) 
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Th« aquations for ths highest order wave node (mode 1) are 

% = ( r <7/ % + I /7/ ^ i\ t U ut (ty s -T s ) oo) 


• ^ 

= ( I n , t-JL, ) L, + 


U U, <-h-f s }n) 


We will now assume that the term x 2 represents thermal forcing, denoted 
by x*J while the orograhic term will be denoted by H*. Also letting 



(*£/?< W + *£/ ) = v h 


(32a) 


(32b) 


We have the following set of equations for the wave: 

<1, = %% + r z v, + h* 

“t; = Tz % *■ k v , 


(33) 


( 34 ) 
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Th« steady state solution (cf. (18) and (6)) Is (Ivan by: 

fi = - H'oht' r 7 )- 
v = 

I ik a " 


(35) 


(36) 


The resonance condition is then 



(37) 


Whereas in the baiotropic case yO is the resonance condition (from (13) 
and 32a). Since for a westerly positive stream flow i z is negative, a 
higher value of f g is needed to obtain the resonance than in the baro- 
tropic rase. That is, the zonal wind required to obtain resonance in 
the baroclinic case is higher (for negative t z corresponding to wind, 
increasing with height). 


A series of numerical experiments with the numerical model were 
performed to examine the nonlinear resonance case. The experiments 
consisted of selecting values for mean wind and shear t ? and finding 
the steady state ^ by averaging over a 30-day period. The results 
are shown in fig. 5.3, which displays the response of the wave f for 
va-ious values of and T; . The results are in agreement with equation 
(37); we note that as shear Increases, a correspondingly high value of 
^7 is needed to cause a resonant response. 
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fig 4.1 Schematic illustrating multiple equilibria obtained by solving 
equations 5.6 and 5.7 simultaneously. 
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2 

fig 4.2 Amplitude of the stationary component of as a function of 

0 

the zonal flow for the barotropic ca9«. Dashed line la 
linear case, solid line is nonlinear result. 



5.0 Cone lus Iona 


The model described in section 3 has been verified by demonstrating 
that the model yields reasonable results. In the absence of any forcing, 
and with reasonable initial conditions, we find free Rossby waves 
propagating with the correct phase speed. When orographic forcing is 
added, a stationary component develops. .The noristatl onary component has* 
an instantaneous amplitude and phase speed which changes with time. 

This is to be expected from the interactions with orography. When 
diabatic heating is added, the response of the stationary temperature 
field is roughly linear. 

The phase plane plots show an approach to equilibrium as one would 
expect from constant forcing without feedback. Further judgements, 
based on the limited amount of runs performed thus far and on the lack 
of verisimilitude of the thermal forcing, would be unwarranted at this 
time. 

The barotropic and baroclinic resonance experiments are intended as 
preliminary steps in understanding the mechanisms underlying the 
existence of steady, or quasi-steady, states in the atmosphere. The 
theory of resonance for the cases covered :.n the existing literature 
(zonal flow and two waves), are verified and extended by the model. 

The effect of shear on the resonance curves is also investigated. 

For positive shear (t) the zonal wind required to obtain resonance is 
very high. 

Further work in this direction will proceed along the lines of 
investigating more rea. _tic mode structures both in zonal flow and 
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In Che wave atructuree. Considering first Che barocropic case, Che 
addition of more waves raises che possibility of wave-wa\ e-wave inter- 
actions, some of which may also produce resonance leading to stationary 
situations. A guide to this possibility is the numerical results of 
Egger (1978, 1979) who has simulated blocking situations by wave-wave 
interactions in a barotropic, forced, channel mode. The extension of 
this work to spherical coordinates now underway. Further experiments 
with both barotropic and baroclinic analogs of Egger' s experiments is 
planned. Also, analytical work based on these experiments is beinpc 
carried out. 

The results of this study should have an important bearing on both 
the atmospheric response to imposed forcing and feedback mechanisms, 
but also on the nature and existence of multiple equilibria. But these 
two topics are intimately connected, as we have shown. 

Thus we are addressing the fundamental cause of planetary-wave 
generation and their (nonlinear) interaction with other waves and the 
zonal flow to form occasional, smaller scale "blocks" in the atmospheric 
circulation. The role of multiple equilibria in this process involving 
several scales is of great importance. 
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6. Multiple Flow Equilibria and Chaotic Behavior 

We have undertaken an Investigation of the occurrence of multiple stable 
equilibrium states in simple but heuristi cally valuable models of large-scale 
atmospheric flows. Our aim 1 b to expand upon the recent results along these 
lines of Charney and his colleagues (1,2). They have shown that multiple stable 
equilibria do occur in highly truncated spectral models of both barotropic and 
baroclinic flows which are topographically or thermally forced at the lower 
boundary. By numerical integration of an associated grid-point model for the 
topographic forcing, having many more degrees of freedom than the spectral model, 
they have verified the presence of these equilibria. The equilibria can be 
viewed beuristically as being capable of undergoing transitions among themselves 

due to small-scale perturbations. From this viewpoint, one of the transitions 

* 

bears a striking resemblance to the often observed shift of the zonal westerlies 
to persistent blocking patterns. 

Our work centers on several of the difficult mathematical issues involved 
in analyzing the dynamics of spectral models like those considered by Charney 
and his colleagues. Ultimately, one would like to characterize the attractor sets 
of these dissipative, nonlinear dyn; ^lical systems. The determination of equilib- 
rium states entails solving relatively large systems of polynomial equations that 
express the algebraic conditions lor an equilibrium. The number of equations 
and unknowns equals the number of degrees of f reedotj of the spectral model. The 
only known systematic solution method for finding all of the Isolated solutions of 
a polynomial system is the so-called homotopy continuation method (3). We have 
developed a computer implementation of the Li-Yorke homotopy continuation algorithm 
(4) to provide an appropriate numerical solver for our investigation. Our aim is 
to uoc the computer capability to find all of the equilibrium states for Chamey-like 


models, especially for cases witn relatively large numbers of degrees of freedom. 

Our solver has been tested on a Charney model having six degrees of freedom. 
We have thusly rediscovered all the equilibria indicated by Charney and DeVore 
in (1) for selected model parameters. In (1) it was shown that the problem of 
solving the particular polynomial systems that arose was reducible to sampler 
subproblems which enabled the authors to find all of the equilibria in a finite 
region of the dynamical phase space. We are as yet uncertain if other equilibria 

exist outside this region and are pursuing the question. The next application of 

I 

our solver will be to the baroclinic models introduced in (2) where twelve degrees 
of freedom are involved. In this case it is almost certain that the equilibria 
already found do not exhaust the totality. 

Other efforts under consideration include (a) applying the solver to the 
global baroclinic model of Rodenhuis, et al. and (b) applying recently discussed 
techniques (5) for estimating the dimension of the dynamical attractor set to the 
most meteorologically interesting of the spectral models studied. 
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VIII. A Simulation of the January Standing Wave 
Pattern Including the Effects of Transient Eddies 

by 

J. D. Opsteegh 


A. D. Vernekar 


Abstract 


A steady-state linear, two- lava 1 primitive aquation 
modal is usad to simulate the January standing wave pattern 
as a response to mountain, diabatic and transient eddy effects. 
The model equations are linearized around an observed zonal 
mean state which is a function of latitude and pressure. The 
mountain effect is the vertical velocity field resulting 
from zonal flow over the surface topography. The diabatic 
heating is calculated using parameterized forms of the 
heating processes. The transient eddy effects, that is the 
flux convergence of momentum and heat by transient eddies, 
are computed from observations. Separate responses of the 
model are computed for each of the three forcing functions. 

The amplitude of the response to diabatic heating is 
small compared to observed values. The vertical struct' re 
is highly baroclinic. At the .upper- level, the phase of the 
waves is approximately in agreement with the observations. 

The amplitude of the response to mountain forcing is comparable 
with observations. The wavelength of the response in the 
Pacific sector is shorter than observed. The vertical 
structure is equivalent barotropic. The combined response to 
diabatic heating and mountain forcing is dominated by the 
contribution from the mountains. The phase shows some agree- 
ment with the observations, but the Aleutian low is located 
too far to the west and an unrealistic high appears to the 
west of the dateline. 


The amplitude of the response to transient eddy effects 
is comparable to the observations in middle and low latitudes. 
At high latitudes the amplitudes are much too large. The 
assumption of linearity is not valid for strong forcing at 
high latitudes where the zonal wind is very weak. The vertical 
atructure of the response is almost equivalent barotropic. 

A comparison of the responses to mountain and transient 
eddy effects show an interesting phase relationship. The 
troughs produced by the transient forcing are found in the 
lee of the troughs produced by the mountains (very close to 
the ridge) indicating that transient forcing is organized 
by the mountain effects. 

The combined model response to all three forcing functions 
shows a good agreement with observations except at very high 
latitudes . 


A Numerical Simulation of the Influence 
of the Hadley and Ferrell Circulation 
on Forced Stationary Planetary Waves 


by 

J. D. Opsteegh 


ABSTRACT 


Thi influence of the Hedley end Perrell circulation on planetary waves* 
resulting from a local heat source in middle latitudes* the subtropics 
and the tropics is investigated with a linear steady-state primitive equation 
model. The model has been linearised around the observed sonally averaged 
January mean state. This mean state includes the mean meridional circulation 
(MMC). Both simulations with and and without the MMC have been performed. 

Zt is shown that incorporation of the MMC causes a shift in the resonance 
wavelength towards waves with lower sonal wavenumber. The inclusion of a 
small amount of friction has a very significant influence on the quasi -resonant 
waves. However* it appears that the solution is not sensitive to the particular 
value of the friction coefficient. 

As a consequence of tropical heating a strong Valker circulation can be 
observed with large meridional gradients in zonal momentum. This means that 
the advection of zonal momentum with the mean meridional wind must be an im- 
portant process and therefore the inclusion of the MMC has a significant 
influence on the tropical Walker circulation. However* the remote response 
at middle and high latitudes is only slightly damped and the phase of the 
mid-latitude planetary waves does not change. 

Heating in the sub-tropics generates a strong jet in the tropics. How- 
ever the structure of this jet is completely different from a Walker ciruclation. 
' it is much more geostrophically tied to the pressure pattern. The Hadley 
circulation has again a large influence on the strength and position of this 
jet. But in addition, also the mid-latitude pressure pattern changes sig- 
nificantly. 
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X. Infrared Radiative Transfer Through A Regular 
Array of Cuboidal Clouds 


by 

Harshvardhan 
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A. INTRODUCTION 

There is s need for an adequate radiative parameterization of clouds 
in clinate and general circulation models. To data, all models have 
utilised plane parallel theory to construct parameterisation schemes 
although a significant portion of the total cloud field is broken. 

The particular case of monochromatic radiative transfer through isolated 
cuboidal clouds has been studied in both solar and infrared wavelenths 
(McKee and Cox, 1974» Davies, 1978| Liou end Ou, 1979, and others). 

However, in any broken cloud field, each element can not be considered 
independently as there is cloud -cloud interaction and for incident solar 
radiation, there may be partial shading of neighboring cloud elements. 
Ellingson and Kolczynski (1980) have used the technique of Avaste (1969) 
to compute heating rates in the atmosphere in the presence of an array 
of black cylindrical clouds. 

The present study considers a regular array of identical cuboidal 
clouds, as in Aida (1977), overlying a non-reflecting surface. Only 
infrared radiative transfer is considered and the theory of transmission 
through the array and emission from the cloud field is first developed 
for black isothermal clouds. The work is then extended to lCpn radiation 
interacting with a cloud field having prescribed radiative properties. 

For this, the two stream solution to the isolated cloud problem (Harshvardha n 
et al., 1980) is used in conjunction with the black cloud results. 
Computations are made for both fluxes and radiances as a function of 


the cloud fraction 
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I RADIATIVE TRANSFER THROUGH A BLACK ARRAY 


Consider an idealised broken cloud field in the form of an extended 
regular array of cuboida. A portion of this array ia Illustrated in 
Fig. 1 , from which it follows that the area of the top and hottoa face 
is Af — s 2 and of each side face, A } » ss*. Me define the aspect ratio 
of each element as 


*• 

a » - 
s 




(l) 


If the spacing between each element is equal in the x- and y- dirwetions, 
then the fraction of the field covered by clouds when viewed normally, 
the cloud fraction, is 

f 




( 2 ) 


For an isothermal cuboid, the monochromatic power emitted by the top 
face is itboAd and aach sida face, itBqAi where B 0 (T 0 ,1) is the Planck 
function. Half the energy emitted by the side faces escapes to space 
when the. cloud is Isolated. When the cloud is embedded in a field as 
in Fig. 1, a fraction, F, of the energy diffusely emitted by each side 
face is intercepted by the faces of neighboring clouds. In this case, 
it can be shown that the total power emitted to space by each isothermal 
cuboid is 


P - TO 0 A 0 [l + 2a Cl— F) jj (3) 

We define the power ratio analogous to omittance in plane parallel theory 

as 


*• - War 1 + 2 * tl - p > ,4 > 

We also define the effective amittanco of a broken cloud array as the 
power emitted by a region containing the array in the absence of ground 
emission to the power that would be emitted by the region if it were 
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entirely covered by e block surface et the mm teapereture. Zf there are 
N Identical clouds of tops as in tiw pres mt nodal, then it is evident 
that the effective emittanee is 

1 * ST*, * f *> <*> 

because f - Ni^/Ag, where Ag is the total area of the region. The effective 
transmittance through an array of black clouds is given by energy balance, 

T « 1-E (6) 

Thus, the power emitted to space from a broken cloud array at T # over a 
black ground at temperature, T,, is 

P out “ rt »V + V (7> 

and the flux per unit area of the field is given by 

M - —Hi. . TJBjE + TTBjT (8) 

g 

The above equation is similar to that obtained by the plane parallel theory 
with E playing the role of the cloud fraction. We therefore introduce the 
concept of effective cloud fraction, f # « E, such that 

M - VB 0 f # + ITBj <l-f # ) (9) 

1. Effective Cloud Fraction 

The computation of t m essentially involves the computation of F as 
they are related through eqs. (4) and (5). Alternately f # may be obtained 
by measuring the transmission of diffuse light through an opaque array. 

We have attempted both methods, the former by using the theory of angle or 
configuration factors (Sparrow and Cess, 1978) and the latter by setting up 
an array of black blocks on a light table and measuring the flux with a 
Tektronix Model J 16 digital photometer with a cosine corrected remote 
illuminance probe. The probe was moved around and the mean of a number of 
readings was used* to compute f 0 for each configuration. The theory is expected 
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to give good results os long ss shading is not too sorious a problem taut 
as f inorsasss it becomes lnoraasingly difficult to estimats F. Ths rssults 
ars givsn In Pig. 2 with a fit to ths data and thsory givsn by 

, . 6 do) 

a l+2af (1+0.1 5f) 

Ths curvss follow .no gsnoral pattorn of ths rssults prsssntsd in Ohvril 
and Epik (1978) and Piling son and Kolccynski (1980) . It nay ba no tad that 
as f+0, f # ‘*(l+2a)f, which is ths corrsct limit for an isolated cuboid. 


2. Pluses 

A consequence of ths finits nature of each cloud in ths array is to 
increase the effective cloud fraction because f^ is always greater than f . 
This implies that energy loss to space from a broken cloud array is less 
than that computed using the cloud fraction without considering the sides 
of the clouds. The effective black body amission temperature }f a broken 
cloud array at 255°K overlying a black ground at 290°K ia shown in Pig. 3. 
The curve marked a ■ 0 is the flat plate or plane parallel case. Zt is 
evident that errors of 5° - 10°K can be made for clouds of aspect ratio 
unity, which is a common value for fair weather cumulus. 

Furthermore, there is a change in the cooling rate across the cloud 
layer because of the radiation intercepted by the sides of the cloud. For 
a black isothermal cloud in a non-participating atmosphere, it can be shown 
that 


& 

finite - £ e 

(11) 

iSfT 

f 


pi ) 

infinite 



Thus# the cooling rate in the atmospheric window could be two to five times 
that obtained by the flat plate assumption. 
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3. Radiances 

Zf the array ia oanpoaad of black iaotharmal blocks, tha radianca in 

' any direction depends on the field of view obscured by the clouds. Whan 

viewed from on top, only the cloud fraction obscures the ground, whereas 

from any lower observation angle, the sides of the cloud also obscure the 

ground. This effect is shown in rig. 4 in which the radiance at 4*0° 

for an isothermal regular array of cuboids of aspect ratio unity is given 

0 

as a function of zenith angle. The clouds are again at 255 K while the ground 
ie at 290°K. A plot of the brightness temperature when ■ 45° is also shown 
for f - 0.1 and 0.2 illustrating the effect of viewing direction. It is seen 
that the radiance drops off rapidly as more of the clouds and less ground is 
visible. This suggests that the brightness temperature measured by a satellite 
instrument can not be related to the cloud top temperature for broken cloud 
fields without considering the geometry of the array and the viewing angle. 
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C. MONOCHROMATIC RADIATIVE TRANSFER THROUGH A NON-BLACK ARRAY 

Let ut now consider a regular cloud array as in Pig. 1 but relax the 
assumption that tha cloud faces are black emitters. He shall consider 10pn 
radiation and assume the clouds to be water droplets with a Cl size distribution. 
Furthermore, the clouds are isothermal and the ground is non-reflecting. He 
shall try and use the results obtained for the black array to solve for fluxes 
and radiances. 


1. Isolated Cloud 

The radiation field in an isolated cuboidal cloud has already been obtained 
using the two stream approximation by Harshvardhan et al., 1980. It is sufficient 
to note here that fluxes leaving the faces of a cuboid may be expressed in terms 
of a function, Lo, which’is obtained from the solution of a diffusion type 
equation. If the flux leaving the face is designated M, then at the X • ^ 
face, M - M(y,z) and so on. We can also obtain the mean flux leaving 
each face, M. If the cuboid has equal sides, s, then it is obvious 
that (MJjjaSi * (M) ^ *. Using this information we may construct the boundary 

2 yr* j 

conditions to solve the array problem. 


2 . Conditions 

In the isolated cloud problem it is assumed that the side faces see the 
ground whereas if the cloud is in an array, the side faces receive flux from 
neighboring clouds as well as the position of the ground that is visible. 

If the diffuse emittance from the side faces was independent of position, 

mm 

such that M(y,z) - M, then the boundary condition on the side face can be 
shown to be 

(K) s =■ FN + (1-F) Bj 

*"5 2 


(12) 
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where F is the angle factor batviaan tha faca In question and all the other 
faces visible from it# and D, is tha ground amission. 

However , even if we assume tha exiting flux to be diffusa we can not 
assume that it is uniform and tha error introduced by eq. (12) will be con- 
siderable as the cloud fraction, f, increases. In the limit, f+1, we know 
that the input to the side face is exactly the output of the neighboring 
face, point for point. He have therefore chosen to weigh the boundary 
conditions as follows 


(M) • 

x- 5 


0 (M) s + 
X— j 


(i-4>) + (1 -f) 


(13) 


where 4>=*0 when f<f* 

and 1 - wh en f>f* 

f 

with f* * | 


The problem can then be solved as for the isolated cloud using the values 
of F obtained in the previous section. It v,\y be shown that the upward 
flux from this array is 

2a (1-f) (l+0,15f ) f M + fMt + U ‘ f> B, (14) 

_ l+2af (1+0. 15f) l+2af (1+0. 15f) 

where K is the flux out of each side face and Mf is the flux out of the top 

face. 


3. Results 

Computations have been made for a-1 and the radiative model mentioned 
earlier for which the single scattering albedo, u>« 0.64 and asymmetry parameter, 
g ■ 0.86 at X* lOym. Fig. 5 shows the flux escaping to space from a regular 
array of cubes of optical depths varying from 0.5 to 20. Results are presented 
as an effective black body temperature and may be compared with the curve marked. 
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•■1 in rig. 3 to not* th* departure from the black caee. Xt ii evident 
that an optical depth of 10 or greater way be approximated quite well 
by black cuboids. 

Thie ie shown further in rig. 6 in which the radii' n at$ - 0 is 

plotted for various cloud fractions for a regular array of cubes >f 

optical depth 10x10x10. Th* ground and cloud temperatures are as before 

and the plot may be compared with rig. 4 which is the black case. For 

0 

zenith angles less than 60 , the two plots are very close indicating that 
finite cloud fields may be approximated by totally opaque diffusa emitters 
for purposes of modeling the radiative field exiting the array. 

D. Summary 

Finite cloud fields have been modeled by a regular array of cuboids. 

It is shown that the role of the sides of the clouds is very important in 
computing radiative quantities. Using an approximate technique to solve 
che isolated cloud problems, results have been obtained for radiative transfer 
through a participating array of cuboids. One conclusion of this study is that 
clouds may be modeled by black casque emitters if th* optical dimensions 
exceed 10 on each side, which is true in the infrared for virtually all 
cloud fields. 
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FIGURE CAPTIONS 


1. Schematic showing regular array of cuboids. 

2. Effective cloud fraction, f e , as a function of the cloud fraction, f, 
for various aspect ratior. 

3. Upward flux from array expressed as effective black body temperature for 
T e » 255°K and Tj - 290°K. Aspect ratios are marked on individual curves. 

4. Radiance vs. zenith angle expressed as brightness temperature for two 

azimuth angles, <£» o° ( — ) and <J>-45°( ); a=*l, Tj, »255°K, ■ 290°K. 

Cloud fractions marked on individual curves. 

5. Same as Fig. 3 but for non-black clouds of unit aspect ratio. Radiative 
properties given in text. Optical depths marked on individual curves. 

6. Same as Fig. 4 but for non-black clouds of unit aspect ratio at $=*0°. 
Optical dimensions are 10x10x10. Radiative properties given in text. 




Figure 1 


Schematic showing regular array of cuboids. 
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Figure 2 


Effective cloud fraction, f , as a function of the cloud fraction, 

C 

f, for various aspect ratios. 


1 




BLACK ROPY TEMPERATURE °K 



CLOUD FRACTION, f 

Figure 3 

Upward flux from array expressed as effective black body 
temperature for T 0 = 255°K and Tj = 290°K. Aspect ratios are 
marked on individual curves. 
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Figure 4 


Radiance vs. zenith angle expressed as brightness temperature 
for two azimuth angles, $ = 0°(-) and <t> = 45°( — ); a = 1, 

T q = 255°K, = 290°K. Cloud fractions marked on individual 


curves. 
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CLOUD FRACTION, f 

Figure 5 

Same as Fig. 3 but for non-black clouds of unit aspect ratio. 
Radiative properties given in text. Optical depths marked on 
individual curves. 
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Figure 6 


Same as Fig. 4 but for non-black clouds of unit aspect ratio 
at 4> = 0°. Optical dimensions are 10x10x10. Radiative 
properties given in text. 
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1. Introduction 

Doppler radars can now measure directly the details of the horizontal 
wind components. While great uncertainty still exists about the vertical 
components (Miller, 1978), It appears that under certain circumstances 
It Is also possible to deduce that component as well (Carbone et al., 1980; 
Clarke et al., 1980; and Gal-Chen, 1979). The radars cannot, however, 
measure directly the three dimensional structure of thermodynamical 
variables such as density and pressure fluctuations. These latter 
variables are nevertheless essential in studies of convective systems. 

The ideal solution would be to develop instrumentation capable of directly 
measuring thermodynamical variables with a spatial and time coverage 
comparable to that of Doppler radars. Unfortunately, the technology 
for that kind of network does not exist and one must use indirect methods. 
Recognizing the above cited limitations, several Investigators (Hane and 
Scott, 1978; Leise, 1978; Gal-Chen, 1978; and Chong et al., 1980) have 
proposed indirect techniques whereby the combined use of the relevant 
hydrodynamical equations and the observed wine will permit a unique 
determination of the density and pressure fluctuations. 

In this paper the basic methodology for obtaining density and pressure 
fluctuations is outlined. This is followed by a discussion of possible 
approaches to verification and error analysis. The last section contains 
some demonstrations how the method works with real data. Temperature and 
pressure fluctuation deduced from Doppler radar observations of the boundary 
layer (Gal-Chen and Kropfli work in progress) are presented. The results 
of Chong et al. (1980) and Gaillard and Gill et (1980) who used Doppler wind 
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data to estimate the thermodynamical structure of weak convective cell 
generated by the passage of a front are also discussed. 

2. The Basic Methodology 

Because the method has been described In great details elsewhere 
(Ga'i-Chen, 1978), only a brief description will be given here. The 
governing hydrodynamlcal eqs. for conservation of momentum In x-y-z 
directions (z-the vertical; x-y-the horizontal directions) are written as; 


ap/ax » -poDu/Dt + fj = 

F 

(1) 

3p/3y a -p 0 Dv/Dt + f 2 = 

G 

(2) 

pg + ap/az » -poDw/Dt + 

f 3 »« 

(3) 


where p-pressure fluctuations; p-density fluctuations from the basic 
density-poi u, v, w-velocities in the x-y-z directions respectively; 
t-tim?.; fj (1 * 1, 2, 3) - forces other than pressure gradients that can be 
specified in terms of the observed kinematics (e.g., turbulent friction, 
coriolis force); g-gravlty; Du/Dt, Dv/Dt, Dw/Dt - accelerations In the 
x-y-z directions respectively. The right hand side of Eqs. (1), (2), and 
(3) are presumably known or can be calculated from the observed kinematics 
and are therefore denoted symbolically as F, G, and H. 

Inspection of (1) and (2) reveals that the system is overdetermined, 
it will have a solution If and only if 

aF/ay a aG/ax (4) 

If the measurements were error free and the modeling assumptions (e.g., 
the formulations of turbulent frictions) are exact, then one Is assured 
that (4) is satisfied. In practice (4) is not satisfied and thus (1) 
and (2) do not have a solution in the usual sense. The system can 
nevertheless be solved in the least square sense, i.e.. 
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Hid p/ax - F) 2 ♦ (3p/3y - 6) 2 ] dxdy - Min (5) 

This .Is a standard variational problem (Courant and Hilbert, 1953, 

pp. 164-274). The resulting Euler Eq. Is a Poisson Eq. for the pressure 

fluctuations 

3 2 p/3x 2 + 3 2 p/3y 2 » 3F/3x + 36/3y (6) 

with the Neumann boundary conditions 

(ap/3x)n x + (3p/3y)n y - Fn x ♦ Gn y (7) 

Here n x and n y are the direction cosines of the normal to the boundary. 
Eq. (6) with the boundary conditions (7) can be solved at each horizontal 
level where kinematics measurements are available. 

Inspection of (7) reveals that If a solution to (6) exists, It Is 
not uniquely determined. Instead, If p(x, y, z Q ) Is a solution at a 
particular level z * z Q , the p(x, y, z Q ) + c(z Q ) Is also a solution. To 
remove this spurious component one can subtract from the pressure 
its horizontal average. Thus, the method allows the estimation of the 
deviation of the pressure from its horizontal average . The vertical 
profile of the horizontally averaged pressure need to be estimated by 
other means (Gal-Chen, 1978). 

To deduce the density fluctuatlons-p , one substitutes the calculated 
pressure fluctuations In Eq. (3). Because only the deviations of the 
pressure from its horizontal average are given, the same is true for 
the density; consequently , one modifies the calculated H (Eq. 3) by 
subtracting from it its horizontal average. 



3. Possible Approaches to Verification 
■a. Statement of the Problem 

Because the method outlined In the previous section Is Indirect* 

verifications of the derived density and pressure fluctuations are 

critical. Before discussing verification techniques. It is useful to 

outline possible error sources. These could be the following: 

(l.a) The given kinematics, while qualitatively correct Is not 

accurate enough to allow calculations of derivative of velocity product 

2 2 2 2 

such as suw/az, 3vw/az, 3w /3z, v u, v v, v w, etc. These quantities are 
needed in order to estimate F, G, and H (Eqs. (1), (2), and (3)). 

(2. a) The scanning rate is not fast enough to capture the temporal 
evolution. This again, may lead to a significant error in the estimation 
of F, G, and H (Gal-Chen, 1978). 

(3. a) The data is correct, but does not satisfy the particular 
modeling assumptions. For instance, in the case of a three dimensional 
cloud model, the turbulent friction parameterization may be wrong and 
if turbulent friction contribute substantially to the F, G, H estimation 
(Eqs. (1), (2), and (3)) a serious error may result. Another common 
source of error is the use of finite differences to estimate spatial 
and temporal derivatives. At least five points/wavelength are needed 
to get a fair estimate of the first derivative. 

While obviously direct verifications are the most desirable, they 

are very difficult to come by, parti culary, with respect to pressure 

. _ . 

fluctuations. Consequently, we are concerned in this review on three ' 
indirect verifications. They are: numerical simulations; momentum 


checking; and time continuity and physical plausibility. These verifi- 
cation methods are only partially objective and a lot remains to be done 
before one can confidently claim that the technique Is verified. 

m»» \ • 

b. Numerical Simulations as a Guide to Verifications 

To test the viability of the numerical procedure, and also to get 
an order of magnitude estimate of the sensitivity of the method to 
observational errors, numerical models have proven to useful. Gal-Chen 
(1978) and Hane, Wllhelmson, and Gal-Chen (1981) have used wind data 
generated by numerical models as a replacement for observations. From 
this data, they deduced the density and pressure fluctuations and 
compared them to the "exact" density and pressure. The sensitivity of 
the method is assessed by Inserting simulated errors into the "observed" 
kinematics, and examining the resulting "error" In the temperature and 
pressure field. The simulations Indicated that errors of +25 cm/sec for 
boundary layer studies and *50 cm/sec for convective storm studies can 
be tolerated. It was also found that a scanning rate of 2 mln/volume 
for boundary Uyer observations and 4 mln/volume for convective storm 
studies Is acceptable. It must be stressed however, that the numerical 
simulations do not take into account many other sources of systematic 
errors such as pulse volume averaging or incorrect lower/upper boundary 
conditions of the vertical velocity (Carbone et al., 1980). They also 
by definition cannot evaluate errors which are due to Incorrect modeling 
assumptions. 
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To simulate boundary layer observations Gal-Chen has used Deardorff 
(1974a) model. For severe storm simulations, Hane et al. have used the 
. Klemp-WI lhelmson (1978) cloud model. With no errors added the procedure 
recovers exactly the vertical velocities, buoyancy, and pressure. To 
test the sensitivity of the procedure to observational errors we have 
contaminated the data as follows: (a) addition of random noise to the 

"observed" horizontal velocities +15 cm/sec for the boundary layer 
case and +_50 cm/sec for the severe storm case, (b) We have taken 
"observations" which are 2 min apart in the boundary layer case and 4 
min apart in the severe storm case. This simulates the fact that radar 
measurements are non-simultaneous. In both models the actual time step 
(i.e., interval between successive time integration) is the order of 
10 sec. Thus, the models, generate data every 10 sec, but we "observed" 
it every 2 or 4 minutes. The results for the boundary layer case are 
summarized in Fig. 1 and for the severe storm case in Fig. 2. In 
Fig. 1 four curves are plotted as a function of z: 1) <!e'|^ here 8* 

is the potential temperature deviation from the horizontal average, | | 

Is a symbol for absolute value and is a symbol for horizontal average; 
2 ) < C ( 1 6 ' 1 - <|e' |>) 2 > 1/2 i.e.* the standard deviation associated with 

3 ) <Ue'|>= <|e* - e' obs |^> wbere e ' 0 bs the " observed " 
potential temperature deviation, | A8 ' J is the "observed" absolute value 
of the "error"; and 4) <( 1 Ae' | • <| Ae* |>) 2 > 1/2 , the standard deviation , 
of the absolute value of the error. In Fig. 2 we have also plotted 
the temperature statistics as a function of z. Solid line is the 
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retrieved temperature, broken line Is the standard deviation cf the 
temperature. Short dashed lines (near axis) Is the measure of the error 
(average and standard deviation). From these simulations we have concluded 
that our procedure may be a viable one. 
c. Momentum Checking 

It has already been mentioned (section 2) that only under error free 
conditions the deduced horizontal pressure gradients will exactly balance 
the right hand side of the horizontal momentum Eqs. (Eqs. (1), (2)); rather 
they are a least square fit to these equations. The quantity E r defined 
by, 

//[{a P/ax - F) 2 + (9 P /ay - G) 2 ] dxdy _ (8) 

//CF 2 + G 2 ] dxdy " r * 

could be used as a measure of how good the retrieved pressure gradients 
are. When E r = 0 the fit is perfect. If F, G are generated by white 
noise then one can show that E f = 0.5. On the other hand, if F, G are 
perfect (i.e., 3F/3y = 3G/3x) but Eq. (6) is solved by second order 
finite diference approximations rather than exactly, E r = 0.25. Overall 
when E p >_ 0.5 we suspect that the retrieved pressure gradients are 
virtually useless. It is also worth noting that low E r values are 
necessary, hut not sufficient indicators of good quality retrieved 
pressure gradients. For instance one can generate artificially a 
kinematic flow field which is also potential (i.e., v x u = 0) under 
that condition Eq. (4) is automatically satisfied and consequently 
the fit is perfect; but incorrect. 


Gal -Chen (1978) finds typical E r values In the order of .1. These 
values are obtained for simulated boundary layer flows with 100 m 
horizontal resolutions. The flow fields are perturbed by random white 
noise of theorder +15 cm/sec. The scanning rate is of the order of 
2 Min/Vol. Gal -Chen and Kropfll, on the other hand, find typical 
E r values of the order of .35. The wind data that they have used is 
based on dual Poppler radar observations of the planetary boundary layer. 
The spatial resolution is 250 m in the horizontal and 200 m in the 
vertical. The scanning rate is about 90 $ec/Vol. Table 1 sunmarizes 
these results for 5 vertical levels and for selected volume scans. 

Because time tendencies 3/St are needed for the calculations of F and G 
(Eqs. (1) and (2)) one scan is actually a combination of two scans. Thus 
scan 1 is a combi nation of scan 1 and 2, etc. Because there is a gap of 
100 sec between scan 2 and 3, it is not possible to retrieve the pressure 
for scan 2. Taking into account the courser resolutions of the radars 
(compared to the simulated flow) the results are very encouraging. 

Smaller E r values could most probably be obtained for more organized 
convective systems (e.g., severe storms, squall line) provided that one 
can be assured about the quality of the retrieved vertical velocities. 

Clearly, momentum checking is an important simple consistency check 
that should be used prior to any scientific interpretations of the 
retrieved pressure. 



d. Time Continuity and Physical Plausibility as a Guide for 
* Veri fl cation? 

This kind of "verifl cation" is necessarily qualitative and subjective 

and must be pursued with caution and skepticism. In Figs. 3 and 4 we 

have plotted the pressure fluctuation deviations from its horizontal 

averages for two scans which are approximately 200 sec apart. The wind 

data which forms the basis for the pressure retrieval is deduced from a 

dual Doppler* radar measurement taken in September 1978 as part of 

PROJECT PHEONIX (Kropfli and Hilderbrand, 1980). The height where they 

are plotted is 500 m above the ground, roughly half the height of the 

mixed layer for that particular day. The boundary layer at that time 

seemed to be weakly unstable. Because the larger scale features of 

wind field have shown some continuity and persistence from scan to scan, 

it is reasonable to expect the pressure fluctuations to show similar 

behavior. This is indeed the case. The magnitude of the pressure 

-2 -2 

fluctuations is 0.5 - 2 pascals (.5 x 10 - 2 x 10 mb) in agreement with 

values obtained by Deardorff (1974b) for numerically simulated planetary 
boundary layer flows. The corresponding E r values (Eq. 8) are 0.33 
for the first scan (Fig. 3) and 0.26 for the second scan. 

The low E r values, the time continuity of the pressure field, and 

_ o 

the fact that values of the order of 10 mb are expected for a weakly 
unstable planetary boundary layer, gives some confidence in the 
viability of the method. Direct verifications of such low pressure 
fluctuations is impossible with present technology. On the other hand. 


pressure fluctuations play a crucial role in the dynamics of unstable 
boundary layer (see e.g., Zeman, 1981) and great uncertainty exists 
about the proper modeling of the pressure terms in the eqs, for the 
second moments statistics of a turbulent fluid. 

As another illustration of the potential of the technique, 
we present the results of Chong et al. They compute the horizontal 
gradients of various terms In the vertial Eq. of motion. Fig. 5 presents 
these components in a horizontal plane which is approximately 2500 m 
above the ground. The upper part displays the horizontal gradient of the 
vertical acceleration [a(Dw/Dt)/3x, 3(Dw/Dt)/3y)]; the middle part is 
the horizontal gradient of the perturbation potential temperature 
[30/3x, 8&/3y]; the lower part is the horizontal gradient of the 
vertical pressure force [-3/3X (3p/3z), -3/3y (3p/3z)]. For comparison 
purposes these gradients are normalized and presented as °K/km. Two 
distinct relative maxima in the vertical acceleration are apparent, one 
at the left and one at the right. The left updraft has positive 
bouyancy associated with it, and a vertical pressure gradient opposite 
in sign to the bouyancy force. The right updraft has negative bouyancy 
and again a vertical pressure gradient of a opposite sign. In the former 
case, the updraft is associated with thermally unstable air which enters 
the cell from the west, approximately 1500 m above the ground (note the 
display is for air 2500 m above ground). Thus, in that case, the rising 
motion is due to thermal instability with pressure gradient as a retarding 
force. In the latter case, the air enters the cell at an even lower 
altitude (=1000 m) where it is thermally stable . In that case the rising 
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motion is due to pressure gradient forces (probably Induced by surface 
convergence) with bouyancy as a retarding force. The role of the non- 
hydrostatic pressure as revealed by this study, corresponds closely to that 
suggested by conceptual models (Cotton, 1975) and seems highly plausible. 

Gaillard and Glllet also use dual Doppler radars to study convection 
forced by the passage of cold fronts. Their retrieved pressure at. 
the lower levels appears to agree qualitatively with surface measurements. 

It curves out cyclonlcally at upper levels as is to be expected. 

4. Summary and Conclusions 

It appears that the combined use of the equations of fluid dynamics 
and the three dimensional kinematics field derived from multiple Doppler 
radar observations also permits deduction of pressure and density 
fluctuations. In order for the method to work several conditions must 
be satisfied: 

(a) high quality wind data that will allow calculations of spatial 
and temporal first and second derivatives, 

(b) sufficiently fast scanning rate to capture the temporal 
evolutions of the spatial scale resolved by the radars, and 

(c) a good fluid dynamical mode.*l of the phenomenon under study. 

Studies using real data are encouraging. It appears that pressure 

fluctuations with 30% accuracy can be retrieved. Retrieved temperature 
and pressure demonstrate the importance of non-hydrostatic pressure gradients 
in convective dynamics and appear to display continuity and coherence. 

Firm verification of these results, as well as further technique developments 



(particularly In retrievals of vertical velocities) should be pursued 

with' vigor. Almost everything In this vast and difficult problem still 
remains to be done. 
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TABLE CAPTIONS 


- Table 1: Values of E r (Eq. 8) for selected scans and vertical levels. 
It took approximately 100 s to scan a volume. The vertical distance 
between consecutive scans Is 200 m. E f Is a measure of the average 
relative error of the retrieved pressure gradients. 



1 

.31 

.28 

.44 

.37 

2 

.35 

.29 

.35 

.30 

3 

.33 

.26 

.38 

.32 

4 

.36 

.45 

.32 

.32 

5 

.44 

.23 

.37 

.42 
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FIGURE CAPTIONS 

Figure 1. Statistical analysis of the calculated virtual potential 
temperature deviations e' ■ e - <e> for the boundary layer simulation. 

Curve A Is a vertical profile of <|e'|>; curve B 1s^(|e'| - 
i.e., the standard deviation associated with A; curve C Is the average 
"observation error" as' |>; curve D Is the standard deviation associated 
with C. 

Figure 2. Same as Figure 1, but for the severe storm simulations. Solid 
line Is for the retrieval temperature, broken line is the standard 
deviation of the temperature, and short dashed lines (near axis) Is the 
measure of the error. 

Figure 3. Pressure deviations (from its horizontal average) for level 3 
(=<500 m above the ground) and scan 1 (corresponding to the first 100 s 
of the experiment). The horizontal resolution Is 250 m. The values of the 
pressure deviations are given In Pascals (1 mb * 100 Pa). 

Figure 4. Same as Figure 3, but for scan 3. Scan 3 and scan 1 are 
approximately 200 s apart. Note that larger scale features do exhibit 
continuity and persistence from scan to scan. 

Figure 5. Horizontal gradient In vertical acceleration of air motion. 
Horizontal gradient in potential temperature perturbation e‘. Horizontal 
gradient in vertical pressure force. In o der for the thr^e gradients to be 
comparable, they are all normalized and represented in °K km ^ (Source: 

Chong et al., 1980) 
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